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Preface 



A scientist often encounters established ideas that were once the subject 
of debate, sometimes controversy. Often, we use those ideas with no 
knowledge of their historical development, nor of the assumptions on 
which they are based. We rarely stop to ponder the validity of an 
established idea. This is not surprising as this is how we have been 
building our edifice of physical theories, by standing on the shoulders of 
giants, to paraphrase Isaac Newton. 

Yet established ideas and theories need to be challenged and revisited 
when new data or new theories that contradict or shed new light on 
them, become available. We need not be afraid of new information that 
risk overturning accepted ideas. After all, this is how new paradigms 
arise and how progress is achieved. 

The question of whether stars are gaseous or liquid is one debate 
that most scientists are oblivious to. Yet this was a subject of vigorous 
debate in the late 19th and early 20th centuries, with well-known physi- 
cists lining up behind both sides of the question. Larissa Borissova and 
Dmitri Rabounski provide a summary of the history of this debate and 
a personal perspective on how they were pulled into it. 

Recent evidence for liquid stars, in particular the extensive research 
performed by Pierre-Marie Robitaille who has proposed the liquid 
metallic hydrogen model of the Sun* leads us to revisit this question. 
Interestingly enough, stellar plasmas are modelled using Magnetohydro- 
dynamics, i.e. magnetic fluid dynamics, a combination of Maxwell's 
equations of electromagnetism and the Navier-Stokes equations of fluid 
mechanicst Magnetohydrodynamics is also used to model liquid metals. 
This is an indication that the theory of liquid stars is highly plausible 
as an explanation of solar and stellar astrophysical data. 

My interest in this research area stems from the astrophysical re- 
search I performed on stellar atmospheres of Wolf-Rayet stars at the 
University of Ottawa's Department of Physics for my thesis on "Laser 
Action in C IV, N V and O VI Plasmas Cooled by Adiabatic Expansion" . 
Wolf-Rayet stars exhibit mass loss and an expanding stellar atmosphere. 
This results in population inversion of certain atomic transitions due to 

* Robitaille P.-M. A high temperature liquid plasma model of the Sun. Progress 
in Physics, 2007, vol.3, issue 1, 70-81. 

tTajima T. and Shibata K. Plasma Astrophysics. Perseus Publishing, Cam- 
bridge, 2002; Kulsrud R. M. Plasma Physics for Astrophysics. Princeton University 
Press, Princeton, 2005. 
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the rapid cooling of the expanding plasma and the recombination of the 
free electrons into higher excited ionic states, and laser action in the 
corresponding emission lines. This physical mechanism has been pro- 
posed as the explanation for the prominent spectral lines observed in 
the spectra of Wolf-Rayet stars. 

In this book, Larissa Borissova and Dmitri Rabounski provide a 
general relativistic theory of the internal constitution of liquid stars, 
a model that was lacking till now. This they accomplish by using a 
mathematical formalism first introduced by Abraham Zelmanov for cal- 
culating physically observable quantities in a four-dimensional pseudo- 
Ricmannian space, known as the theory of chronometric invariants. 
This mathematical formalism allows to calculate physically observable 
(chronometrically invariant) tensors of any rank, based on operators of 
projection onto the time line and the spatial section of the observer. The 
basic idea is that physically observable quantities obtained by an ob- 
server should be the result of a projection of four-dimensional quantities 
onto the time line and onto the spatial section of the observer. 

This analysis allows them to propose a classification of stars based on 
three main types: regular stars (covering white dwarfs to super- giants), 
of which Wolf-Rayet stars are a subtype, neutron stars and pulsars, and 
collapsars (i.e. black holes). Their theory also provides a model of the 
internal constitution of our solar system. It provides an explanation for 
the presence of the asteroid belt, the general structure of the planets 
inside and outside that orbit, and the net emission of energy by the 
planet Jupiter. 

The ultimate test of any theory of stellar structure is the stellar mass- 
luminosity relation which is the main empirical relation of observational 
astrophysics. Using their theory, the authors can calculate the pressure 
inside stars as a function of radius, including the central pressure. As 
pointed out by the authors, the temperature of the incompressible liquid 
star does not depend on the pressure, only on the source of stellar energy 
(as opposed to a gas, in particular as given by the equation of state of 
an ideal gas). The authors match the calculated energy production of 
the suggested mechanism of thermonuclear fusion of the light atomic 
nuclei in the Hilbert core (the "inner sun" ) of the stars to the empirical 
mass-luminosity relation of observational astrophysics, to determine the 
density of the liquid stellar substance in the Hilbert core. 

Pulsars and neutron stars are found to be stars whose physical ra- 
dius is close to the radius of their Hilbert core. They are modelled by 
introducing an electromagnetic field in the theory due to their rotation 
and gravitation. Electromagnetic radiation is found to be emitted only 
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from the poles of those stars, along the axis of rotation of the stars. 

This book represents a solid contribution to our understanding of 
stellar structure from a general relativistic perspective. It provides a 
general relativistic underpinning to the theory of liquid stars. It raises 
new ideas on the constitution of stars and planetary systems, and pro- 
poses a new approach to stellar structure and evolution which is bound 
to help us better understand stellar astrophysics. 

Ottawa, September 2, 2013 Pierre Millette 

Astrophysics research on stellar at- 
mospheres, Department of Physics, 
University of Ottawa (alumnus) 



Chapter 1 
Problem Statement 



§1.1 Introducing a new theory of the internal constitution of 
stars 

In this book, we suggest a new mathematical theory of the internal 
constitution of stars, and the sources of stellar energy. The theory is 
based on the common consideration of a star and its field according to 
the General Theory of Relativity. 

This is an alternative to the conventional theory of stars which was 
introduced in the 1920's by Arthur Eddington [1] and others in the 
framework of classical mechanics and thermodynamics. 

As is known, the conventional theory resulted in the model of gaseous 
stars: stars are considered as gaseous spheres, consisting of mostly hy- 
drogen and a very inhomogeneous interior so that the hydrogen of the 
extremely hot and dense central region is carried into a process of energy 
release. It assumes, after Hans Bethe [2], that this exothermic process 
is thermonuclear fusion producing helium from hydrogen. Two other 
versions of the gaseous model differ from Eddington's theory in details. 
Edward Milne [3] conjectured two (or more) different states of matter 
inside a star. Nikolai Kozyrev [4] arrived at the peculiar picture of low 
density and temperature inside stars, and a non-nuclear source of stellar 
energy. 

Another theory of the internal constitution of stars was much pop- 
ulated in the 1920-1930's due to James Jeans [5,6]. This is the model 
of liquid stars. The public discussion between Jeans, who defended the 
liquid model, and Eddington, the follower of the gaseous model, was 
fixed in the dozens of short communications they published in the sci- 
entific journals. Indeed, Eddington eventually won. Despite the many 
astrophysical evidences of liquid stars, Jeans' theory did not possess a 
solid mathematical basis. His theory was based on observational anal- 
ysis and arguments rather than a mathematical model. In contrast, 
the theory of gaseous stars was mathematically well-grounded by Ed- 
dington. In particular, the mathematical model of gaseous stars gives 
a theoretical deduction of the mass-luminosity relation, which is one 
of the main relations of observational astrophysics*. This is a "trump 

*The most comprehensive deduction of the mass-luminosity relation in the frame- 
work of the model of gaseous stars is given in Part I of Kozyrev's paper [4] . 
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card": once the gaseous model predicts the mass-luminosity relation, 
the theory is usually claimed to be true in general while all its incon- 
sistencies with observational analysis are merely some "difficulties" to 
be resolved in the future. Thus the model of gaseous stars became the 
conventional model for decades to come, until the present time. 

We now have to make an important note. As is known, the core of 
the mathematical theory of the internal constitution of stars consists of 
two equations: the equation of mechanical equilibrium and the equa- 
tion of heat equilibrium. The mechanical equilibrium means that the 
weight of any unit volume of the stellar matter is put into equilibrium 
with the pressure from within the star. The heat equilibrium (energy 
balance) means that the energy produced within any unit volume of the 
stellar matter is put into equilibrium with the energy flow (radiations, 
convection, or heat conductivity) from within the star onto its surface. 
These two main equations of the theory come from general physics, and 
they do not depend on whether the stars are out of gas, liquid, or some- 
thing else. Only then, by introducing the equation of state of ideal gas 
(and many other particular assumptions) into the main equations, the 
conventional theory yields gaseous stars and other conclusions including 
the mass- luminosity relation. 

Jeans' theory of liquid stars cannot follow this way. The equation 
of state of ideal liquid, provided by classical physics, is so simple that 
it contains not the characteristics of stellar matter which are necessary 
for further deduction by means of the equations of equilibrium. 

Instead of all these considerations of classical mechanics and ther- 
modynamics, we suggest an absolutely different approach to the prob- 
lem. It is based on the simultaneous consideration of a star and its 
field according to the General Theory of Relativity. We consider liquid 
stars: this matches certain new observational evidences for the state of 
condensed matter inside stars; in particular, that the Sun consists of 
high-temperature liquid metallic hydrogen (see [7-10]). 

In the framework of the General Theory of Relativity, the structure, 
matter, and field of such a star are characterized by Schwarzschild's 
metric of a sphere filled with incompressible liquid. The recent theo- 
retical result obtained by one of us [11, 12] showed that, if the Sun is 
represented as a liquid sphere according to the metric, the Sun's field 
has a space breaking (discontinuity) in the asteroid strip (this implies 
that the space breaking impedes the substance to be formed as a planet 
in this orbit). We are therefore sure, hereby, of following the right path. 

We deduce Einstein's field equations in the form that models stars 
as liquid spheres. This is a particular form of field equations, which 
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may or may not satisfy the particular space metric. Therefore, we then 
prove that the obtained particular form of the field equations satisfies 
Schwarzschild's metric of a liquid sphere. 

Then, on the basis of the obtained energy- momentum tensor of a per- 
fect liquid (as contained on the right-hand side of the field equations), 
we deduce the conservation law for the liquid substance of regular stars. 
Solving the equations of energy-momentum conservation, we obtain the 
pressure and density of the liquid substance inside the stars. We then 
obtain the formula for the luminosity of stars in the framework of the 
liquid model. We study how this theoretical formula can be compatible 
with the mass-luminosity relation (which is the main empirical rela- 
tion of observational astrophysics). As a result, we obtain the physical 
characteristics of the mechanism that produces energy inside the stars. 

Concerning the stellar energy mechanism itself, we conclude that it 
is the conversion of substance into radiation on the surface of the tiny 
central "core" inside each star. The core can have a different density 
than that of the other substance of a star (a liquid sphere is homoge- 
neous inside), and is selected by the collapse surface with the radius 
determined according to the star's mass. Despite almost all the mass 
of the star is located outside the core (the core is not a black hole) , the 
force of gravity approaches infinity on the surface of the core due to the 
inner space breaking of the star's field therein. The super-strong force 
of gravity is sufficient for the transfer of the necessary kinetic energy 
to the lightweight atomic nuclei of the stellar substance, so that the 
process of thermonuclear fusion begins. The energy produced by the 
thermonuclear fusion is that very energy which stars radiate, in other 
wordss, the tiny core of each star is its luminous "inner sun" , while the 
produced stellar energy is then transferred to the physical surface of the 
star due to thermal conductivity (which is regular to liquid media). 

Neutron stars and pulsars, being rapidly rotating objects, consist 
a special type of stars. The structure, matter, and field of such stars 
should be described by another metric, which is that of a rotating liquid 
sphere under special physical conditions (which are particular to neutron 
stars and pulsars) . We introduce such a metric. According to the metric, 
the liquid substance of neutron stars and pulsars is in the same state 
as high-density physical vacuum. We then deduce a particular form of 
Einstein's field equations which satisfies the metric. We show that the 
energy-momentum tensor of the obtained field equations satisfies the 
conservation law only in the case where the energy flow from within the 
object is very anisotropic, and is directed toward the north and south 
poles (while the axis of the magnetic field does not coincide with the axis 
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of rotation of the object). This matches the well-known observational 
data about neutron stars and pulsars. 

This is our plan for the upcoming chapters. As a result, we obtain a 
mathematical theory of liquid stars, and of the sources of stellar energy 
according to the General Theory of Relativity. 

Before proceeding with the steps, in the next §1.2, we survey the 
space-time metrics we use in our theory. Then we introduce a new 
classification of stars. This classification is based on the location of the 
space breaking of a star's field with respect to its surface (the space 
breaking is calculated according to the metric and proper parameters 
of the star). 

At the end of this chapter, §1.3 gives a survey of the important 
mathematical apparatus of physically observable quantities in the space- 
time of the General Theory of Relativity we require for our further 
calculations. 

§1.2 Modelling a star in terms of General Relativity 

Stars are spherical bodies filled with substance and light. Their fields 
are spherically symmetric as well. Therefore, once considering a star 
in terms the General Theory of Relativity, the structure, matter, and 
field of such an object should be given by a spherically symmetric space 
(space-time) metric. 

Among the space-time metrics known due to the General Theory of 
Relativity, three primary metrics describe spherically symmetric fields. 
These are Schwarzschild's metric of a mass-point, Schwarzschild's met- 
ric of a sphere filled with incompressible liquid, and de Sitter's metric 
which describes a spherical distribution of physical vacuum (A-field). 
All these three metrics will be used in our consideration of stars. 

1.2.1 The mass-point metric 

ds 2 = (l - ^) c 2 dt 2 - - r 2 {dO 2 + sin 2 0 # 2 ) (1.1) 

^ r 

was introduced in 1916 by Karl Schwarzschild [13]. The metric describes 
the field of a spherically symmetric massive body to so large a distance 
from it that the physical size of the body is neglected (assuming the body 
is a mass-point). The metric is written in the spherical coordinates r, 
<f), 9, whose origin meets the mass-point. Also, herein r g = J^ 1 is the 
Hilbcrt radius of the massive body* while M is the body's mass (which 
is the mass of the field source) . 



*This is not the same as the physical radius of the body. At a distance of the 
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According to the metric (1.1), such a space does not rotate or deform. 
The gravitational inertial force (see §1.4 for detail) in the space can be 
deduced on the basis of the component goo of the fundamental metric 
tensor. As is seen, goo of the mass-point metric (1.1) has the form 

900 = 1-^- (1.2) 
r 



Differentiating the gravitational potential w = c 2 (l — y/goo) with respect 
to x l , we obtain 

dw c 2 dg 00 ^ ^ 



dx l 2\/goo dx l 
We then substitute it into the general formula for gravitational iner- 
tial force (1.42) while taking the absence of rotation of the space into 
account. We obtain the formulae for the covariant and contravariant 
components of the gravitational inertial force 

«~&7Ta' F ' = -&' (L4) 

r 

As is seen from the formulae, the gravitational inertial force in a mass- 
point space is due Newtonian gravitation, and is reciprocal to the square 
of the distance r from the gravitating mass. 

The curvature of a mass-point space is due to the Newtonian field 
of gravitation, produced by the massive body in the origin of the coor- 
dinates. This is not a constant curvature space; its curvature decreases 
with distance from the massive body (the field source). At an infinitely 
large distance from the body the space is flat. 

1.2.2 A space filled with a spherically symmetric homogeneous dis- 
tribution of physical vacuum (the A-field in Einstein's field equations) 
without any island of mass presented therein is described by de Sitter's 
metric 

ds 2 = (l - c 2 dt 2 - - r 2 (d6 2 + sin 2 9 # 2 ) . (1.5) 

The metric was introduced in 1918 by Willem de Sitter [16] as a static 
model of the Universe. It is assumed that A < 10~ 56 in the cosmos, so 



Hilbert radius from the center of gravity of the massive body (r = r g ), gravitational 
collapse occurs: in a rotation-free space, this is a state by which the component 300 
of the fundamental metric tensor g a p is zero (goo = 0). See §5.1 and §5.2 for details. 
The Hilbert radius was introduced due to David Hilbert (1862—1944) who considered 
it in 1917 [15] on the basis of Schwarzschild's mass-point metric. It is also known 
as the Schwarzschild radius, despite the fact that Karl Schwarzschild (1873-1916) 
never considered gravitational collapse in his papers [13, 14] . 
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physical vacuum has a very low density therein. A modern version of 
the static model of the Universe is presented in [17]. 

The fundamental metric tensor, via its components according to 
de Sitter's metric (1.5), manifests that such a space does not deform or 
rotate. Therefore, the gravitational inertial force (1.42) in such a space 
is only due to goo which is 

5oo = 1 - ^ • (1.6) 
Accordingly, after the same algebra as previously, we obtain 

«-£-^T, f-£r. (1.7) 

O 1 _ All o 

3 

This is a non-Newtonian gravitational force, which is proportional to 
distance r: the force (A-force) grows along with the distance at which it 
acts. If A < 0 (the observable density of vacuum is positive), this is an 
attraction force. If A > 0 (the observable density of vacuum is negative), 
this is a repulsion force. See Chapter 5 of our book [18] for details. 

The curvature of a de Sitter space is due to the non-Newtonian 
gravitational field produced by physical vacuum (A-field), which homo- 
geneously fills the space. The curvature is the same everywhere within 
the space. This is a constant curvature space, in other words. 

1.2.3 The metric of a sphere filled with incompressible liquid was 
originally introduced in 1916 by Karl Schwarzschild [14] in a truncated 
form containing substantial limitations. He artificially pre-imposed the 
limitations during the deduction of the metric in order to set the field 
free of breaking*. The true metric of a sphere filled with incompressible 
liquid remained unknown until 2009, when one of us deduced it in the 
most complete form [11,12], which is free of any limitations and thus 
takes space breaking into account. The model of stars as liquid spheres 
plays a key role in our theory. We therefore consider the metric of a 
sphere filled with incompressible liquid in the complete form 



ds 1 




(1.8) 



*Actually, once a limitation is pre-imposed on the metric, the geometry of the 
metric space is artificially truncated. In this sense, the metric Schwarzschild intro- 
duced in 1916 is not the genuine metric of the space of a liquid sphere. 
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as that deduced in the papers [11,12]. Herein, a = const is the physical 
radius of the liquid sphere, while r g — 2G J i is the Hilbert radius calcu- 
lated according to the mass M of the liquid (which is the field source). 
The complete deduction of the metric, containing all the necessary de- 
tails, will be presented in §2.1 of the book wherein we will suggest the 
metric for regular stars. 

The metric (1.8) is written for distances r <a: this is the "internal 
metric" of a sphere filled with incompressible liquid. At the surface of 
the sphere (r = a) the metric coincides with the mass-point metric. Also, 
as was proven in [12] (this deduction will be repeated in §2.1 of the 
book), the "outer metric" of the sphere (r>a) is as well the same as 
the mass-point metric: the external field of a liquid sphere is the same 
as the Newtonian gravitational field of a mass-point. 

As is seen from the metric of a liquid sphere (1.8), such a space does 
not deform or rotate. Therefore, according to the definition of the grav- 
itational inertial force (1.42), the force in such a space is only due to 
5oo- In the metric (1.8), we have 




5oo = 7 h l-^- A /1--5* • (1-9) 



After the same algebra as previously, we obtain 

Fi=-^ r = 7 i=-= 1 (1.10) 



F i = - (?r f . v ~ — °; (i.ii) 

a 3^-v^¥ 

Since r <a inside the sphere, Fi < 0 therein. Hence, this is a force of 
attraction. Its numerical value is proportional to distance r. The force 
is zero at the center of the sphere, and reaches its ultimate-high value 
on the surface. 

It is possible to show that the curvature of such a space, being due 
to the aforementioned field of attraction, increases with distance from 
the center of the liquid sphere onto its surface, in other words, the 
space inside a liquid sphere is not a constant curvature space. (We will 
provide the proof and discuss both the four-dimensional curvature and 
the three-dimensional observable curvature of the space in §2.3.) 
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1.2.4 Wc now suggest a new modelling of stars in terms of the 
General Theory of Relativity. 

Consider stars as spherical bodies consisting of liquid. In the frame- 
work of the liquid model, the internal structure, matter, and field of 
a star is described by the metric of a sphere filled with incompressible 
liquid. This is formula (1.8). As was shown above, the force of gravita- 
tion increases therein proportionally to the distance from the center of 
the star. The external field of the star is described by the mass-point 
metric (1.1). In the external field, the regular Newtonian gravitational 
force acts. The force is reciprocal to the distance from the star. 

The field of a liquid sphere, as such, is not continuous everywhere. 
According to the external metric (1.1) and internal metric (1.8) of a 
liquid sphere, its field has space breaking which appears at two distances 
from the center of the liquid sphere. Due to this fact, we now introduce a 
new classification of stars according to the General Theory of Relativity. 
We hereby explain how to build it. 

The space breaking occurs due to the violation of the signature pre- 
scription conditions of the space metric. It means that the space has 
a singularity in that region (surface or volume) wherein at least one 
of the signature conditions is violated. The signature conditions for a 

sign-alternating diagonal metric (H ) as that of the four-dimensional 

pseudo-Riemannian space (which is the basic space-time of the General 
Theory of Relativity) have the form 



The first three are known as the weak signature conditions. The fourth 
is known as the strong signature condition. If one or all weak signature 
conditions are violated, while the strong signature condition is true, this 
is a removable singularity. If the strong signature condition is violated, 
the space-time has unremovable singularity: in this case the solution is 
regularly failed from consideration, because it "has no physical sense". 
Actually, one could not see the physical meaning therein. However, it 
is very meaningful mathematically. We therefore will take any space 
singularity (space breaking) under consideration. 

Consider now the space of a liquid sphere. According to the exter- 
nal metric (1.1) of the sphere, the first signature condition is violated 



5oo > 0 

5oo 5n < 0 

5oo 5n 522 > 0 

5 = 5oo 5n 522 533 < 0 




(1.12) 



§1.2 Modelling a star in terms of General Relativity 



17 



(.9oo = 0) at the distance r = r g from the center: 



5oo = 1 - - = 0 
r 

5oo5n = -K 0 
5oo5n522 = r 2 > 0 
g = -r 4 sin 2 6 < 0 



(1.13) 



The internal metric (1.8) of the sphere manifests that the second, third, 
and fourth signature conditions are violated at the distance 



(1.14) 



r = r br = \ ~ 




from the center of the sphere, where the aforementioned three functions 
approach infinity*: 

900 = 2(i-a) >0 

5oo5ii -> -oo 



5oo.9ii 522 



oo 



.9 = .900.911522 533 



-oo 



(1.15) 



This means that the field of a liquid sphere has space breaking at two 
distances from the center: 

1. The first space breaking occurs on the surface, spherically covering 
the center of gravity of the liquid sphere at the distance of the Hil- 
bert radius r = r g . This is a surface of gravitational collapse (ac- 
cording to the condition g 0 o = 0 of this space breaking) . In other 
words, while a liquid sphere itself may not be a collapsar, it always 
contains a "core" which is selected from the other liquid substance 
by the surface of gravitational collapse. In the case where the liq- 
uid sphere is a star (as in the said model of liquid stars) , each star 
contains such a core. The core is much smaller than the physical 
radius of regular stars: while the radius of the collapsed core 
(Hilbert radius) of the Sun is r g = 2.9 xlO 5 cm (2.9 km), the phys- 
ical radius of the Sun is 7.0xl0 10 cm (700,000 km). We therefore 
refer to it as the inner space breaking; 



*As is known, a function has a breaking when approaching infinity. 
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2. The second space breaking occurs on the spherical surface cover- 
ing the liquid body at the distance r^r — \/ a3 / r g- The distance is 
much larger than the physical radius of regular stars. Thus this is 
the outer space breaking (contrary to the inner space breaking at 
the Hilbert radius). For example, the second (outer) space break- 
ing of the Sun's field occurs at the distance rb r = 3.4xl0 13 cm = 
= 340,000,000 km = 2.3 AU*. This space breaking is located in- 
side the asteroid strip, close to the orbit of the maximal concentra- 
tion of the asteroids (the asteroid strip is situated from 2.1 AU to 
4.3 AU from the Sun). This implies that the space breaking does 
not permit the substance to be formed into a common physical 
body (such as a planet) in this orbit. 
If the physical radius a of a liquid star is the same as the Hilbert ra- 
dius r g — ™ , it is a gravitational collapsar. In this case (r g = a), the 
internal metric of the liquid sphere (1.8) takes the form 

ds 2 = Ul-^\ c 2 dt 2 - -^-j - r 2 (d6 2 + sm 2 9 # 2 ) . (1.16) 

1 - 1? 

This metric, under the particular condition a 2 = j > 0 (thus A > 0), has 
the same form as de Sitter's metric (1.5), 

ds 2 = (l- c 2 dt 2 - - r 2 {de 2 + sin 2 fl # 2 ) , (1.17) 

V 3 ) i - Alf 

3 

which describes a spherical distribution of physical vacuum (the A- field). 
This means that such a collapsed object, which is a liquid sphere in the 
state of gravitational collapse, consists of liquid whose state is close to 
the state of high-density physical vacuum. 

As a result, the new liquid model allows us to introduce a new clas- 
sification of stars according to the location of the space breaking of a 
star's field with respect to the physical surface of the star: 

Type I: Regular stars including the Sun 

The collapsed core (Hilbert radius r g ) of a regular star is many or- 
ders less than the physical radius a of the star (r g <C a) . The outer 
space breaking r^ T of the star's field is located far away from the 
star, in the cosmos (rj, r ^> a). This, of course, is the list of almost 
all the visible stars: super-giants, the Sun, brown dwarfs, and even 
white dwarfs. We will consider regular stars in Chapter 2; 

*1 AU = 1.49 X 10 13 cm (Astronomical Unit) is the average distance between the 
Sun and the Earth. 



vj Djecr 


A/Toco H/T 

iviass ivi , 


Radius 


Hilbert radius 




Space breaking 


Tbr 


lype 




gram 


a, cm 


r s , cm 


a 


rv, cm 


a 




Red super-giant* 


4.0 xlO 34 


7.0 xlO 13 


5.9 xlO 6 


8.4 xlO" 8 


2.4xl0 17 


3.4 xlO 3 


I 


White super-giant* 


3.4 xlO 34 


4.8 xlO 12 


5.0 xlO 6 


1.0 xlO -6 


4.7xl0 15 


9.8 xlO 2 


I 


Sun 


2.0 xlO 33 


7.0 xlO 10 


2.9 xlO 5 


4.1 xlO -6 


3.4 xlO 13 


4.9 xlO 2 


I 


Jupiter (proto-star) 


1.9 xlO 30 


7.1 xlO 9 


2.8 xlO 2 


4.0 xlO -8 


3.4 xlO 13 


4.8 xlO 3 


I 


White dwarf* 


2.0 xlO 33 


6.4 xlO 8 


3.0 xlO 5 


4.7xl0 -4 


2.9 xlO 10 


0.45 xlO 2 


I 


Red dwarfs 


6.7 xlO 32 


2.3xl0 10 


9.9 xlO 4 


4.3 xl0~ 6 


1.1 xlO 13 


4.8 xlO 2 


I 


Brown dwarfs 


1.5 xlO 32 


7.0 xlO 9 


2.2 xlO 4 


3.1 xlO" 6 


4.0 xlO 14 


5.7 xlO 4 


I 


Wolf-Rayet stars 


1.0 xlO 35 


1.4xl0 12 


1.5xl0 7 


1.1 xlO -5 


4.3 xlO 14 


3.1 xlO 2 


la 


Neutron stars 


2.6 xlO 33 


l.OxlO 6 


3.9 xlO 5 


0.39 


1.6 xlO 6 


1.6 


II 


Pulsar 5 


3.9 xlO 33 


1.6xl0 6 


5.8 xlO 5 


0.36 


2.7xl0 6 


1.7 


II 


Black holes 


various 


various 


various 


1 


1 


1 


III 



*Betelgeuse. tRigel. tSirius B. § Radio-pulsar J1903+0327. 



Table 1.1: Classification of stars according to the General Theory of Relativity. The classification is presented 
with the numerical values of the parameters we calculated for the typical members of the families of stars. 
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Type Ia: Wolf-Rayet stars 

They are almost the same as regular stars, except that the power- 
ful stellar wind consisting of the particles of the stellar substance, 
which are permanently erupted from the star, should be taken 
into account (it is the property characterizing Wolf-Rayet stars). 
Stellar wind will be considered in Chapter 3; 

Type II: Neutron stars and pulsars 

For such a star, the radius of the collapsed core is close to the 
physical radius of the star (r g < a) but does not reach it (other- 
wise the star would be invisible for observation) . The outer space 
breaking rv of the star's field is located in the outer cosmos, and 
is also close to the physical surface of the star but does not reach it 
( r 6r<;a)- Also, stars of this Type II rotate at high speeds which 
are close to relativistic velocities. As a result, the metric and 
energy-momentum tensor of such a star differ from those of regu- 
lar stars. These are neutron stars and pulsars. We will focus on 
these stars in Chapter 4; 

Type III: Black holes 

The Hilbert radius r g (radius of the inner space breaking) and 
the radius of the outer space breaking rf, r of a such an object 
meet each other on its physical surface (r g = rv =a). These are 
gravitational collapsars: the condition of gravitational collapse 
(5oo = 0) occurs in the physical surface of such an object, so all of 
its mass is concentrated within the collapsed surface. Black holes 
will be under focus in Chapter 5 of the book. 

This classification is presented in Table 1.1, with the numerical values of 
the parameters calculated for the typical members of the known families 
of stars. 

The new model of liquid stars according to the General Theory of 
Relativity, surveyed in the classification of stars, will be a subject to 
develop in the upcoming chapters. 

§1.3 physically observable quantities 

Before considering stars in terms of the General Theory of Relativity, we 
shall outline a theory of physically observable quantities in curved four- 
dimensional space (space-time) . A comprehensive expisition of the said 
physically observable quantities has already been given in the respective 
chapters of our books [18, 19]. We now give the necessary theoretical 
basics of the theory of physically observable quantities according to [19], 
with some amendments which are required for the current study. 
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In order to build a descriptive picture of any physical theory, we 
need to express the results through real physical quantities, which can be 
measured in experiments (physically observable quantities). In the Gen- 
eral Theory of Relativity, this problem is not a trivial one at all, because 
we are looking at objects in a four-dimensional space-time, and so we 
have to determine which components of the associated four-dimensional 
tensor quantities are truly physically observable. 

Here is the problem in a nutshell. All equations in the General 
Theory of Relativity are cast in generally covariant form, which does 
not depend on our choice of the frame of reference. The equations, as 
well as the variables they contain, are four-dimensional. Thus, we ask, 
which of those four-dimensional variables are truly observable in real 
physical experiments, i.e. which components are true physically observ- 
able quantities? Intuitively we might, at first glance, easily assume that 
the three-dimensional components of a four-dimensional tensor consti- 
tute a physically observable quantity. Yet, at the same time, we cannot 
be absolutely sure that what we simply observe are truly the three- 
dimensional components per se, if not more complicated variables which 
depend on other factors, e.g. on the properties of the physical standards 
of the space of reference. 

As is known, a four-dimensional vector (a lst-rank tensor) has as few 
as 4 components (1 time component and 3 spatial components). A 2nd- 
rank tensor, e.g. a rotation or deformation tensor, has 16 components: 1 
time component, 9 spatial components, and 6 mixed (time-space) com- 
ponents. Now, are the mixed components truly physically observable 
quantities? Tensors of higher ranks have even more components; for 
instance the Riemann-Christoffel curvature tensor has 256 components, 
so the problem of the heuristic recognition of genuine physically ob- 
servable components becomes far more complicated. Besides, there is 
an obstacle related to the recognition of the observable components of 
covariant tensors (in which indices occupy the lower position) and of 
mixed-type tensors, which have both lower and upper indices. 

We see that the recognition of physically observable quantities in the 
General Theory of Relativity is not a trivial problem. Ideally we would 
like to have a mathematical technique to calculate physically observable 
quantities for tensors of any given ranks unambiguously. 

Numerous attempts to develop such a mathematical method were 
made in the 1930's by some of the most outstanding researchers of that 
time. The goal was nearly attained by Landau and Lifshitz in their fa- 
mous The Classical Theory of Fields [20] , first published in Russian in 
1939. Aside for discussing the problem of physically observable quanti- 
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ties itself, in §84 of their book, they introduced the interval of physically 
observable time along with the phsyically observable three-dimensional 
interval, which depend on the physical properties (physical standards) 
of the space of reference of an observer. But all such attempts made in 
the 1930's were very limited to just solving certain particular problems. 
None of them led to a versatile mathematical apparatus. 

A most complete mathematical apparatus for calculating physically 
observable quantities in a four-dimensional pseudo-Riemannian space 
was first introduced by Abraham Zelmanov and is known as the theory of 
chronometric invariants, or the chronometrically invariant formalism. 
It was first presented in 1944 in his Ph.D. thesis [21] — then in his 
condensed papers of 1956-1957 [22,23]. 

The essence of Zelmanov's mathematical apparatus of physically ob- 
servable quantities (chronometric invariants), designed especially for 
the four-dimensional, curved, non-uniform pseudo-Riemannian space 
(space-time), is as follows. 

At any point of the space-time we can place a three-dimensional 
spatial section x° = ct = const (three-dimensional space) orthogonal to 
a given time line x l = const. If a spatial section is everywhere orthogonal 
to the time lines, which pierce it at each point, such a space is referred 
to as holonomic. Otherwise, if the spatial section is non-orthogonal 
everywhere to the aforementioned time lines, the space is referred to as 
non-holonomic. 

Possible frames of reference of a real observer include a coordinate 
net spanned over a real physical body (the reference body of the ob- 
server, which is located near him) and a real clock located at each point 
of the coordinate net. Both the coordinate net and clock represent a set 
of real references to which the observer refers his observations. There- 
fore, physically observable quantities registered by an observer should 
be the result of a projection of four-dimensional quantities onto the time 
line and onto the spatial section of the observer. 

The operator of projection onto the time line of an observer is the 
world- vector of four-dimensional velocity 



of his reference body with respect to him. This world- vector is tangential 
to the world-line of the observer at each point of his world-trajectory, 
so this is a unit-length vector 




(1.18) 



dx a dx 13 



g a/ 3 dx a dx@ 
ds~ 2 



= +1. 



(1.19) 



9aj3 



ds ds 
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The operator of projection onto the spatial section of the observer (his 
local three-dimensional space) is determined as a four-dimensional sym- 
metric tensor h a p, which is 

h a f} = - 9afS + b a bfs 

fret? = _ g aP + 5 a fe /3 

hZ = -g£ + b a bf i 

The world-vector b a and the world-tensor h a p are orthogonal to each 
other. Mathematically this means that their common contraction is zero 
(h a p b a = 0, h a ^b a = 0, b a = 0, bP a b a = 0). So, the main properties of 
the operators of projection onto the time line and the spatial section of 
the observer are commonly expressed, obviously, as follows: 

b a b a = +l, h^b a = 0. (1.21) 

If the observer rests with respect to his reference object (such a case 
is known as the accompanying frame of reference), then b l — 0 in his 
reference frame. The coordinate nets of the same spatial section are 
connected to each other through the transformations 

x° = x° (x°, x\ x 2 , x 3 ) 

x 1 = x 1 (a; 1 , x 2 , x 3 ) , = 0 

where the third equation displays the fact that the spatial coordinates in 
the tilde-marked net are independent of the time of the non-tilded net, 
which is equivalent to a coordinate net where the lines of time are fixed 
x l = const at any point. The transformation of the spatial coordinates 
is nothing but a transition from one coordinate net to another within 
the same spatial section. The transformation of time means changing 
the whole set of clocks, so this is a transition to another spatial section 
(another three-dimensional space of reference). In practice this means 
replacement of one reference body with all of its physical references with 
another reference body that has its own physical references. But when 
using different references, the observer will obtain different results (other 
observable quantities). Therefore, the physically observable projections 
in an accompanying frame of reference should be invariant with respect 
to the transformation of time, which implies invariance with respect 
to the transformations (1.22). In other words, such quantities should 
possess the property of chronometric invariance. 



(1.20) 
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We therefore refer to the physically observable quantities determined 
in an accompanying frame of reference as chronometrically invariant 
quantities, or chronometric invariants in short. 

The tensor h a p of projection, being considered in the space of a 
frame of reference accompanying an observer, possesses all properties 
attributed to the fundamental metric tensor, namely 



hfh 



bib" 




(1.23) 



where 5f is the unit three-dimensional tensor*. Therefore, in the accom- 
panying frame of reference the three-dimensional tensor hit can lift or 
lower indices in chronometrically invariant quantities. 

So in the accompanying frame of reference the main properties of 
the operators of projection are 



b a b a = +l, 



hlb a =0. 



hfhl = 5* 



(1.24) 



Calculate the components of the operators of projection in the accom- 
panying frame of reference. The component b° comes from the obvious 
condition b a b a =g a pb a bP = 1, which in the accompanying frame of ref- 
erence {b % =0) is b a b a =goob°b° = 1. This component, in common with 
the remaining components of b a , is 



b Q = 



1 



bo = go a b a 



'9oo ■ 



b l = 0 
h = g ia b° 



9zo 



(1.25) 



while the components of h a p are 
^00 = 0, h 00 - ~ 00 

h 0 i = 0 , 
h i0 = 0 , 

, . 90i90k ui u 

h ik = - 9ik H , /i 



h 0i = _ r 0i 



h io = _ r i0 



■9 
9° 

9* 



1 

.9oo 



5oo 



hi 


= 0 


hi 


= % = 




_ gio 




9oo 


K 


= -9l 



> . (1.26) 



*This tensor 5^ is the three-dimensional part of the four-dimensional unit ten- 
sor <5? , which can be used to replace indices in four-dimensional quantities. 
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Zelmanov created a comprehensive mathematical method for the cal- 
culation of the chronometrically invariant (physically observable) pro- 
jections of any generally covariant (four-dimensional) tensor quantity, 
and set it forth as a theorem (we refer to it as Zelmanov 's theorem) : 

Zelmanov 's theorem 

Assume that Qqo"o are the components of the four-dimensional 
tensor Qqq 'q of the r-th rank, in which all upper indices are not 
zero, while all m lower indices are zero. Then, the quantities 

T* k - p = ( 5 oor ? Q:;c' (1.27) 

constitute a chronometrically invariant three-dimensional contra- 
variant tensor of (r — m)-th rank. Hence the tensor T lk --- p is a re- 
sult of m-fold projection onto the time line by the indices a, f3 . . . a 
and onto the spatial section by r — m the indices ji, v . . . p of the 
initial tensor Q^" p a - 

According to the theorem, the chronometrically invariant (physically 
observable) projections of a four-dimensional vector Q a are 

b a Q a = -^, h l a Q a = Q\ (1.28) 



while the chr.inv. -projections of a symmetric tensor of the 2nd rank Q a/3 
are the following quantities: 

6 a 6 /3 Qa^ = — , h ia bPQ a p = -Q^=, h l a h\Q a P = Q lk . (1.29) 
5oo v.9oo 

The chr.inv. -projections of a four-dimensional coordinate interval 
dx a are the interval of the physically observable time 

dr = y^dt+ 9 °' dx\ (1.30) 

and the interval of the observable coordinates dx 1 which are the same as 
the spatial coordinates. The physically observable velocity of a particle 
is the three-dimensional chr.inv. -vector 

v J = ^-, ViV* = h ik v\ k = v 2 , (1.31) 

CIT 

which at isotropic trajectories becomes the three-dimensional chr.inv.- 
vector of the physically observable velocity of light 
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Chronometrically projecting the covariant or contravariant funda- 
mental metric tensor onto the spatial section of an accompanying frame 
of reference (b l = 0) 

hf hi g al 3 = g lk - hb k = - h ik 

h i h p 9 aP = 9 lk -b l b k = g lk = - h lk 

we obtain that the chr.inv. -quantity 

hik = - gik + bib k (1-34) 

is the chr.inv. -metric tensor (the observable metric tensor), using which 
we can lift and lower indices of any three-dimensional chr.inv. -tensorial 
object in the accompanying frame of reference. The contravariant and 
mixed components of the observable metric tensor are, obviously, 

h ik = -g ik , hi = -gi=S i k . (1.35) 

Expressing g a p through the definition of h a p — — g a p + b a bp, we ob- 
tain the formula for the four-dimensional interval 

ds 2 = b a bp dx a dx p - h aP dx a dxP, (1.36) 

expressed through the operators of projection b a and h a p. In this for- 
mula b a dx a — cdr, so the first term is b a bp dx a dx@ = c 2 dr 2 . The second 
term h a ^dx a dx^ =da 2 in the accompanying frame of reference is the 
square of the observable three-dimensional interval* 

do 2 = h lk dx l dx k . (1.37) 

Thus, the four-dimensional interval, represented through the physically 
observable quantities, is 

ds 2 = c 2 dr 2 -do 2 . (1.38) 

The main physically observable properties attributed to the accom- 
panying space of reference were deduced by Zclmanov in the framework 
of the theory, in particular — proceeding from the property of non- 
commutativity (non-zero difference between the mixed 2nd derivatives 
with respect to the coordinates) 

dx l dt dtdx* c 2 1 dt' 1 ' ' 



* This is due to the fact that h a p in the accompanying frame of reference possesses 
all properties of the fundamental metric tensor g a p. 
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*d 2 



*d 2 



— A — 

c z at 



dx 1 dx k dx k dx 1 
where the chr.inv. -operators of differentiation of Zelmanov are 
*d 1 d *d d g 0i d 



dt 



/ffoo dt 



dx 1 



dx i g 00 dx° ' 



(1.40) 



(1.41) 



The first two physically observable properties are characterized by 
the following three-dimensional chr.inv. -quantities: the vector of the 
gravitational inertial force Fi and the antisymmetric tensor of the an- 
gular velocities of rotation of the space of reference A ik which are 



Ft = 



1 



dw 



/goo \dx l 



dt 



1 



Aik= 2\d^~dt^ + ~^ FiVk ~ FkVi) - 



2 c 2 



(1.42) 



(1.43) 



Here w and Vi characterize the body of reference and the reference space. 
These are the gravitational potential 



w = c 2 (1 - y/itoS) , 
and the linear velocity of rotation of the space 



1-4 
c z 



'goo 



Vi 



goi 
s/goo 



cg 0i y/go~o~ 



Vi = h ik v , 



v 2 = v k v k = h lk v l v k 



(1.44) 



(1.45) 



We note that w and Vi do not possess the property of chronometric 
invariance, despite Vi = h ik v k can be obtained as for a chr.inv. -quantity, 
through lowering the index by the chr.inv. -metric tensor hi k . 

Zelmanov also found that the chr.inv. -quantities Fi and Ai k are 
linked to each other by two identities (Zelmanov's identities) 



*3A ik 1 f*dF k *dFA 
dt 2\ dx 1 dx k ) 



*dA kn 
dx 1 



*dA m . 
dx k 



dt 

*dA lk 
dx m 



0. 



(1.46) 



(F l A km + F k A mi + F m A lk ) = 0 . (1.47) 



In the framework of quasi-Newtonian approximation, i.e. in a weak 
gravitational field at velocities much lower than the velocity of light and 
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in the absence of rotation of space, Fi becomes a regular non-relativistic 
gravitational force Fj = ^ . 

Zclmanov also proved the following theorem setting up the condition 
of holonomity of space: 

Zelmanov's theorem on holonomity of space 

Identical equality of the tensor to zero in a four-dimensional 
region of space (space-time) is the necessary and sufficient condi- 
tion for the spatial sections to be everywhere orthogonal to the 
time lines in this region. 

in other words, the necessary and sufficient condition of holonomity of 
a space should be achieved by equating to zero the tensor . Naturally, 
if the spatial sections are everywhere orthogonal to the time lines (in 
such a case the space is holonomic), the quantities goi are zero. Since 
goi = 0, we have Vi = 0 and A^ — 0. Therefore, we will also refer to the 
tensor An- as the space non-holonomity tensor. 

If the conditions Fi — 0 and A^ = 0 are met in common somewhere 
in the space, the conditions goo = 1 and goi = 0 are as well true therein. 
In such a region, according to (1.30), dr — dt: the difference between 
the coordinate time t and the physically observable time r disappears 
in the absence of gravitational helds and rotation of the space. In other 
words, according to the theory of chronometric invariants, the difference 
between the coordinate time t and the physically observable time r 
originates in both gravitation and rotation attributed to the space of 
reference of the observer (the local space of the Earth in the case of an 
Earth-bound observer), or in each of these physical factors separately. 

On the other hand, it is doubtful to find such a region of the Universe 
wherein gravitational fields or rotation of the background space would 
be absent in clear. Therefore, in practice the physically observable time 
r and the coordinate time t differ from each other. This means that the 
real space of our Universe is non-holonomic, and is filled with a grav- 
itational field, while a holonomic space free of gravitation can be only 
a local approximation to it. 

The condition of holonomity of a space (space-time) is linked directly 
to the problem of integrability of time in it. The formula for the interval 
of the physically observable time (1.30) has no integrating multiplier. 
In other words, this formula cannot be reduced to the form 

dr = Adt, (1.48) 

where the multiplier A depends on only t and x % : in a non-holonomic 
space the formula (1.30) has non-zero second term, depending on the 
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coordinate interval dx l and g^. In a holonomic space A^ = 0, so goi = 0. 
In such a case, the second term of (1.30) is zero, while the first term is 
the elementary interval of time dt with an integrating multiplier 



Hence time is integrable in a holonomic space (An- = 0), while it cannot 
be integrated in the case where the space is non- holonomic (An- ^ 0). In 
the case where time is integrable (a holonomic space), we can synchro- 
nize the clocks in two distantly located points of the space by moving 
a control clock along the path between these two points. In the case 
where time cannot be integrated (a non- holonomic space) , synchroniza- 
tion of clocks in two distant points is impossible in principle: the larger 
the distance between these two points is, the more the deviation of time 
on these clocks is. 

The space of our planet, the Earth, is non-holonomic due to the 
daily rotation of it around the Earth's axis. Hence two clocks located at 
different points of the surface of the Earth should manifest a deviation 
between the intervals of time registered on each of them. The larger 
the distance between these clocks is, the larger the deviation of the 
physically observable time (expected to be registered on them) is. This 
effect was surely verified by the well-known Hafele-Keating experiment 
[24-27] concerned with displacing standard atomic clocks by an airplane 
around the terrestrial globe, where rotation of the Earth's space sensibly 
changed the measured time. During a flight along the Earth's rotation, 
the observer's space on board of the airplane had more rotation than 
the space of the observer who stayed fixed on the ground. During a 
flight against the Earth's rotation it was vice versa. An atomic clock 
on board of the airplane showed a significant deviation of the observed 
time depending on the velocity of rotation of space. 

Because synchronization of clocks at different locations on the sur- 
face of the Earth is a highly important problem in marine navigation 
and also aviation, in an early time de-synchronization corrections were 
introduced as tables of the empirically obtained corrections which take 
the Earth's rotation into account. Now, thanks to the theory of chrono- 
metric invariants, we know the origin of these corrections, and are able 
to calculate them on the basis of the General Theory of Relativity. 




(1.49) 



so we are allowed to write the integral 




(1.50) 
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In addition to gravitation and rotation, the reference body can de- 
form, changing its coordinate nets with time. This fact should also be 
taken into account in measurements. This can be done by introducing 
into the equations the three-dimensional symmetric chr.inv. -tensor of 
the rate of deformation of the space of reference 



= 1 *dh lk 
lk ~ 2 dt 



D ik 
D = 



1 *dh ik 
2~ 



dt 



> . 



h* k D,. 



dt ' 



h = det || h ik | 



(1.51) 



The regular Christoffel symbols of the 2nd rank (r" ) and the 1st 
rank (r^ ;CT ) 

'dg^a dg u<T dg^ 



dxf* dx° 



are linked to the respective chr.inv. -Christoffel symbols 



■ h l 



dx k 



*8h 



km 



*dh 



)k 



dxi 



dx r ' 



(1.52) 



(1.53) 



which are determined similarly to T™ v . The only difference is that here, 
instead of the fundamental metric tensor g Q( 3, the chr.inv. -metric ten- 
sor hik is used. The Christoffel symbols characterize the property of 
inhomogeneity of space. 

The components of the regular Christoffel symbols are linked to the 
other chr.inv. -chractersitics of the accompanying space of reference by 
the following relations: 

; 9ok Tqo " 



Dl + At = 



>9oo 



Ok 



goo 



(1.54) 



ia fe/3 pm 
9 9 1 a/3 



2 yk 
c 1 00 

5oo 
= h iq h ks &Z 



(1.55) 



(1.56) 



Wc now express the chr.inv.-Christoffcl symbols through the chr.inv. - 
propcrties of the accompanying space of reference. Expressing the com- 
ponents g a/3 and the first derivatives from g a p through F i: An-, D^, w, 
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and Vi , after some algebra we obtain 
1 / w \ dw 



Toco 



c 2 / St ' 
1 /. w\ 2 1 <9w 

r 1 (i w ) dw 

c 2 V c 2 J dx % ' 

1 / w \ / 1 \ 1 <9w 

1/, w \l~^ l ( dvi dvi\ 1 _ x" 

c i 1 - ?) h - 2 W + & ) + 2? ^ + ^. 



1 / <9v 7 9«i 
1 



2 c 2 



Vk(F t Vj +FjVi) 



■F k ViVj , 



r° - - — - 

1 oo - c3 



1 dw 



1 c 2 

■ft-M 1 -?)'* 4 - 



r° - — 

0i ~ c 2 



c 2 

rg, = i(i-5)(i3f + ^ + l« i ^) ) 



r u . = - 



x Vj (£>? + A;?)+ Vi (£>; + ^) + i,; iVj F 



(1.57) 
(1.58) 
(1.59) 
(1.60) 

(1.61) 

(1.62) 
(1.63) 
(1.64) 
(1.65) 

(1.66) 



(1.67) 
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Zelmanov also deduced formulae for the chr.inv. -projections of the 
Ricmann-Christoffel curvature tensor. He followed the same procedure 
by which the Ricmann-Christoffel tensor was built proceeding from the 
non-commutativity of the second derivatives of an arbitrary vector Q a 
taken in the given space. Taking the non-commutativity of the second 
chr.inv. -derivatives of an arbitrary vector 

*Vi*V fc Q, -*V k *ViQi = 2 ^ k * d g t l ~ "ih-Qj > (!- 69 ) 
where the chr.inv. -covariant differential from the vector is 

*V k Q l dx k = dQ l + A l kl Q k dx\ 
he obtained the chr.inv. -tensor 

*8A 3 *3A j 

TT - j _ ^ ^ fci I A m A 3 A m A 3 (1 71 

lkl ~~dx~ k IW a km kl lm ' ^ > 

which is like Schouten's tensor in the theory of non-holonomic mani- 
folds [28]. The tensor H' lki , differs from the Riemann-Christoffel tensor 
Rp^s". due to the presence of space rotation A ik in the formula (1.69). 
Nevertheless, its generalization gives the chr.inv. -tensor 

Cikij = ^ {Hikij — Hjkn + Hkiji — Hujk) , (1-72) 

which possesses all the algebraic properties of the Riemann-Christoffel 
tensor in this three-dimensional space. Therefore, Zelmanov called Cikij 
the chr.inv. -curvature tensor, which actually is the tensor of the phys- 
ically observable curvature of the three-dimensional spatial section of 
the observer. Its contraction step-by-step 

Ckj — C ki j. — h im Ckimj i C = Cj — h lj Cij (1-73) 

gives the chr.inv. -scalar C which is the observable three-dimensional 
curvature of this space. 

The tensor Hikij is connected with the curvature tensor Cikij by 

Hikij — Cikij + (ZAkiDji + AijDki + AjkDu + 

+ AkiDij+A H Djk). (1.74) 

The contracted tensors Hik = H lk l and Cik = C lk l are connected as 

Hk =Cik + ^ (AkjDf + A X jD{ + AkiD) . (1.75) 



(1.70) 
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In a particular case where the space does not rotate, FLxuj and Cikij are 
the same. This is as well true for Hi k and Cik- In this particular case, 
the tensor Cik = h^Cukj has the form 



The Riemann-Christoffel tensor R a /3-yS, being a two-pair symmetric 
tensor (its paired indices are non-symmetric inside each pair, while the 
pairs are symmetric with respect to each other), has three chr.inv.- 
projections according to formula (1.29) of the chronometrically invariant 
formalism. They are as follows: 

xlk = _ c 2 ^ yijk = _ c ^1 1 z ijkl = c 2 R^ kl . (1.77) 

Substituting the necessary components of the Riemann-Christoffel ten- 
sor R a f3yS into the formulae for its chr.inv. -projections (1-77), and by 
lowering indices, Zelmanov obtained the formulae 

X ij = ^-(Dl + At)(D jl +A jl ) + 

+ (*V i F j +*V j F i )-±F i F j , (1.78) 

Y ijk = *V t (D jk + A jk ) - *Vj (D ik + A lk ) + -| A tJ F k , (1.79) 

Zikij — DikDij — DuDk-j + A ik Aij — 

- A u A kj + 2A i:j A kl - c 2 C lk i 3 , (1.80) 

where Y(ij k ) = Yijk + Yjki + Ykij = 0 just like in the Riemann-Christoffel 
tensor. Contraction of the observable spatial projection Z ik ij step-by- 
step as Zu — h kj Zikij and Z = h ll Zu gives 

Z u = D ik D k - D U D + A lk A k + 2A lk A k l - c 2 C a , (1.81) 

Z = K l Z a = D lk D' k -D 2 - A lk A lk - c 2 C . (1.82) 

At the end of our survey of the chronometrically invariant formalism, 
consider Einstein's field equations* 

R a p ~ ^ gapR = - xT al } + \g a p . (1.83) 



*The left-hand side of the field equations (1.83) is often referred to as the Einstein 
tensor G a p = R a p ~ \ 9a/3R> i n notation G a p = — xT a ^ + \g a p- 
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The field equations, except for the fundamental metric tensor g a p, in- 
clude: R a fj — R' aa p. is Ricci's tensor (the second-rank symmetric tensor 
coming from contraction of the Riemann-ChristofFel curvature tensor), 
R = g a/3 R a fj is the curvature scalar, h= = 18.6 xl0~ 28 cm/gram 
is Einstein's constant of gravitation, G — 6.672 xlCP 8 cm 3 /gramxsec 2 
is Gauss' constant of gravitation, T a p is the energy- momentum tensor 
of the matter distributed in the space, and A [cm -2 ] that describes 
physical vacuum (see §5.2 of the book [18]). 

Landau and Lifshitz [20] use x = instead of x = as used by 
Zelmanov. To understand the reason, set x= as in our study, and 
consider the chr.inv. -projections of the energy-momentum tensor 

p=^, J l = ^=, U lk = c 2 T k , (1.84) 

which come with formula (1.29) as the projections of any second-rank 
symmetric tensor. They have the following physical meaning: p is the 
observable density of mass, J % is the observable density of momentum, 
and U lk is the observable stress-tensor. Ricci's tensor has dimension 
cm~ 2 . This means that the scalar chr.inv. -projection of the field equa- 
tions, = — - T °° + A, and the quantity * To ° = 87rC *p nave same 
' goo goo ' ^ ■' goo c 2 

dimension which is cm~ 2 . Hence, the energy-momentum tensor T a p 
has the same dimension as mass density (gram/cm 3 ). Therefore, once 
we would use xr = on the right-hand side of the field equations, we 
would use not the energy-momentum tensor T a p but rather c 2 T a p. 

The chr.inv. -projections of Einstein's equations (1.83) are calculated 
as those of a second-rank tensor (1.29). They have the form (we refer 
to them as the chr.inv. -Einstein equations) 

*8D / 1 \ 

— + DuD* + AjiA lj + (JV, - - 2 Fjj F 3 = 

= -^(pc 2 + U)+Xc 2 , (1.85) 

*Vj (h ij D - D 13 - A 13 ) + -| FjA ij = xJ\ (1.86) 

^ - (Dij + A l3 ) (Dl + A£) + DD lk - DijDl + 
+ SAijAi + \ (*^F k + *V fe F,) - \ FF k - c 2 C lk = 

= I (pc 2 h ik + 2U ik - Uh ik ) + \c 2 h lk , (1.87) 
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where U = h lk Uik is the trace of the stress-tensor ?7 ife . 

Also, the energy-momentum tensor T a p of distributed matter should 
satisfy the law of conservation which is 

V a T aa =0. (1.88) 

The chr.inv. -projections of the conservation law are calculated as those 
of a first-rank tensor (1.28). We refer to them as the conservation law 
equations. The equations have the form 

^ + Dp + 1 Dij + * V, J* - 1 Fi X = 0 , 

*r) J k ~ 
— + DJ k + 2 (D k + A k ) J 1 + * V, U lk - P F k = 0 , 

where the chr.inv. -operator *Vi = *Vi — -^Fi is created on the basis of 
the chr.inv. -differential operator (see Notations) . 

Given these definitions, we can find how any geometric object of a 
given four-dimensional pseudo-Riemannian space (space-time) is consti- 
tuted from the viewpoint of any observer whose location is this space. 
For instance, having any equation obtained in the generally covariant 
tensor analysis, we can calculate the chr.inv. -projections of it onto the 
time line and onto the spatial section of any particular body of refer- 
ence, then formulate the respective chr.inv. -projections in terms of the 
physically observable properties of the reference space. This way we will 
arrive at fully qualified equations containing only quantities measurable 
in practice. 

Thus, we now have all the necessary mathematical "equipment" re- 
quired for our further development of the mathematical theory of the 
internal constitution of stars, and of the sources of stellar energy, ac- 
cording to the General Theory of Relativity. 



(1.89) 
(1.90) 



Chapter 2 
Regular Stars and the Sun 



§2.1 Introducing the space metric of a regular star. Ein- 
stein's field equations in the form satisfying the metric 

In this chapter, we introduce the new mathematical theory of liquid 
stars being applied to regular stars. This means Type I of stars in 
terms of the new classification we have just introduced according to the 
General Theory of Relativity (see §1.2, and Table 1.1 therein). It covers 
the widest variety of stars, which includes super-giants, sun-like stars 
(including the Sun), dwarfs, and, white dwarfs*. 

The structure, matter, and field of a liquid star are characterized 
by Schwarzschild's metric of a sphere filled with incompressible liquid. 
The metric was originally introduced in 1916 by Karl Schwarzschild [14]. 
He, however, introduced it in a truncated form containing substantial 
limitations: he artificially pre-imposed these limitations during the de- 
duction in order to set the field free of breaking, thus resulting in the 
geometry of the metric space artificially truncated. In other words, the 
metric introduced by Karl Schwarzschild is not quite the genuine metric 
of the space of a liquid sphere. The true metric of a sphere filled with 
incompressible liquid, which is free of the said limitations, thus takes 
space breaking into account, as deduced in 2009 by one of us [11, 12]. 
We now repeat the deduction here, according to the most detailed ex- 
planation [12], along with some recent amendments and comments. 

Consider an empty space that houses a spherical island which is a 
liquid. The structure, matter, and field of such an massive island should 
be characterized by a space metric which possesses spherical symmetry. 
As is known, all spherically symmetric metrics have the following general 
form 

ds 2 = e v c 2 dt 2 - e x dr 2 - r 2 {d6 2 + sin 2 #d0 2 ) , (2.1) 

where e v and e A are functions of r and t. 

The matter and field of the spherical island (which is a liquid) should 
satisfy Einstein's field equations (1.83), which in the case under consid- 



*White dwarfs are considered separately in the framework of Eddington's theory 
of gaseous stars. 
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eration have the A-field neglected, i.e. 




(2.2) 



where R a p is Ricci's curvature tensor, R is the curvature scalar, x = 
= = 18.6 xlO -28 cm/gram is Einstein's constant of gravitation, and 
T Q/ 3 is the energy- momentum tensor of the distributed matter (liquid). 
The energy-momentum tensor (i.e. the distributed matter) should sat- 
isfy the conservation law 



where V CT is the four-dimensional symbol of covariant differentiation (see 
Notations). 

Einstein's field equations connect the components of the fundamen- 
tal metric tensor, the space curvature, and distributed matter according 
to Riemannian geometry In other words, the invariant square form of 
Ricmannian metric, ds 2 = g a p dx a dx l3 = inv, in common with Einstein's 
field equations characterize Riemannian spaces (the spaces whose geom- 
etry is Riemannian) . Concerning the General Theory of Relativity, this 
means as follows. Let us have a Riemannian space having a metric ds 2 , 
and suggest that matter be distributed in it (thus we suggest a particular 
formula for the energy- momentum tensor T a p). Then, the components 
of the fundamental metric tensor g a p (known from the formula of the 
metric ds 2 ) and the components of the suggested energy- momentum ten- 
sor, being commonly substituted into (respectively) the left-hand side 
and the right-hand side of Einstein's field equations should transform 
the equations into identities. 

This is the way how, on the basis of the general formulae of a spheri- 
cally symmetric metric (2.1), to deduce the metric of a sphere filled with 
liquid. We take the energy-momentum tensor of a perfect liquid, then 
substitute its components into the right-hand side of the field equations. 
Then we take the components of the fundamental metric tensor from 
the general spherically symmetric metric (2.1) in their general (non- 
particular) form containing the coefficients e v and e A . We substitute 
the components into the left-hand side of the field equations. Then we 
look which form of the coefficients e" and e A makes the left-hand side of 
the field equations the same as the right-hand side (thus transforming 
the field equations into identities). Finally, we substitute the obtained 
particular formulae for the coefficients e v and e A back into the general 
formula of spherically symmetric metrics. Voila! The metric of a sphere 
filled with perfect liquid has been obtained. 



v„r ff = 



o, 



(2.3) 
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One might as well just ask, why did Schwarzschild himself not do 
just that? Instead, why did he follow another complicated way, full 
of assumptions and suppositions? Well. . . Let us come back to our 
deduction. 

As is known, the energy-momentum tensor of a perfect liquid (which 
is incompressible and non- viscous) has the form 

T a ? = (po + Jr) U a U? - J g*?, (2.4) 

where p = pa = const is the density of the liquid (which is constant), p 
is the pressure, while 

dr a 

U a = — , U a U a = 1 (2.5) 
as 

is the four-dimensional velocity of the liquid flow with respect to the ob- 
server (his reference space coincides with the space of the liquid sphere, 
with the origin of the coordinates located at the center). 

Hence forth we express the field equations in component notations 
with the physically observable properties of the space selected. 

We see that 

9oo = e v , g 0 i = 0 

9n = - e A > .922 = - r 2 , g 33 = -r 2 sin 2 9 

in the metric of spherically symmetric spaces (2.1). According to the 
chronometrically invariant formalism (see §1.3), the gravitational po- 
tential in such a space has the following formulation 

w = c 2 (l-e*). (2.7) 

Because goi = 0 in the metric, the space does not rotate. Therefore, the 
linear velocity of the rotation is Vi = 0 as well. Hence the chr.inv. -tensor 
of the angular velocity of space rotation is zero 

^=K^-S) + ^ <F *- F " )=o ' < 2 ' 8 > 

while the chr.inv. -vector of gravitational inertial force has the form 

where the prime denotes differentiation along the radial coordinate r. 
With these, the chr.inv. -metric tensor hik of the space has the non- 
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zero components 

hu = e\ h 22 = r 2 , h 33 = r 2 sin 2 *? , (2.10) 

/i n -e- A , h 22 = ±, h 33 = —^-, (2.11) 
r 1 r 2 sin 2 9 

h = dct || h ik || = eV sin 2 9 . (2.12) 

Thus, the chr.inv. -tensor of the rate of space deformation, Dn~, has only 
the following non-zero components: 

Oii^e^, r> u = ^e- A -5, £=^e"S, (2.13) 

where the upper dot denotes differentiation along the time coordinate 
t. The chr.inv. -Christoffcl symbols, which characterize space inhomo- 
geneity, are calculated according to their definition given in §1.3 with 
the components of the chr.inv.-metric tensor hik- After some algebra, 



we obtain formulae for the non-zero components of A ijim 
y 

An,i = ye A , A 22 ,i=-r, A 33 ,i = - r sin 2 0 , (2.14) 

Ai2, 2 = r, A 33 , 2 = - r 2 sin 0 cos 9 , (2.15) 

Ai 3 ^3 = r sin 2 0 , A 23 , 3 = r 2 sin 0 cos 9 , (2.16) 
and formulae for the non-zero components of A^- 

Ah = y, A 22 = -re- A , A33 = - r sin 2 0e~ A , (2.17) 

A 2 2 = ^, A 2 3 = - sinflcosfl, (2.18) 

A? 3 = -, A^cotfl. (2.19) 



As was shown in §1.3, in a rotation-free space the second-rank chr.inv. - 
curvature tensor Cik = h^Cukj, which is the physically observable cur- 
vature tensor, has the form (1.76). A rotation-free space is the case 
under consideration. After some algebra, we obtain the non-zero com- 
ponents of Cik for the spherically symmetric metric (2.1). They are 
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Let us calculate the chr.inv. -projections of the energy-momentum 
tensor of a perfect liquid (2.4) according to the associated projections* 
which are (1.84). With 6* = 0 and 6°=— L= (1.25) which characterize 
an accompanying frame of reference (in the case under consideration, 
the observer accompanies the liquid sphere), we obtain 

p=^=p 0 , r = -^==0, U lk =c 2 T tk =ph tk . (2.21) 

900 y^goo 

whcrcfrom we also have, for U = h tk Ui k , 

U = 3p. (2.22) 

The obtained condition J 1 — 0 means that the liquid is free of flow, while 
U lk =ph lk means that the observer's reference frame accompanies the 
liquid medium. 

The chr.inv. -Einstein equations (1.85-1.87) in a rotation-free space 
now take the simplified form 

+ D Jt D« + ('V, - 1 F,) F' = -% (p 0 c 2 + U) , (2.23) 

*Vj {h ij D — D ij ) = 0 , (2.24) 
*dD ik „ „, „ „, 1 



- DijDl + DD lk - D tJ Di + - (*V t F k + *V fe F,) - 



Qj. "«7"fc 1 ""IK "»J"fc 1 2 

- ^ FiF k - c 2 C lk = £ ( Po c 2 h lk + 2U lk - Uh ik ) , (2.25) 

where *Vi is the symbol of chr.inv. -differentiation (see Notations). The 
chr.inv. -equations of the conservation law (1.89-1.90) also simplify to 

Dpo + ^DijU** =0, (2.26) 

*V, J7 ife - po^ fe = 0 , (2.27) 

where *Vi = *Vi — ^F l (see Notations). 

Substitute, into the chr.inv. -Einstein equations (2.23-2.25), the chr. 
inv. -characteristics of the space we have obtained above for the spheri- 
cally symmetric metric (2.1) and also the obtained chr. inv. -components 



*They are the observable density of mass, p, the observable density of momen- 
tum, J*, and the observable stress-tensor U tk of the liquid. 
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of the energy-momentum tensor of a perfect liquid. After some algebra, 
we obtain the chr.inv. -Einstein equations (2.23-2.25) in component no- 
tation (the third tensorial equation splits into three, where the second 
and third equations remain the same) 



A - 



Az> A 2 



2v' V 2 
V ^- L - 

r 2 

= - x(p 0 c 2 + 3p) e A , 



AV 
2 



— e 

r 



= 0, 



\v A 2 



— c 



A2 



AV K) 
2 2 



+ 



2c 2 A' 



c 2 (A' - i/') 



2 c 2 
^(1 



= ^(p 0 c 2 -p) e A , 
- A ) =^(p 0 c 2 -p). 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



The second equation manifests that A = 0 in this case. This means that 
the inner space of the liquid sphere does not deform: with A = 0, we 
have Dn = 0, D n = 0, and D = 0 according to (2.13). Taking this cir- 
cumstance into account, as well as the stationarity of A, we reduce the 
field equations (2.28-2.31) to the final form 



c 2 e~ A 



AV 2v' V 2 

I 1- 

2 r 



2c 2 A' 



AV 

— + 2 



2 

O') 2 



x(pac 2 +3p) e A , 
= x(p 0 c 2 -p)e x , 
° 2 { ~ - V>) e- X + 2 4(l- e~ X ) = » (p Q c 2 p) . 



(2.32) 
(2.33) 
(2.34) 



To solve the field equations (2.32-2.34), we need a formula for the 
pressure p. To find the formula, we now deal with the conservation equa- 
tions (2.26-2.27). However, due to the absence of space deformation in 
the case under consideration (Dm =0), the chr.inv. -scalar conservation 
equation (2.26) vanishes. Only the chr.inv. -vectorial conservation equa- 
tion (2.27) remains. It takes the form 



*Vi(pti k )- (po + ^)F k =0. 



(2.35) 
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Since * V, h lk — 0 is true always (as well as V a g aa = 0 for the fundamen- 
tal metric tensor), the remaining conservation equation (2.35) reads 

hlk -B~( Pa + i) Fk= °- (2 - 36) 

Because = j n a rotation- free space, this formula reduces to a 
non-trivial equation which has the form 

p'e- x + (p 0 c 2 +p)^ e- x = 0, (2.37) 

where p' = i>'= e A ^0. Dividing both parts of (2.37) by e" A , we 
obtain 

which is a plain differential equation with separable variables. It easily 
integrates as 

p 0 c 2 +p = Be~%, B = const. (2.39) 

Thus, we obtain the pressure p as a function of v, which is 

p = Be~% - p 0 c 2 . (2.40) 

In looking for an r-dependent function p(r), we integrate the field 
equations (2.32-2.34). Summarizing (2.32) and (2.33), we find 

c2 ( A, + ^=^-f. (2.41) 
r 

Express v' herefrom, then substitute the result into (2.34). We obtain 

V A ' + ( eA - !) - ^e- A -f = x (p 0 c 2 - p) e\ (2.42) 

Substituting p from (2.40) into (2.42), we obtain the following differen- 
tial equation with respect to A: 

A' + e — - - xp 0 re x = 0. (2.43) 
r 

We introduce a new variable y = e x . Thus A' = Substituting into 
this equation y and y', we obtain the Bernoulli equation (see Kamke [29] , 
Part III, Chapter I, §1.34) 

y' + f(r)y 2 + g(r)y = 0, (2.44) 
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where 

f(r) = --xpor, g{r) = --. (2.45) 



It has the following solution: 

1 -*<'>/ {2M) 



V 

where 

E(r) = e Sa{r)dr . (2.47) 
Integrating (2.47), we obtain E{r) which is 

E(r) = e^ 1 ^ = e ln * = ^ , L = const > 0, (2.48) 

thus we obtain - = e~ x which is 

v 

e- x = - ( L(l-„ Por )dr = l-?W- + Q, Q = const. (2.49) 
r J L \r J 3 r 

To find Q, we re-write equation (2.42) as 

^(^-^)+^ = xPo- (2-50) 

This equation has a singularity at the point r = 0, where the numerical 
value of the right-hand side term of the equation (the density of the 
liquid) grows up to infinity by r — > 0, i.e. at the center of the sphere. 
This is a contradiction with the initially assumed condition p 0 = const, 
which is specific to incompressible liquids. As a matter of fact, this 
contradiction should not be in the theory. We remove this contradiction 
(and the singularity) by re-writting (2.50) in the form 



e- x (1 - r\') = — (re- x ) = 1 - x Po r 2 . (2.51) 



d_ 

dr 

After integration, we obtain 

re -* = r _ + A , A = const. (2.52) 

Because A = 0 at the central point r — 0, it should be zero at any other 
point as well. Dividing this equation by r ^ 0, we obtain 

e"* = 1 - ^ • (2.53) 
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Comparing this solution with the value e _A obtained earlier (2.49), 
we see that they meet each other if Q = 0. Besides, we should suggest 
that e A ° = 1 at the central point r = 0, consequently Ao = 0. 

Thus we have obtained the components h 11 =e~ x and hn~e x of 
the chr.inv. -metric tensor hik in the form expressed through the radial 
coordinate r, i.e. 

ft 11 =e"* = 1-^, h 11= e x = (2.54) 

3 -, xpor 2 



3 

Hence forth, we should introduce a boundary condition on the sur- 
face of the sphere. We have r = a on the surface, where a is the radius 
of the sphere. Thus 

= 1 - ^ . (2.55) 

On the other hand, the solution of this function is also the mass-point 
solution in emptiness. Hence, we have 

(2.56) 

where M is the mass of the sphere. Comparing both these formulae of 
e~ Aa , and taking into account that Einstein's constant of gravitation is 
x=^f-, we find 

M = = P0 V, (2.57) 

where V = ^p- is the volume of the sphere. Hence, we have obtained 
the regular relation between the mass and the volume of a homogeneous 
sphere. 

Our next step is the look for the solution e~ A outside the sphere, 
where r > a. Since outside the sphere the density of matter (the liquid) 
is po = 0, we obtain, after integration of (2.51), 

re~ x = [ T dr- (\ Po r 2 dr = r - . (2.58) 

Jo Jo 3 

We obtain, from this formula, that 

e- A = 1 - . (2.59) 

Taking (2.55) and (2.56) into account, we arrive at the mass-point so- 
lution in emptiness 

e- = l-H^. (2.60) 
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To obtain v we use equation (2.41). Substituting 

2xpgr 

X'= 2 



1 - 



(2.61) 



and the obtained formula of e A into (2.41), we obtain, after transforma- 
tions, 

2xp 0 r 2 

(2.62) 



V + 



>cB re 2 



0. 



1 



xppi"' 



3 "3 

We introduce a new variable = y. Thus, v' = — -%-. Substituting 
these into (2.62), we obtain the Bernoulli equation 



, xB ry 2 



y' + 



3 y 



2 c 2 



l - 



1 - 



where 



f(r) 



xB 



xpgr 



2c 2 



1 - 



xpgr A 



g(r) 



1 - 



(2.63) 



(2.64) 



Thus, we have the integral 

/. xpgr 
g{r)dr = - / 3 - - 



ln./V, 



xp 0 r z 



where 



E(r) = N i 



N = const, (2.65) 



(2.66) 



In the region where the signature condition hu = e x > 0 is satisfied, 
we have 



i-^->o, 



therefore we use the modulus of the function here. 
Next, we look for - = which is 



(2.67) 
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We obtain, after integration, 



We now find the constants B and K . To find B, we re- write the formula 
of p by the condition that p = 0 on the surface of the sphere (r — a). 
Thus, we obtain 

B = p 0 c 2 e^, (2.70) 

where is the value of the function on the surface. As a result, 
we have 



2 



e f = l^i e ^ __ + ^. . (2.71) 



To find if, we take the value of ea on the surface (r = a) 



xp 0 

We obtain, from this formula, that 



= — + JfJl . (2.72) 



T7~~ ' . • (2-73) 
'l 



The quantity means the numerical value of by r = a (i.e. on 
the surface of the sphere). Therefore, we can apply it to the mass-point 
solution in emptiness at r~a, i.e. 



Taking the formulae of e~r , (2.55) and (2.56), into account, we obtain 



ef = i | 8 



1 - 


xpor 2 \ 
3 




1 


>fp 0 a 2 
3 / 





This formula on the surface (r = a) meets the mass-point solution in 
emptiness: = jl-*^ = yjl -^- 
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Thus the metric of the space of a sphere filled with perfect liquid is, 
since the formulae of v and A have already been obtained, as follows: 



ds 2 = M 3 Ji-^-Ji-^l Jde- 



. 2 



4 \ V 3 V 3 

-r 2 (^ 2 +sin 2 0# 2 ). (2.76) 



dr 2 



1 



3 

Taking (2.55) and (2.56) into account, we re-write the formula (2.76) as 



4 



1 - 



r 2 (d0 2 • sin 2 0# 2 ). (2.77) 



2GMr 2 



Finally, since 2< ^ 2 M = r g is the Hilbert radius calculated according to 
the mass M of the liquid sphere, while taking the obtained formula of 
e~ into account, we re- write the metric in the final form 



, 2 

1 



<l»- = - I 3 , / f - — ■ 1 ^-fj c 2 dt 2 

dr 2 



-r 2 (d6 2 + sin 2 <9# 2 ) . (2.78) 
1 ^ 

This is the final formula for the "inner" metric of the space of a sphere 
filled with perfect liquid. As is seen, the "inner" metric completely 
coincides with the mass-point metric in emptiness on the surface of the 
liquid sphere (r = a). 

Hence forth, we obtain the space metric outside the liquid sphere 
(r > a). We have already obtained the "external" solution for e~ A (2.59), 
which coincides with the "external" mass-point solution for this function 
(2.60). Outside the sphere, B = 0 (2.39). Hence, (2.41) takes the form 

A' + z/ = 0, (2.79) 

where, according to (2.60), 

c 2 r 2 i _ 2GM y ' 
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Substituting (2.80) into (2.79) then integrating the resulting equation, 
we obtain 



v = In ( 1 - —j— )+P, P = const, (2.81) 

thus 

c 2 r J 



e" = P(l-H^y (2.82) 



Since this function is also 



- = l-^, (2.83) 



on the surface (r = a) of the liquid sphere, we obtain P=l. Having 
the obtained formulae for c v (2.83) and e A (2.60) substituted into the 
spherically symmetric metric (2.1), we obtain that the "outer" space of 
a sphere filled with perfect liquid is described by the mass-point metric 
in emptiness (1.1), which is 



ds 2 



(l _ ^ c 2 dt 2 - - r 2 {d9 2 + sm 2 6 d<j> 2 ) . (2.84) 



§2.2 The outer space breaking of the Sun's field matches 
the asteroid strip 

Herein we suggest a new model of the Solar System according to the 
General Theory of Relativity. Namely — the Sun and the planets will 
be considered as liquid spheres according to the metric of a liquid sphere 
(2.78) we have obtained in the foregoing. The metric was also shown 
in formula (1.8), in §1.2 wherein we surveyed the problem statement of 
the modelling of a star in terms of the General Theory of Relativity. As 
was also proven in the previous §2.1, the outer space of a liquid sphere 
is described by the mass-point metric in emptiness (1.1). 

Note that herein we do not discuss whether the internal planets can 
be represented as liquid spheres or not. Astrophysicists and geologists 
may simply appeal to the magma, because it is in the state of liquid 
stone. However, the jovian planets (Jupiter, Saturn, Uranus, and Nep- 
tune), according to their density and other parameters, can surely be 
considered as stars. Herein, we only limit ourselves to the theoretical 
modelling of the Sun and the planets, without an analysis of their origin 
or other astrophysical factors. In detail, we focus on the location of the 
"inner" and "outer" space breaking of their fields: the space breaking 
of the field within and outside the physical body (liquid sphere). Then 
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we compare the obtained result with the observed distribution of the 
planets within the Solar System. 

Our approach to the Solar System is simple. As is known, given 
a four-dimensional Riemannian space with a sign-alternating diagonal 

metric (H ), the breaking occurs in that region (point or surface) of 

the space wherein at least one of the four signature conditions 



.9oo > 0 

.900.911 < 0 

5oo 3nff22 > 0 

.9 = 3oo 3n 322 333 < 0 



(2.85) 



is violated. The space (space-time) of the General Theory of Relativ- 
ity is one of this type of space. We therefore consider the signature 
conditions in the space within and outside the liquid Sun. 

2.2.1 In the "inner" space metric of a liquid sphere (2.78), while 
taking into account that 

W = 2GM = r, 
3r c 2 r r 

therein* the fundamental metric tensor has the following non-zero com- 
ponents: 

3oo = t 



(2.87) 



3n = — = r , (2-88) 

1 a3 1 3 

322 = ~r 2 , (2.89) 
333 = -r 2 sin 2 fl. (2.90) 

We obtain, from those components, that at a distance from the center 
of the sphere which is 

,3 . .. 

(2.91) 





*See formulae (2.59) and (2.60) in §2.1. 
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the second, third, and fourth signature conditions are violated* 

» = ?(l-^)>0 

goo 9n -> -oo 
.9oo 9i 1.922 -> oo 
.9 = 9oo 9U 922 933 -> -oo 

This means that the field of the liquid spherical body has space breaking 
on the spherical surface covering the body at the distance r^ r = y/a 3 /r g 
from its center. 

The Hilbcrt radius r 9 = 2< ^ 2 M (the radius of gravitational collapse) 
calculated for regular physical bodies is many orders less than their 
physical sizes. Hence, a^>r g for a regular spherical liquid body (thus 
the body is not a collapsar). Therefore, — \/a 3 /r g ^> a: the spherical 
surface of the space breaking of the field is located far away from the 
physical surface of the liquid body (the field source) , and hence far away 
from the inner field. In other words, the inner field and liquid substance 
of the body produce breaking in the outer space of the body. 

What does the outer space breaking of the field mean from the phys- 
ical viewpoint? Has this space breaking a real action on a physical body 
appearing in it, or is it only a mathematical fiction? As will be shown 
in the next §2.3, the space (space-time) of a liquid sphere possesses 
space breaking in its four-dimensional curvature tensor Rap-yS by the 
condition r = rv. Namely, — the component i?oioi (2.113), which is 
the four-dimensional curvature of the space in the (r-i)-direction 0101, 
possesses breaking at the distance r = rt, r from the center of the liquid 
sphere (the curvature function becomes infinite, i?oioi ^ oo, on the sur- 
face of the radius r — r^r)- Because the four-curvature determines the 
gravitational field which fills the space (and vice versa), the breaking at 
r = rv implies breaking in the gravitational field of the liquid sphere. 

This is the physical sense of the outer space breaking of the field of 
a liquid sphere. 

2.2.2 The outer field of a liquid sphere is due to the same liquid 
substance, which fills the sphere and produces the field within the sphere 
itself (its inner field). According to the "outer" space metric (2.84), 
we see that the fundamental metric tensor of the outer space has the 



(2.92) 



* Namely — these three functions approach infinity. As is known, a function has 
such breaking when approaching infinity. 
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following non-zero components 



9oo = 1 - , (2.93) 
5ii = — , (2-94) 

r 

922 = ~r 2 , (2.95) 
g 33 = -r 2 sin 2 8. (2.96) 

We see that at a distance of 

2GM 

r = r g = — (2.97) 

from the center of the body, the first signature condition (g 00 > 0) is 
violated 

5oo = 1 - - = 0 



r 

5oo.9ii = -1 < 0 
5oo5n522 = r 2 > 0 
g = -r 4 sin 2 6 < 0 



(2.98) 



In other words, the outer field of a liquid sphere produces space breaking 
deep within the sphere itself, in its inner space close to the center. For 
example, the calculated Hilbcrt radius r g = 2 Jf is only 2.9 km for the 
Sun, while for the Earth it is nothing but only 0.88 cm. 

2.2.3 Concerning regular stars, and the Sun in particular, the 
aforementioned findings imply the following (as per our new model of 
liquid stars according to the General Theory of Relativity) : 

1. At the center of each star, a small core exists. The core is sep- 
arated from the other mass of the star by the said inner space 
breaking in the star's field, at the distance of the Hilbert radius 
r g from the center. The inner space breaking means, physically, 
that the liquid substance of the star has a singularity on the sur- 
face of the Hilbert radius r g from the center, thus the small core 
is separated from the major mass (the physical sense of the phe- 
nomenon will be more clear from the example of the outer space 
breaking in the field of the Sun); 

2. The field of each star has an outer space breaking surrounding the 
star by a spherical surface. This "bubble" has a very large radius 
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of = y 7 a 3 /r g , which is many orders larger than the physical 
radius a of the star. Physically, the outer space breaking impedes 
the near substance such as small stones or dust orbiting the star 
to be formed as a planet in the orbit of the radius rv ■ 

Let us now calculate the radius of the outer space breaking of the 
Sun's field by formula (2.91), which is r& r = ^a 3 /r g = '*cp 0 . Substi- 
tute the Sun's density po = lAl gram/cm 3 , or the mass M = 2.0xl0 33 
gram and the radius a = 6.95><10 10 cm. We obtain 

r br = 3.4 x 10 13 cm = 340, 000, 000 km = 2.3 AU, (2.99) 

where 1 AU = 1.49 xlO 13 cm (Astronomical Unit) is the average distance 
between the Sun and the Earth. We obtain that the spherical surface 
(bubble) of the outer space breaking of the Sun's field is located within 
the asteroid strip, very close to the orbit of the maximal concentration 
of asteroids (as is known, the asteroid strip is located, approximately, 
2.1 to 4.3 AU from the Sun). 

This truly amazing finding brings us to a conclusion that the internal 
constitution of the Solar System can be calculated according to the 
liquid model. Namely, — we consider the Sun and the planets as liquid 
spheres, then we calculate the outer space breaking rv in the field of 
each of the cosmic bodies. The results of the calculation are collected 
altogether in Table 2.1. 

These results associated with the planets and the Sun, according to 
Table 2.1, lead to the next conclusions: 

1) The outer space breaking of the Sun's field is located within the 
asteroid strip, near the maximal concentration of asteroids; 

2) The internal planets of the Solar System (Mars, the Earth, Venus, 
and Mercury) are located within the "bubble" of the outer space 
breaking of the Sun's field; 

3) The "bubbles" of the outer space breaking of the field of each of 
the internal planets are as well located within the "bubble" of the 
outer space breaking of the Sun's field; 

4) The outer space breaking of the fields of Mars and the Earth 
reaches the asteroid strip; 

5) The outer space breaking of Mars' field is located at 2.9 AU from 
the Sun. It is within the asteroid strip near the orbit of Phaeton, 
the hypothetical planet which was once orbiting the Sun according 
to the Titius-Bode law at r — 2.8 AU, and whose distraction in 
the ancient time gave birth to the asteroid strip; 



Object 


Mass M, 


Density p 0 , 


Radius 


Hilbert radius 


Orbit, 


Space breaking 


Location of r br 




gram 


gram / cm 


a 1 cm 


r g , cm 


ATT 


ATT 

r br , AU 


from the Sun, AU 


Sun 


1.98 xlO 33 


1.41 


6.95 xlO 10 


2.9 xlO 5 


— 


2.3 


2.3 


Internal planets 














Mercury 


2.21 xlO 26 


4.10 


2.36 x 10 s 


0.03 


0.39 


1.3 


-0.9-1.7 


Venus 


4.93 xlO 27 


0.10 


D.19 x 1U 


U./3 


U. 12, 


1 o 

1.2 


—0.0- 1.9 


Earth 


5.97 xlO 27 


5.52 


6.38 x 10 s 


0.88 


1.00 


1.1 


-0.1-2.1 


Mars 


6.45 xlO 26 


3.80 


3.44 x 10 s 


0.10 


1.52 


1.4 


0.1-2.9 


Asteroid strip 










2.5* 






Jovian planets 














Jupiter 


1.90 xlO 30 


1.38 


7.11 xlO 9 


280 


5.20 


2.3 


2.9-7.5 


Saturn 


5.68 xlO 29 


0.72 


6.00 xlO 9 


84 


9.54 


3.2 


6.3-12.7 


Uranus 


8.72 xlO 28 


1.30 


2.55 xlO 9 


13 


19.2 


2.4 


16.8-21.6 


Neptune 


1.03 xlO 29 


1.20 


2.74 xlO 9 


15 


30.1 


2.4 


27.7-32.5 


Pluto 


1.31 xlO 25 


2.00 


1.20 x 10 s 


0.002 


39.5 


1.9 


37.6-41.4 


Kuiper belt 










30-100 







*The maximal concentration of the asteroids of the asteroid strip is registered at ~2.5AU from the Sun, while the asteroid 
strip continues from 2.1 to 4.3 AU (approximately). 



Table 2.1: The internal constitution of the Solar System according to the General Theory of Relativity. 
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6) The "bubble" of the outer space breaking of Jupiter's field meets, 
from its inner side, that of Mars at r = 2.9 AU from the Sun (this 
is the case of the "parade of the planets"). It is very near 2.8 AU, 
which is the theoretical orbit of Phaeton according to the Titius- 
Bode law; 

7) The "bubbles" of the outer space breaking of the field of the other 
jovian planets (Saturn, Uranus, and Neptune) are located within 
the inner boundary of the Kuiper belt (the strip of the aphelia of 
the comets orbiting the Sun); 

8) The outer space breaking of Neptune's field meets, from the outer 
side of the "bubble" , the inner boundary of the Kuiper belt; 

9) The "bubble" of the outer space breaking of the field of Pluto is 
completely located within the Kuiper belt. 

The fact that the outer space breaking of the Sun's field is located 
within the asteroid strip, near the maximal concentration of asteroids, 
allows us to say: yes, the space breaking considered in this study has 
a real physical meaning. Probably, the Sun's space breaking impedes 
asteroids to be joined into a common physical body (one refers to it 
as Phaeton). Alternatively, if Phaeton was an already existing planet 
orbiting the Sun near the "space breaking orbit" in the past, the force 
of gravitation of another massive cosmic body, emerging near the Solar 
System in the ancient ages (for example, another star passing near it), 
has displaced Phaeton to the "space breaking orbit" near it, thus leading 
to the distraction of Phaeton's body. 

Thus the internal constitution of the Solar System is formed by the 
geometric structure of the Sun's field according to Riemannian geometry 
as manifest within the laws of the General Theory of Relativity. 

§2.3 The geometric sense of the outer space breaking 

Consider the properties of the curvature of the space of a liquid sphere. 
First, let us calculate the components of the chr.inv. -curvature tensor 
Cikij , which is the physically observable curvature tensor of the space. 

In a rotation- free space (Aik = 0), which is the space of a liquid sphere 
under consideration, Cikij — Hikij according to the definition of the ten- 
sor Hikij (1-74). Therefore, we calculate Cikij = Hikij — hjmH^^ 1 by the 
formula of Hfe™ (1.71), wherein we substitute the respective chr.inv.- 
Christoffel symbols A* fc (2.17-2.19) already obtained for the metric of 
a liquid sphere (2.78). After some algebra, we obtain that the chr.inv.- 
curvature tensor Cikij in the space of a liquid sphere has the following 
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non-zero components: 



C1212 — #1212 — ^— , (2.100) 

3 1 _ xpor 2 
1 3 

r - TT - ^Po r 2 sin 2 (9 

W313 — il 1313 — -, [Z.WL) 

3 1 _ xppH 
1 3 

C 232 3 - #2323 = - ^f- r 4 sin 2 9 . (2.102) 

We see that, in the space of a liquid sphere, the non-zero components 
of the observable space curvature tensor Cikij satisfy the condition 

Cikij = {hkihij — huhkj) , (2.103) 

where the negative constant — is the observable three-dimensional 
curvature of the space in the respective two-dimensional direction. This 
means that the three-dimensional space of a non-rotating liquid sphere 
has a constant negative curvature. Calculating the observable curvature 
scalar C = h lk dk, where the non-zero components of dk are 

Cn = -^-^, (2-104) 

3 

C33 2*cp 0 r 2 

22 = ~ 2~~n = o ' (2.105) 

sin 0 3 



we obtain 



C = -2xp 0 = const < 0. (2.106) 



Hence, according to (2.103), the chr.inv. -curvature tensor Cikij is ex- 
pressed through the observable curvature scalar C as 

C 

Cikij = -q (hkihij — huhkj) ■ (2.107) 

Thus, the observable three-dimensional space of a non-rotating liquid 
sphere is a constant negative curvature space. Therefore, the curvature 
radius 5ft of the three-dimensional space is imaginary. It is formulated 
through the observable curvature scalar C by the relation 

C = -2xp 0 = ^, (2.108) 
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thus we obtain, finally, 

$t=—!— . (2.109) 

Let us calculate the components of the full Riemann-Christoffel cur- 
vature tensor 

jj _ 1 / dg aS dgp 7 dg fjS dg aj \ 

aMS ~ 2 {dxPxi dx a dx s dx a dx-< dx^dx 5 J 

+ g aT {^as,tr^pi,T - r/3i )(7 r a7)T ) . (2.110) 

According to the metric of a liquid sphere (2.78), we have g^ = — and 
Fik,j = — Aikj ■ Thus, calculating the non-zero components of T a p t §, 



3^/1-^-0" 



x/A)'' " ' " " ' 



roi.o = -roo.i = — p; / 1 ( 2 - m ) 



3 

xp 0 r ! 



n,i 



(1 - ^) 

and substituting these into (2.110), we obtain 



(2.112) 



-Roioi = - -p5 , i (2.113) 

v 1 - ^ 

-R1212 = —7-^ = — C1212 , (2.114) 



3 


1 >cp 0 r 2 
1 3 


_ xpo 


r 2 sin 2 0 


3 


1 3 


_ xpo 
3 


r 4 sin 2 6 = — 



i?i3i 3 ^^^r = -C 131 3, (2.115) 



i?2 323 = r 4 sin 2 0 = - C 232 3 ■ (2.116) 

We see that the component i?oioi> determining the four-dimensional 
curvature in the (r-i)-direction 0101, does not satisfy the condition of 
four-dimensional constant curvature spaces which is 

Ra/3-fS = Q (gp 7 9aS - 3/3<55a 7 ) , Q = COHSt. (2.117) 

Therefore, we arrive at the following conclusion: 
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• The four- dimensional space of a non-rotating liquid sphere is not 
a constant curvature space. This is in contrast to the observable 
three-dimensional space of the liquid sphere which, as was proven 
above, is a constant negative curvature space. 
We see also, from the formulae for C1212 (2.100) and C1313 (2.101), 
that the three-dimensional observable curvature Cikij possesses space 
breaking 



In other words, the three-dimensional chr.inv. -curvature Cikij and the 
four-dimensional Riemannian curvature Rap^s have space breaking by 
the condition r — r^r- Concerning the model of liquid stars, this means: 
• In the field of each star, the three-dimensional observable space 
curvature Cikij and the four-dimensional Riemannian curvature 
Ra/3-yS have common space breaking on the spherical surface at 
the distance r = r\, r = y/Z/xpo = ^Ja? jr g from the star. 
This is the geometric sense of the outer space breaking of the field of 
a star (in the framework of the liquid model under consideration) . 

§2.4 The force of gravity acting inside a liquid star 

In a rotation-free space, according to the definition of the gravitational 
inertial force (1-42), the force is only due to goo (which is determined 
by the gravitational potential w) . Let us calculate that force. Since the 
gravitational potential is w = c 2 (l — ^/goo), we obtain 




Cl212^-00, C1313 -> -00 (2.118) 



by the condition r = rv = y/'i/xpo = y/a 3 /r g . By the same condition 
r = Tb Tl according to the formula for i?oioi (2-113), we have 

#0101^-00. (2.119) 




<9w 



c 2 <9#oo 



(2.120) 




In the "inner" metric of a non-rotating liquid sphere (2.76) 



- 9n °- 4 ^V 1 3 V 1 3^ J 



(2.121) 



or, in the same metric written in the other form (2.78), 




) 



2 



3oo 



(2.122) 



58 



Chapter 2 Regular Stars and the Sun 



We therefore obtain that the force acting inside it is 

xp 0 c 2 r 1 

1_ ^ 7 ; / \ / (2.126) 

(aVi-^-V 1 "^) V 1 "^ 

or, in the other form, 

f.'- ^r, ^ = . = ■ (2-125) 



° 3 (V^¥-\A r ¥)^ 



«3 



F' = - ^f £ ^ _ . (2.126) 

This is a force of attraction: since r < a inside the sphere, F\ < 0 
therein. The force is proportional to distance r. Its numerical value is 
zero at the center of the sphere, then it increases with distance upto its 
ultimate-high value on the surface of the star (where r = a) 

_ xp 0 c 2 a 1 c 2 r g 1 
(ri) r =a — 7- 7 — 7TY 7~ ' (2.127) 



1 



(FW-^P = -^. (2.128) 



§2.5 Solving the conservation law equations: pressure and 
density inside the stars 

Consider now the pressure p and density po inside a regular liquid star. 
A formula connecting pressure and density inside a medium is the equa- 
tion of state. It follows as a solution of the conservation law equations. 

We have now already obtained almost all that is needed for the 
formula. In §2.1, we solved the conservation law equations with the 
energy-momentum tensor of a perfect liquid (2.4), which points to the 
substance of liquid stars. After substitution of the physically observable 
components (2.21) of the energy-momentum tensor, the general equa- 
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tions of the conservation law (1.89-1.90) take the particular form (2.26- 
2.27). In a non-deforming (static) space such as the space of a regular 
star, only the vectorial conservation equation remains non- vanishing. It 
has the form (2.36). The equation is solved as the formula (2.40) 

p = Be~% -p 0 c 2 . (2.129) 

Now, substituting the already found integration constant B (2.70) 
and function (2.75) into p (2.129), we obtain the final solution con- 
necting pressure p and density po inside a regular star 



„2 



p = p 0 c 2 ^^ =^- v , = . (2.130) 

Find the pressure in the near-surface layer of a star. The constant 
x= 18.6 xlO^ 28 cm/g is a very small value, while po—lA gram/cm 3 for 
the Sun, the yellow dwarf, and is much less than that for larger stars. 
Therefore, xp^a 2 is much smaller than 1 for even very large stars. For 
instance, for Betclgcusc, which is one of the largest red super-giants: 
M = 4.0xl0 34 gram, a = 7.0xl0 13 cm, po = 2.8xl0~ 8 gram/cm 3 . In this 
case, we have xp 0 a 2 = 2.6 xlO -7 . As a result, we have, for the values of 
r and a, 

^ _ xp Q r 2 _ 1 _ xpor 2 
3 6 

Thus, after some algebra, we obtain the approximate formula for the 
pressure p inside a regular star. It is 

Let h = a — r be the distance from the surface of the sphere to the point 
of measurement. Because h<^ r in the near-surface layer, we have 

a 2 -r 2 = (a-r) (a + r) = h (2a + ft) w 2 aft. (2.132) 

Thus, from (2.131), we obtain the regular formula for the pressure in 
the near-surface layer 

p = p 0 gh, (2.133) 
where ^f- — q is the free-fall acceleration in the star's field near its 

a 2 a 

surface. 
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Object 


Mass M, 


Radius 


Density p 0 , 


Pressure p 0 , 




gram 


a, cm 


gram / cm 


dynes/cm 2 


Red super-giant* 


4.0 xlO 34 


7.0 xlO 13 


2.8xl0~ 8 


5.3 x 10 s 


White super-giant* 


3.4 xlO 34 


4.8 xlO 12 


7.3 xlO" 5 


1.7xl0 10 


Sun 


2.0 xlO 33 


7.0 xlO 10 


1.4 


1.3xl0 15 


Jupiter (proto-star) 


1.9 xlO 30 


7.1 xlO 9 


1.3 


1.2xl0 15 


Red dwarfs 


6.7 xlO 32 


2.3xl0 10 


13 


1.3xl0 16 


Brown dwarf* 


4.1 xlO 31 


7.0 xlO 9 


29 


5.7xl0 15 


White dwarf § 


2.0 xlO 33 


6.4 x 10 s 


1.8xl0 6 


1.9 xlO 23 



*Betelgeuse. t Rigel. * Corot-Exo-3. § Sirius B. 



Table 2.2: The main characteristics of the regular stars. 



The pressure in the central region of a regular star can easily be 
found by assuming r = 0 in the general formula (2.130). Denoting the 
central pressure as p 0 =p r =o, we obtain 



1 , A *P° a2 2 2 2 

^ = ^ 7 7 (2-134) 



Since h = , we can also re- write this formula in the form 

3GM 2 , n . 

Table 2.2 gives the numerical values of the central pressure p 0 we 
have calculated according to this formula for the typical members of the 
known families of regular stars. 

We see that, according to our model of liquid stars, the pressure in 
the central region of Betelgeuse, which is one of largest stars, is only 
0.53 atmosphere (1 atm=10 6 dynes/cm 2 ). The smaller the size of a 
star is, the higher the pressure inside it becomes. The pressure in the 
central region of Rigel, a white super-giant whose radius is 14.6 times 
less than that of Betelgeuse, is 1.7 xlO 4 atm. In dwarfs such as the Sun, 
the central pressure is ~ 10 9 atm. However, in white dwarfs, the central 
pressure reaches 10 17 atm. 

Note that the temperature of condensed matter does not depend on 
pressure. The incompressible liquid of stars is a sort of condensed mat- 
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ter. Therefore, temperatures inside stars depend solely on the formula 
of that particular mechanism which produces stellar energy. 

This note is important for the understanding of the physical condi- 
tions inside stars, and of the sources of stellar energy. 

§2.6 The stellar energy mechanism according to the liquid 
star model and the mass-luminosity relation 

First, we make the transition to the dimensionless characteristics of 
stars, which are expressed in fractions of the respective characteristics 
of the Sun: 

M=^-, a=—, p=—, ... etc. (2.136) 

M Q a Q p Q 

where M = p 0 a 3 for a liquid sphere*. For the luminosity L of a star, 
that is the energy emitted from the entire surface of the star into the 
outer cosmos per one second, we have 

L=^-. (2.137) 

With this representation of the characteristics of stars, the analysis be- 
comes much simpler. This is because only the essential factors remain 
in the formulae while all constant coefficients vanish. 

Let us study what mechanism producing stellar energy can now be 
suggested due to the General Theory of Relativity, so that its produc- 
tivity satisfies the observed luminosity of stars. In other words, to be 
the real mechanism that generates energy in stars, the calculated en- 
ergy production of the suggested mechanism should match the mass- 
luminosity relation which is the main empirical relation of observational 
astrophysics. 

Consider thus the space metric of a liquid star. As we know already, 
the space of a liquid star has two primary regions which are described 
by different metrics: 

1) The internal space metric of the star (the metric of a liquid sphere) 
is valid from the center of the star to its surface. Except on the 
singular spherical surface of the tiny radius, r g = ^M- around the 
center of the star (see below) . The internal metric is also valid on 
the singular spherical surface of the radius n, r — \/a 3 / r g=\/3/ xpo, 
in the far cosmos: the metric produces a breaking of the space 
curvature at this distance from the star; 



*A liquid star has the same density p = po = const in its entire volume, so its 
mass is M= ^-irpoa 3 . In fractions of the Sun's mass, it is M = poa 3 . 
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2) The external metric of the star (the mass-point metric) is valid 
from the surface of the star to infinity. Except on the singular 
spherical surface of the radius, 7V = */ a 3 jr g = / >cpo which cov- 
ers the star distantly from its surface in the cosmos (see above). 
The external metric is also valid deep inside the star, on the sin- 
gular spherical surface of the tiny radius r g — 2G j^ from the center 
of the star: on this spherical surface, the star's gravitational field 
possesses space breaking produced due to the external metric. 

As was shown in §2.3, the outer space breaking in the far cosmos 
only implies space curvature breaking. One can show, on the basis of 
§2.3, that it does not result in an anomaly with respect to the acting 
force of gravitation. 

However, we now show that the force of gravitation has a very strong 
anomaly on the singular spherical surface of the inner space breaking. 
Actually, inside the star at the Hilbert radius r g from its center the 
external space metric is valid (while the internal metric is valid both 
inside the Hilbert radius and outside it). Therefore, all calculations 
for the inner singular surface are processed with the external metric 
(mass-point metric). This is despite the fact that the singular surface 
is located deep within the star near its center. 

According to the fundamental metric tensor of the external metric 
of the star (1.1), the chronometrically invariant (physically observable) 
vector of the force of gravitation F, has the form (1.4). On the singular 
spherical surface of the Hilbert radius r = r g , deep inside the star, the 
observable force of gravity (1.4) reaches an infinitely large magnitude 



1 

1 



F i = -^f (2-138) 



i.e. the gravitational field possesses space breaking on the surface. 

Due to its infinitely large magnitude there, the force of gravity, by 
definition, is sufficient for the transfer of the necessary kinetic energy 
to the lightweight atomic nuclei of the stellar substance, so that the 
process of thermonuclear fusion begins. The energy released in the 
thermonuclear fusion is the energy that the stars radiate. 

The singular spherical surface of the Hilbert radius r g — 2G J^ sur- 
rounds the geometric center of every star. This means that at the center 
of each star a luminous "inner sun" is located. The "inner sun" is tiny 
compared to the size of the star. For example, the Hilbert radius of the 
Sun is only 2.9 km while the physical radius of the Sun is 700,000 km. 
Therefore, the zone where thermonuclear fusion is processed is not only 
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the surface layer of the radius r g but all the volume of the "inner sun" . 
In other words, the "inner sun" of the radius r g is the very place where 
thermonuclear fusion produces helium from hydrogen, thus providing 
energy for the luminosity of the star. The energy is then transmitted 
from the "inner sun" of the star to its surface due to heat conductivity 
(the conventional transfer of heat in liquids); then it is radiated from 
the surface into the cosmos. 

Since the "inner sun" of a star has a radius equal to the Hilbcrt 
radius r g , we will further refer to it as the luminous Hilbert core of a 
star, or merely — the Hilbert core. 

The luminosity of a star that shines due to the suggested mechanism 
of stellar energy depends only on two factors: the volume of the Hilbcrt 
core V= |7rr-g where energy is released, and the density p g of the stellar 
substance therein (which can differ from the density p 0 of the main mass 
of the star, see the explanation in the next page) . In terms of the dimen- 
sionlcss characteristics of stars, it is 



Recall that the suggested mechanism of stellar energy does not de- 
pend on the pressure in the central region of the star: the super-strong 
force of gravity (2.138) that acts therein provides the conditions neces- 
sary for thermonuclear fusion. But the productivity of the mechanism 
depends on the density of the stellar substance in the Hilbert core. 

Calculate the density of the Hilbert core so that the suggested mech- 
anism of stellar energy satisfies the observed mass-luminosity relation. 

Proceed from the facts of observational astronomy. It shows the 
mass-luminosity relation L = M 2 6 for the stars whose masses are in the 
range between 0.2 M Q and 0.5M Q , L = M 4,5 for masses between 0.5M Q 
and 2M©, L = M 3 6 in the range between 2Mq and 10 Mq, and L — M 
for masses much heavier than 10 Mq. See Table 2.3. 

These empirical data of observational astronomy match with our 
theoretical formula for the luminosity of stars L (2.139), if the stellar 
substance of the Hilbert core (wherein stellar energy is released) has the 
density as shown in Table 2.4. 

On the basis of the function p g = M y according to Table 2.4, we are 
able to know how dense the Hilbert core of a star is compared to the 
main mass of the star (known from astronomical observations). We can 
thus calculate, for some typical stars, the following ratio: 





(2.139) 



Pg 



My 



(2.140) 
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Observed mass-luminosity 
relation L = M X 


Scale of the stellar masses, in fractions 
of the Sun's mass Mq 


L = M 2 - 6 


M = 0.2 . .. 0.5 


L = M 4 - 5 


M = 0.5 ... 2 


L = M 3 - 6 


M = 2 ... 10 


L = M 


M > 10 


Table 2.3: The observed mass-luminosity relation L — M x . 


Density of the Hilbert 
core p g 


Scale of the stellar masses, in fractions 
of the Sun's mass Mq 


p g = M 04 


M = 0.2 . .. 0.5 


P 9 = M 1 ' 5 


M = 0.5 ... 2 


P 9 = M 0 ' 6 


M = 2 ... 10 


p g = M- 2 


M > 10 



Table 2.4: Density of the substance inside the Hilbert core. 



Object 


Mass M 


Density p 0 


Ratio p g / Pq 


Betelgeuse (red super-giant) 


20 


2.0xl0" 8 


1.3xl0 9 


Rigel (white super-giant) 


17 


5.2xl0" 5 


6.7xl0 7 


Jupiter (proto-star) 


9.5 xlO" 4 


0.9 


0.069 


Red dwarfs 


0.34 


9 


0.072 


Brown dwarf (Corot-Exo-3) 


0.021 


21 


0.010 


White dwarf (Sirius B) 


1 


1.3xl0 6 


7.7xl0" 7 



Table 2.5: Ratio p g /po for some typical stars. 
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The results of the calculations are shown in Table 2.5. On the basis 
of the calculated ratio p g /p 0 as shown in Table 2.5, we arrive at the 
following conclusion. The luminous Hilbcrt core of a star — its "inner 
sun" — can have a density that differs from the density of the main 
mass of the star. It depends on the particular type of the star. For 
instance. The stellar substance of the Hilbert core of a giant or super- 
giant is many orders denser than the main substance of the stars. The 
Hilbert core of the star that is similar to the Sun has approximately 
the same density as the star. Concerning the dwarf stars, the Hilbcrt 
core of such a star is more rarefied than the main substance of the star. 
The greater the density of a dwarf star is, the less the density of its 
core becomes compared to the density of the entire star. In such a star 
as the white dwarf, the Hilbert core is many orders of magnitude more 
rarefied than the main substance of the star. 

Respectively, the following question arises. All physical bodies have 
masses, therefore each body should have a Hilbert radius core inside 
itself. Not only stars, but also planets and even individual elementary 
particles should have such a core. Yet, why do they not shine like stars? 

The answer comes from the state of that substance of which these 
physical bodies consist. Stars are made up of liquid substance which 
consists, mostly, of light chemical elements such as hydrogen and helium. 
Therefore, thermonuclear fusion of light atomic nuclei is possible in the 
Hilbcrt core of each star. Due to the fact that the substance is liquid, 
more and more "nuclear fuel" is delivered from the other regions of the 
star to its luminous Hilbcrt core, thus supporting the combustion inside 
the "nuclear boiler", until the time when all the nuclear fuel of the 
star ends. Another case — the planets. They consist of mostly heavy 
elements with only a minor content of hydrogen. Therefore, as soon as 
the "nuclear boiler" of the Hilbcrt core has finished all the reserve of 
the hydrogen fuel in the central region of the planet, it stops producing 
energy but still remains to exist at the center of the planet, in a latent 
state. 

Astronomers know that the energy emitted by Jupiter exceeds the 
solar energy absorbed by the entire surface of the planet. The same is 
as well true for Saturn. This means, according to our theory, that the 
Hilbcrt core of each of the planets still processes hydrogen into helium 
thereby releasing nuclear energy. 

Concerning individual elementary particles such as protons, neurons, 
and electrons: as is known, they are stable and indifferent for a long 
time as long as they do not interact with other particles. In fact, this 
means that the Hilbcrt core of the proton (as well as of the neutron and 
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the electron) does not interact with the main mass of the particle. Why 
does this happen? We can only guess that either the substance that is 
inside the particles is in the super-solid state, or there is a layer of the 
very strong vacuum between the core and the rest mass. On the other 
hand, the Hilbert core of the proton (and that of the neutron) has a tiny 
radius of (r g ) p = =^ = 2.48 xl(T 52 cm, while the Hilbert core of the 
electron has even a smaller radius of (r 9 ) e = 2 " 6 = 1.35 xlCP 55 cm. As 
has been stated by Albert Einstein already, the geometric laws (space- 
time geometry) of the General Theory of Relativity are true, probably, 
upto the scale of elementary particles. Within the sub-nuclear scale, 
probably, another geometry works thus stating its own laws which differ 
from the laws of the General Theory of Relativity. Therefore, we cannot 
presently state something definite about the physical conditions and 
processes inside elementary particles. 

But as for the regular world of stars and the planets, experimental 
physics and observational astronomy show that Einstein's theory is cor- 
rect and works on these scales with high accuracy. Therefore, all our 
conclusions about the internal constitution of stars, and the mechanism 
that generates energy in stars should be taken into account. 

The particular details of the suggested mechanism of stellar energy 
are a special theme that is out of the scope of this book (which is mostly 
on the internal constitution of stars). 

§2.7 Conclusion 

All the theoretical conclusions about the source of stellar energy, and 
about the internal constitution of stars that are presented in this Chap- 
ter have been obtained in the framework of our model of liquid stars. 
Our model is based on the presentation about stars as space-time ob- 
jects, according to the General Theory of Relativity. Below, we list the 
most important of the conclusions we have thus arrived at: 

1. The field of each star possesses space breaking which surrounds 

the star by a spherical surface. The "bubble" of the outer space 

breaking of the field has a radius of 



which is many orders larger than the physical radius a of the 
star. The three-dimensional observable space curvature Cikij and 
the four-dimensional Riemannian curvature R a p~/8 have common 
space breaking on the surface. The outer space breaking impedes 
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the near substance to be formed as a planet in this orbit. The outer 
space breaking of the Sun's field is located within the asteroid 
strip, near the maximal concentration of the asteroids; 

2. The field of each star possesses inner space breaking, inside the 
physical body of the star, on the surface of the Hilbcrt radius 

_ 2GM 

from the center. This means that there is a small core which is 
separated, by the singular surface, from the major mass of the 
star. On the surface of the core, the force of gravity reaches an in- 
finitely large magnitude. The super-strong gravity, by definition, 
is sufficient for the transfer of the necessary kinetic energy to the 
lightweight atomic nuclei of the stellar substance, so that ther- 
monuclear fusion begins. Thus, nuclear energy is released. The 
liquid "nuclear fuel" is delivered from the other regions of the star 
to the core thus supporting the combustion inside the "nuclear 
boiler" ; 

3. Every star has a mass M. Therefore, the luminous core of the 
Hilbert radius r g — "^f- — the "inner sun" — exists in the center 
of every star. We refer to it as the Hilbert core. This is the place in 
which thermonuclear fusion produces helium from hydrogen, thus 
providing energy for the luminosity of the stars. The energy is then 
transmitted from the "inner sun" of the star to its surface due to 
heat conductivity (the conventional transfer of heat in liquids) for 
it then to be radiated into the cosmos; 

4. The Hilbert core is tiny compared to the size of stars. For example, 
for the Sun, r g ~2.9 km; 

5. The observed relation "mass- luminosity" of stars is satisfied if the 
Hilbert core has a density depending on the particular type of the 
star. The Hilbert core of a giant or super-giant should be many 
orders denser than the main substance of the stars. The Hilbert 
core of a star like the Sun should be approximately the same in 
density as the star. In the dwarf star, the Hilbert core should be 
more rarefied than the main substance of the star (the core of the 
white dwarf should be extremely rarefied); 

6. Each planet has a mass. Therefore, the Hilbert core exists in 
the center of every planet. But planets consist of mostly heavy 
elements with only a minor content of hydrogen. As soon as the 
"nuclear boiler" of the Hilbert core has finished all the reserve 
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of the hydrogen fuel in the central region of the planet, it stops 
producing energy but still remains to exist at the center of the 
planet, in a latent state. 



Chapter 3 
Regular Stars. The Description 



§3.1 Problem statement. The internal space metric of a reg- 
ular non-rotating star 

To broadly understand the description of a regular star again, recall 
that in §2.1 we deduced the space metric of a liquid sphere by follow- 
ing the "historical path" as Schwarzschild did it. Namely, — we took 
the spherically symmetric metric in the general form, then applied the 
particular conditions of a sphere hlled with perfect liquid. The sole dif- 
ference from Schwarzschild's deduction was that we did not assume any 
artificial limitations. When following this deduction, we obtained the 
observable characteristics of the space in the implicit form, as an auxil- 
iary result. It was enough to obtain the space metric of a liquid sphere 
in the final form. Now, we express the characteristics in the explicit 
form, through the components of the fundamental metric tensor of the 
space metric of a liquid sphere which we have obtained in Chapter 2. 
Then we will study the equations of motion inside the star so as to take 
the escape velocity into account. 

So, the space metric of a liquid sphere has the form (1.8) 

-_*±--r»(<» 3 + staW). (3-1) 

i ' 1 a 

We calculate the chr.inv. -characteristics of the space, according to their 
definitions given in §1.3 and the respective components of the funda- 
mental metric tensor of the metric (3.1). The chr.inv. -metric tensor hit 
of the metric (3.1) has the following non-zero components 

h u = K— , h 22 = r 2 , h 33 = r 2 sm 2 6, (3.2) 

_ r r g 
1 a 3 
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while its determinant, and the non-zero spatial derivatives of the loga- 
rithm of the determinant have the form 

r 4 sin 2 # 

h = det\\h ik \\ = — , (3.4) 

1 T T B 

1 a 3 

*dhiVh 2 r„r 1 *dlnVh 

—5 = - + ^T 5— . an =cot6. 3.5 

dr r a 3 i _ fz± 39 

1 a 3 

So forth, after algebra according to the chronometrically invariant for- 
malism (sec §1.3 for the definitions of the chr.inv. -quantities), we ob- 
tain the following. The chr.inv. -vector of the gravitational inertial force, 
acting in the space has the form 

2 

r,= "^(^-?-/-*j^-¥' <3 ' 6> 

f . = 4^£? , (3 . 7) 

where r <a since it is inside the sphere. Therefore, F\ < 0 therein (this 
means that it is a force of attraction). The non-zero chr.inv. -Christoffcl 
symbols have the form 



A i _ r s r 1 Ai __^k --rfl - r -ill\ 

Al1 " a 3 !_!£!!' A22 ~ sin 2 fl~ { « 3 J' 1 ' 

a 3 

A 2 2 = A? 3 = -, A^ 3 = -sin0 cos 6, A^ 3 =cot6». (3.9) 

The non-zero components of the chr.inv. -tensor of the three-dimensional 
observable curvature Cikij, and its contraction dk, have the form 

Cisia r„r 2 



r _ -1313 r fl r i _ r fl r 2 

l-1212 - — -2T - — -3 rrr, C2323- — -3- sin V, (3.1U) 

sin 6 a 1 — 8 a 



7.3 



C\\ — 3- 5- , C22 — — — ^3 — • (3-11) 

a 3 -| _ l£l! ' sin 2 6» a 3 v 7 

n 3 



So, we now have all that is needed to consider Einstein's equations 
in the internal field of a regular non-rotating star. 
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§3.2 Einstein's equations in the internal field of a regular 
non-rotating star 

Consider Einstein's field equations in the space of the metric (3.1). As 
is known, the energy-momentum tensor of a perfect liquid has the fol- 
lowing general form (2.4): 

T a ? = (po + J) U a U p - J g a P, (3.12) 

where po = const is the density of the liquid, p is the pressure, while U a is 
the four-dimensional velocity of the flow of the liquid with respect to the 
observer (the unit four-vector, so U a U a =l). The chr.inv. -projections 
of the energy- momentum tensor have the form (2.21) 

p=^= Po , J l = ^==0, U ik =c 2 T lk =ptf k , (3.13) 

where p is the observable density of mass, J 1 is the observable density 
of momentum, while U lk is the observable stress tensor. With these 
formulae, and by taking into account that the space of the particular 
liquid sphere is free of rotation and deformation (Ai k = 0, D ik = 0), the 
chr.inv. -Einstein equations (1.85-1.87) take the form 

*VjF j 1 Fj F 3 =-| ( Po c 2 + U) , (3.14) 

J*=0, (3.15) 

I (*V 4 F fc + *V k F t ) - ^ FF k - c 2 C lk = 

= 77 (poc 2 hik + 2U lk - Uh ik ) , (3.16) 

where *Vi is the symbol of chr.inv. -differentiation (see Notations), while 
U ik = ph ik and U = 3p. 

Substitute, into the Einstein field equations, the formulae for F i7 d kl 
and hi k , which have been calculated for the metric (3.1). We obtain that 
only two equations remain non- vanishing: 
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Multiplying (3.18) by 3 then summarizing the product with (3.17), we 
obtain 

W 2 - ^P- • (3-19) 

a 6 

Substituting this result back into (3.18), we obtain the equation of state* 
for the liquid substance of regular stars 



p = Po c 2 ^^= -— v = . (3.20) 

This formula completely coincides with the formula for pressure p 
(2.130) we have obtained in Chapter 2 as the result of following the 
path of Schwarzschild. 

This formula for pressure p can also be obtained from the conserva- 
tion equations (2.26-2.27). Because the space of the metric (3.1) does 
not deform (hik ^ f (t), hence = 0), the chr.inv. -scalar conservation 
equation (2.26) vanishes. Only the chr.inv. -vectorial conservation equa- 
tion (2.27) remains. It takes the form 

*Vi(pti k )- ( Po + |r)^ fe =0- (3.21) 

Herein, *Vih lk — 0 is true always as well as V a g ac7 = 0 for the funda- 
mental metric tensor. Therefore, and because the chr.inv. -derivative 
with respect to the spatial coordinates coincides with the regular spa- 
tial derivative in the case where the space does not rotate, the remaining 
conservation equation (3.21) has the form 

hik Z-( p ° + ^) Fk=o - (3 - 22) 

Substituting the formulae for h 11 and F 1 we have obtained for the metric 
(3.1), we transform (3.22) into the differential equation 



dp r g r dr 



(3.23) 



This equation can be re-written in the form 




din ( Po c 2 +p) = -din 3Wl - - dl - , (3.24) 



r„r 2 



"The formula connecting pressure and density inside the medium. 
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which is easy to integrate. After integration, we have 



p + Po c 2 = — Q =, (3.25) 



3a /1_ l£l_ 



where the integration constant Q comes from the obvious condition p = 0 
on the surface of the star (where r = a). Then 



Q = 2p 0 c\ll-^ , (3.26) 

thus we obtain, finally, 



'1 _ 1£ 

p + p 0 c 2 = 2p 0 c 2 L 2 (3.27) 

It is easy to see that this solution leads to the same formula for p as 
(3.20), which we have obtained from Einstein's field equations. 

§3.3 The internal space metric of a regular rotating star 

Consider now the metric of a regular liquid star (3.1) with only the 
change that the star rotates, at an angular speed oj, along its equatorial 
axis (the axis <f> in the spherical coordinates r, 8, <fi). In this case, the 
non-rotating metric (3.1) takes the form 

. 2 



1 




1-2 ) 



^ = ^134/1-^-1/1-^ c 2 dt 2 + 



dr 2 

cdt dtp 5 r 2 (d9 2 + sm 2 0 d<f> 2 ) . (3.28) 

-i r r g 

a 3 

It should be noted that we still consider regular stars, namely those 
stars whose Hilbert radius is much smaller than their physical radius. 

According to the metric (3.28), the star's linear velocity of space 
rotation is 

2ur 2 cos6 

Vl =v 2 =0, v 3 = == == ■ (3.29) 



r g r' 
„3 



As is known from the data of observational astronomy, the major- 
ity of stars rotate at linear speeds v < 420 km/sec. Hence, we have 
v 2 /c 2 < 2 x 10~ 6 : most stars rotate slowly compared to the speed of light. 
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According to the space metric of a regular (slowly rotating) regular 
star (3.28), represented as a rotating liquid sphere, we have 

v 2 = h ik v iVk = h 33 v 3 v 3 , h 33 = - g 33 = 1 . (3.30) 

r z sm 0 

So, v 2 /c 2 in the space of the metric (3.28) has the form 
v 2 Auj 2 r 2 cot 2 0 



3a /1_ IiT 



2 • 



(3.31) 



Expanding, in this formula, the radicals into series, after elementary 
transformations we obtain 

Further, we will neglect higher-order terms of the series which are small 
because r 9 C« for the regular stars. Therefore, we have 

, v 2 v 2 oj 2 r 4 cos 2 6 
v 3 =ujr 2 cos9, — = -| = . (3.33) 

The non-zero components of the chr.inv. -metric tensor of the metric 
(3.28) have the form 

hu = ttt = "1 — > ^22 = t^-=^ 2 , (3.34) 

h 11 i _ r 2z± h 22 

, 1 2fl / o; 2 r 2 cot 2 ^ 

^33 = 7^7 = r Sln 0 ^ + ~2 ) > ( 3 - 35 ) 

while the determinant of the chr.inv.-metric tensor /i^ and the non-zero 
spatial derivatives of the logarithm of the determinant have the form 

, ,,, „ r 4 sin 2 9 / oj 2 r 2 cot 2 9\ 
fc = det|MI = -3- 1 + -= , (3.36) 



*d\ny/h 2 



<9r r a 3 i _ l!l£ 

1 „3 



(3.37) 



uj 2 t 2 



* 9} ^=™te\l-^- 2 = |. (3.38) 

96* | c 2 sin 2 0 1+ W V cot2 e 1 
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Respectively, according to the chronometrically invariant formalism 
(see §1.3 for the definitions of the chr.inv. -quantities), we obtain also 
the other chr.inv. -characteristics of the space. The chr.inv. -vector of 
the gravitational inertial force Fi, acting in the space takes the form 



Fi = — 



F 1 



c°r a 



r\ 1 



<0, 



<0, (3.39) 



(3.40) 



which is a non- Newtonian force of attraction. The approximate formula 
for the force is 



F 1 = F 1 w - 



c r g r 
2a 3 



(3.41) 



The chr.inv. -tensor of the angular velocity of space, Aik, has the follow- 
ing non-zero components: 

2ujr cos 9 



A 1Z = 



A 13 = 



a3 (3^-0 T ¥)v / ^ T ¥ 

2lu (l- Ig^cotfl 



(aVT^-Vl-^rain^l 



uj 2 r 2 cot 2 9 



r a r- 



_ a3 (3^"\/^¥)0 T ¥ 



-423 = 



uir sm( 



123 



(3.42) 
, (3.43) 

x (3.44) 

, (3.45) 

(3.46) 
. (3.47) 
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The approximate expressions of these components have the form 

A 13 = -urcose(l + ^-?jpJ , (3.48) 

LJCOtO ( 3r„ 2r Q r 2 uj 2 r 2 cot 2 6 \ 

^ = : — 2T ( 1 + ~T~ ~ — 1 2 - 3.49 

rsm& \ 4a a A c z ) 

cur 2 sind / 3r 0 r 0 r 2 \ . „. 

As » = 2 ( 1+ 4^"^-)' (3 ' 50) 



(23 



l + fe-^-" V f 2g ). (3.51) 



2r 2 sin# V 4a 4a 3 



The non-zero chr.inv.-Christoffcl symbols for the metric, with higher- 
order terms cj 4 r 4 /c 4 neglected, have the form 

A -=¥rW ^— K 1 -^)- <3 - 52) 

a 3 

*i .2. A 2w 2 r 2 cot 2 fl\ /, r„r 2 \ . . 

Ai 3 = -rsin 2 0 1 + = 1--V ' ( 3 - 53 ) 



* 9 1 1 /, w 2 r 2 cot 2 1 . 

(2 2 \ / 2 2 
l-^J, =cotg (1- ^ r ) . (3.55) 

Neglecting the terms r 2 /a 2 and the product w 2 r 2 /c 2 to r 9 /a, we obtain 
the approximate formulae for these components 

A n = ¥' A22 = " r ( 1 - 1 ^)' ^ 
a i • 2 „ ( 2uj 2 r 2 cot 2 9 r Q r 2 , 

A 33 = -rsm 2 ^l + *L_ ) , (3..,,) 

a ? 1 A o 1 ( _ oj 2 r 2 cot 2 L , 

A 2 2 = - , A? 3 = - 1 + , (3.58) 



r r 



w 2 r 2 



A 2 3 = -sin0cos0( 1 — ), A 3 3 =cotg(l- e2 — 2 - \ , (3.59) 
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and the non-zero components of the chr.inv. -curvature tensor Cikij, 
along with the non-zero components of its contraction dk 



C1212 



1 



1 - 



C1313 = r 2 sin 2 9 



3u; 2 cot 2 0 r g 1 + 



V cof 9 



(3.60) 
(3.61) 



C2323 = 



On = 5- 



1 + 



uj 2 r 2 cot 2 



+ 



cot 2 0 + 



sin 



+ 



3w 2 cot 2 0 



: sin 2 0, 



(3.62) 
(3.63) 



C22 — 
C33 = 



2r„r 2 u 2 r 2 cot 2 9 
1 

a 3 c 2 

2r„ 



cot 2 9 + 



sin 4 (9 



a" 



w 2 r 2 cot 2 ( 



3w 2 cot 2 0 w 
x - h 



2 / 1 

cot 2 0 H j 

sin 



2 • 2 
r sin 



(3.64) 



(3.65) 



Because r 9 /o<Cl and w r /c <?C 1 for regular stars, we neglect the 
terms r 2 /a 2 and the product uj 2 r 2 /c 2 to r g /a. As a result, we obtain 
the approximate formulae for the chr.inv. -curvature: 



C1212 : 

C2323 : 

C n = 
C22 = 

C33 = 



C 13 i 3 = r 2 sin 2 0 



3w 2 cot 2 0 r 



9 

3 ; ' 



+ 



a" 



cot 2 9 



1 



sin 4 9 



2 • 2 
r sin t 



2r g 3w 2 cot 2 0 
2r„r 2 Lo 2 r 2 cot 2 1 



+ 



cot 2 0 + 



1 



2r g 4w 2 cot 2 0 



+ 



c 2 sin 0 



sin 0 / 
r 2 sin 2 0. 



(3.66) 

(3.67) 

(3.68) 
(3.69) 

(3.70) 
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§3.4 Einstein's equations in the internal field of a regular 
rotating star 

We now solve Einstein's field equations in the internal space of a rotating 
regular star, i.e. in the space of the metric (3.28). In the absence of 
space rotation (Ai k = 0), this problem reduces to the problem that was 
considered earlier in §5.2 for regular non-rotating star. 

So, consider the chr.inv. -Einstein equations (1.85-1.87) in the space 
of a liquid sphere, which rotates (A ik ^0) but is free of deformation 
(Dik = 0). They take the form 



— Vj A lJ + — Fj A 13 =xJ\ 



2^ 



(3.71) 
(3.72) 



■ (*ViF k + *V k Fi 



■ FiF k — c Cik — 



— {pc 2 h lk + 2U ik - Uh lk ) 



(3.73) 



where *Vi is the symbol of chr.inv. -differentiation (see Notations). The 
chr.inv. -quantities p, J 1 , U lk are the physical observable projections of 
the energy-momentum tensor T a p of the medium that fills the space. 
(We do not specify the energy-momentum tensor just yet.) 

With the obtained components of Ai k and Fi (see §3.2 for detail), 
the chr.inv. -Einstein equations (3.71-3.73) take the form 



2w 2 cot 2 6> 



1 | 3r g ' ir 9 r2 

^ 2a a 3 



+ =r(l + ^- 



2a 



r g r 2 oj 2 r 2 cot 2 8\ 3c 2 r, 
~2a* c*~ 



+ 



>COt# 



3Tg 

4a 



4r ff r 2 



2 a 3 2 
3w 2 r 2 cot 2 6> 



rJ 3 , 



(3.74) 
(3.75) 



<> ( I - ^ 

2 a a 3 

(pc 2 -U)(l- 



2~ 



u 2 r 2 cot 2 i 



+ 2U n 



2 a 



5r g r 
4 a 3 



uj 2 r 2 cot 2 i 



3c 2 r g 
2 a 3 



->rfJi2 



(3.76) 
(3.77) 
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uj 2 r 2 



r n r 2 oj 2 r 2 cot 2 9\ 3c 2 r n r 2 



2a 2a 3 c 2 
(pc 2 -U)r 2 +2U 22 



+ 



2 a 3 



2w 2 cot 2 dl + ^ 



3r„ 3r„r z 



2a 

3 c r g r sin # 



+ ■ 



2a 2a 3 



(3.78) 



r 2 sin 2 9 + 



( pc 2 ~U)r 2 sin 2 <9 + 2[/ ; 



:s:s 



(3.79) 



In the framework of our approximation (r g /a^;l and uj 2 r 2 /c 2 <C 1) 
which holds true for regular stars, we neglect the terms r 2 /a 2 and the 
product uj 2 r 2 / c 2 to r g /a. As a result, the obtained chr.inv. -Einstein 
equations take the much simplified form 



^9 9 n ^ K 

2uj 2 cot 2 9 + — = - 
2 2 


(pc 2 + 




(3.80) 


u> cot 9 3 
sin w 






(3.81) 


2w 2 cot 2 6»- - (pc 2 - 
2 VP 


-^) = 




(3.82) 


uj 2 r cot 9 = —kU\i , 






(3.83) 


-2-2^-^ = 


xC/f , 




(3.84) 


2w 2 cot 2 6»+ — - - 
2 2 


(pc 2 - 


[/) = x£/ 3 3 . 


(3.85) 



Summarizing (3.82), (3.84), and (3.85), then taking into account 
that U\ + U$ + U$ = U, we obtain 



Alo 2 cot 2 9 + uo 2 = ^ (3pc 2 - U) . 
Summarizing the result (3.86) and (3.80), we obtain 

3w 2 cot 2 9 + — — = xpc 2 . 
4 



(3.86) 



(3.87) 



Multiplying (3.80) by 3 then subtracting the result from (3.86), we 
obtain 

(3.88) 



w 2 cor 9+ — = kU. 
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As is seen from (3.87) and (3.88), 

pc 2 = 3U. (3.89) 

The energy-momentum tensor T a /j should also satisfy the conserva- 
tion equations. The chr.inv. -conservation equation in the rotating space 
of the star, which is free of deformation (Da- = 0), have the form 

*-^ + *V t J l -^F l J l = 0 1 (3.90) 
at c z 

*8 J k 1 

-— + 2A'KP + *V, U lk --^Fi U' lk - P F k = 0 . (3.91) 
at c z 

As follows from the scalar conservation equation (3.90), 

= 0 , i.e. p = const. (3.92) 

The vectorial conservation equation (3.91) with the index i — 3 is satis- 
fied identically. The equations with i = 1 and i = 2 take the form 

+ A 33 t/ 33 + ( A}i + _ JL t/n = ^i, (3.93) 

„, 2t s 9C/ 12 dU 22 d\nVh Tr22 

2A - J + — + ^r + — u + 

+ A 2 33 U^ + ( 2A 2 2 + ^ - 1 Fx) C/ 12 = 0. (3.94) 

In the tensorial chr.inv. -Einstein equations that readily express U\ 
(3.82), C/12 (3.83), f/| (3.84), and C/| (3.85), we now take into account 
that pc 2 = 3U (3.89) and the formula for U (3.88). We obtain 

UJ 2 

hU 11 = uj 2 cot 2 9 - — , (3.95) 

TT i9 ^ 2 cot 9 

xU 12 = , (3.96) 

r 

99 uj 2 uj 2 cot 2 9 

^33 = "'cot'fl + f (3 98) 

r 2 sin 0 4r 2 sin 6 
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Substitute these formulae, and also the other necessary quantities, into 
the remaining conservation equations (3.93) and (3.94). After some 
algebra, we see that the equations are satisfied identically. 

So, Einstein's equations and the conservation equations in the form 
considered herein satisfy the metric of the internal space of a regular 
rotating star, that is the space metric (3.28). 

§3.5 The stationary vortex-free electromagnetic field of a reg- 
ular rotating star 

A realistic star bears an electromagnetic field. Therefore, we should 
introduce the electromagnetic field into the theory of liquid stars. Elec- 
trodynamics in terms of the chronometrically invariant formalism was 
introduced in Chapter 3 of our book [18]. We now follow the deduction 
therefrom, then apply it to the present theory of liquid stars. 

So, as is known from the generally covariant formulation of electro- 
dynamics [20] , the energy-momentum tensor of an arbitrary electromag- 
netic field has the form 

T ™ = (-F a a F^ + \ g a0 F^F^ , (3.99) 

where F Q( g is the electromagnetic field tensor known also as the Maxwell 
tensor. The field tensor F a p is the curl of the four-dimensional electro- 
magnetic potential A a 

F a , = V a A,-V,A a = ^-^. (3.100) 

The physically observable chr.inv. -projections of the four-dimensional 
electromagnetic potential A a are the scalar electromagnetic potential ip 
and the vector electromagnetic potential q l : 

<P=-^=, q l = A\ (3.101) 
V3oo 



The electromagnetic field tensor F a p (5.44) has the following physically 
observable components 

Pom = — = fi * , 3.102 

.9oo 8ttc 2 

4n = = -!- e* km E k H* m , (3.103) 

V300 47TC 

Ult = c 2 T lk = Pcm c 2 h lk - -!- {E l E k + H* l H* k ) , (3.104) 

47T 
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where E % is the three-dimensional chr.inv. -electric strength vector, H* % 
is the three-dimensional chr.inv. -magnetic strength pseudo-vector of 
the field, while e lmn is the completely anti-symmetric unitary three- 
dimensional chr.inv. -pseudo-tensor (see [18], for detail). Namely, - 



pHk _ _Ak nj? p _*difi I *8q n _tp 

OX™ C C z 

TT*i _ 1 -imn tt tt _ _ _ 4 

2 m " ' m " ~ &e» fe™ c 



> . (3.105) 



As is seen from the definitions (3.105), the chr.inv. -electric field strength 
and the chr.inv. -magnetic field strength depend on not only the electro- 
magnetic field potentials tp and q l , but also on the characteristics of the 
space that is filled with the field. The physical and geometric character- 
istics of the space that affect the electric and magnetic field strengths 
are, respectively, the gravitational inertial force Fj acting in the space, 
and the angular velocity of the space rotation . 

Herein, in §5.3, we first consider the electromagnetic field which is 
stationary and vortex-free. This means that the electromagnetic field 
potentials are constant (ip = const, qt = const) while the curl of the vec- 
torial potential is zero {qik =0). In other words, in the further consid- 
eration of the electromagnetic field we will assume that 

*dip *dq n *d qi *dq k fo ma, 

in the formulae for the electric and magnetic field strengths (3.105). In 
this case, we have 



E l = -^F\ Ei = -^Fi 

TT*i _^£_ O*^ / / . Q 

ii 1 ' . STL -m j - i £ * 



(3.107) 



where Q,* 1 is the three-dimensional chr.inv. -pseudo- vector of the angular 
velocity of space rotation 

0*i _ 1 imn a q _ 1 \mn /o l no\ 

It is easy to find that, in the internal space of a regular rotating star, 
that is the space of the metric (3.28), we have 

ST 1 = SU 1 = %, n *2 = ^^l i Q* 2 =u;rcot0. (3.109) 
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thus 

= u? Q + cot 2 6»^ . (3.110) 

As is seen from the formulae (3.107), in the stationary vortex- free 
electromagnetic field the electric field strength E % is determined by the 
scalar electromagnetic potential ip and the gravitational inertial force 
F l acting in the space. The magnetic field strength H* 1 is due to the 
vectorial electromagnetic potential q 1 and the angular velocity Q,* 1 of 
rotation of space. 

Using these formulae for E l and H* 1 (3.107), we obtain that the 
physically observable components (3.102-3.104) of the electromagnetic 
field tensor F a p have the form 

/5cm= 2&(^ + ^"^)' (3 ' m) 
JL = ^£ ikm Fk^ m , (3.112) 

,.2 / 771 jpj \ 2 / rpi jpk \ 

It is easy to see that the trace of the electromagnetic field stress 
tensor, which is expressed as U cm ~ hikU^, satisfies the condition 

U cm = ^ (^f- + = PemC 2 . (3.114) 

As follows from the general form of the energy-momentum tensor 
T a fj that satisfies the metric (3.28), the tensor should satisfy the condi- 
tion pc 2 = 3/7 (3.89). This formula differs from p cra c 2 — U cm (3.114) we 
have just obtained for the stationary vortex-free electromagnetic field. 
We therefore should find such a structure of the electromagnetic field 
that makes U em satisfying pc 2 = 3U. 

So, according to the condition pc 2 = 3U (3.89) and the formula 
ijj 2 cot 2 9+ jlj 2 = xU (3.88) obtained from the chr.inv. -Einstein equa- 
tions, the field density inside a regular rotating star should be 

p= 6U *i\ l U= l ±2±}—. (3115) 

xc 2, x 

Therefore, we substitute the required condition 

Uem = • (3.H6) 

X 
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into (3.114) that we obtained for the stationary vortex- free electromag- 
netic field. We obtain 

^?(^b- + ^ n 0 = ^-' (3A17) 

or, expanding Einstein's constant of gravitation x= ^r~, in the equiv- 
alent form 

c 2 n^il* j = c - FjF j . (3.118) 

^ _ 4Gy 2 



We have considered a stationary electromagnetic field: the scalar 
and vectorial electromagnetic potentials remain unchanged (<p = const, 
Qi = const). Therefore, in the formula (3.118), the quantity 

/~i 2 i 

-j"=n, n<- (3.119) 

cr 4 

is a constant dimensionless coefficient depending only on the scalar po- 
tential ip of the electromagnetic field of the star. 

Using the constant n (3.119), we re-write (3.118) as 

77 1 

Jtl.iSl* 3 = - — F t F\ n<-. (3.120) 

1 - 477 4 

Substituting QjCt j (5.13) and F t F l (2.15) into the condition (3.120), we 
can present it in the alternative (expanded) form 

(1+4 cot 2 0) = - 4 " c2r Y . (3.121) 
v ' 1 - 4n a 6 

If n = n max = j and thus <p — i/) max , the angular velocity of the star 
would be to = oo which is non-sense. Therefore, we conclude that n<\ 
for all real stars including the Sun. With n < j , we obtain the upper 
boundary of the numerical value for the scalar electromagnetic potential 
possessed by a real star 



..2 



<p = -^t= < 1.74 xlO 24 
2VG 



n 1 /2 f . m l/2 



(3.122) 



With the condition (3.120), the stationary vortex- free electromag- 
netic field satisfies Einstein's equations and the metric of the internal 
space of a regular rotating star. In other words, with this condition 
(3.120), a stationary rotating regular star is a permanent magnet. 
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Because n — const in the stationary electromagnetic field (wherein 
ip = const, qi — const), the above condition (3.120) allows us to express 
the characteristics of the electromagnetic field through the geometric 
and physical characteristics of the space. In other words, in this partic- 
ular case, we can geometrize the electromagnetic field. 

Substitute the obtained formulae for the gravitational inertial force, 
and for the angular velocity of space rotation into the physically observ- 
able components (3.111-3.113) of the electromagnetic field tensor F a p. 
By taking the relations (3.117) and (3.119) into account, we obtain the 
observable components of F a p in the form 



U, 



li 



u, 



u 



12 



2ttG 



c 2 r 2 r 2 / 1 

—V + w 2 - + cot 2 1 
16a 6 \4 



P - - 

,J om 



nc 2 ujr g cot# 



4ttG a 3 sm9 ' 

,2 



nc 
2ttG 



c 2 r 2 r 2 
16 a 6 



nc 2 u) 2 cot# 
2ttG r 

„2 



nc 



2nGr 2 



9 9 9 

c r g r 
16 a 6 



- OJ 



cot 2 6» 



— - cor 

4 



(3.123) 
(3.124) 
(3.125) 

(3.126) 
(3.127) 



U 33 = 



2irGr 2 sin 2 9 



— V + 7 + cot2 1 
16a 6 \4 



From these equations we obtain, as previously, 

U C m — hik ^cm — Pern C • 



. (3.128) 



(3.129) 



The chr.inv. -Einstein equations (3.123-3.128) can further be simpli- 
fied. In the surface layer of a star (r w a), the first term in the brackets 
is 



c 2 r 2 r 2 



2 2 

c r 
° ' g 



16a 6 " 16a4 " ^ 
Consider the Sun as an example. Its surface layer makes one full rev- 
olution with a period of 27 days, which is equivalent to the angular 
velocity of rotation cj q ~2.7><10~ 6 sec -1 . Therefore, the second term 
in the brackets is 

i u% ~ 1.8 xlO" 12 sec" 2 . (3.131) 
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The first term in the brackets, while taking the Hilbert radius for the 
Sun r gQ =2.9xl0 5 cm, and the Sun's physical radius a Q = 7.0xl0 10 cm 
into account, is ten times less: 

~ 2.0xl(r 13 sec" 2 . (3.132) 

Therefore, we can neglect the first term in the brackets for even slowly 
rotating stars such as the Sun. Thus the chr.inv. -Einstein equations 
(3.123-3.128) take the simplified form 

Pcm = ^ Q+cot 2 ^ , (3.133) 

3 __^^COt0 

J ° m " AirG «3 sin^' (6A6V 
^ = -^Q-cot 2 ^), (3.135) 

UZ = -^^l, (3.136) 

U™= ~ n /f. 2 (\+co^e). (3.138) 
27rGH sin 9 \4 / 

Nevertheless, the first term in the brackets of the chr.inv. -Einstein 
equations (3.123-3.128) can be sufficient in the case of small stars such 
as white dwarfs or brown dwarfs. This is because we have a 6 (the star's 
radius in the 6th power) in the denominator of the term. 

§3.6 Solving Maxwell's equations in the vortex-free electro- 
magnetic field of a regular rotating star 

As is known, the electromagnetic field is described by Maxwell's field 
equations. They consist of two groups of equations. The generally 
covariant formulation of Maxwell's equations in the four-dimensional 
pseudo-Riemannian space has the form [20] 

V a F^ = —f, VaF*^ = 0, (3.139) 
c 

where the first generally covariant equation expresses Group I of Max- 
well's equations, while the second generally covariant equation repre- 
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sents Group II. Herein, F*^ a = e^ aa ^ F a ^ is the pseudo-tensor which 
is dual to the electromagnetic field tensor F a p, while j M is the four- 
dimensional current vector. 

In terms of the chronometrically invariant formalism of General Rel- 
ativity, the generally covariant Maxwell equations (3.139) take the fol- 
lowing form (see Chapter 3 of the book [18] for details): 



*VjEi --H ik A ik =4wp 
c 

*V k H ik - ^F k W k -- 



* dE " n J\ 47r -i 



> I, (3.140) 



*V t H n - 
*V k E* ik 



E* ik A 



■F k E 



*ik 



*dH* 



+ DH* 



> II. (3.141) 



Herein, E* lk = — e lkn E k is the pseudo-tensor which is dual to the electric 
strength tensor Ei, H* 1 = | e lmn H mn is the pseudo- vector which is dual 
to the magnetic strength tensor H mn , while D — h lk D ik is the rate of 
space deformation. (See these definitions in formula (3.105) of §3.5). 

The physically observable charge density p and the physically ob- 
servable current vector f are the respective chr.inv. -projections of the 
four-dimensional current vector j M 



1 

P= — 

c . 



Jo 



'goo 



(3.142) 



Because the space of a liquid sphere under consideration is station- 
ary (the space metric does not depend on time), and also because we 
assume that the electromagnetic field is stationary, all terms containing 
the space deformation tensor Dn~ and the time derivatives of the elec- 
tromagnetic field strengths vanish. In this particular case, the chr.inv.- 
Maxwell equations (3.140-3.141) take the simplified form 



1 



H lk A lh 



4-7T p 



*V k W k - \F k H ik = —f 
c 



1 



E* ik A lk = 0 



II. 



*W k E* lk - —F k E* lk = 0 



(3.143) 



(3.144) 
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Substitute, into the equations (3.143-3.144), the formulae for the 
gravitational inertial force Fj (3.41) and the angular velocity of space 
rotation (3.48-3.51) we have obtained in the framework of the space 
metric of a regular rotating star (3.28). Also, substitute the formulae 
for the electric strength E % and the magnetic strength H* 1 (3.107) of 
the stationary vortex-free electromagnetic field we have obtained in the 
space. To simplify the algebra, take into account the following terms for 
the case if they appear in the brackets commonly with the others. For 
the Sun (w© ~ 2.7 xlCT 6 sec" 1 , r s ©=2.9xl0 5 cm, a© = 7.0xl0 10 cm), 
these terms take the following numerical values: 

^. = 4.1xl0- 6 
a 

^| = 8.5 xlO- 28 cm" 2 ) . (3.145) 
uj 2 

— = 8.1 xlO" 33 cm" 2 
c z 

For other regular stars, these terms take similar numerical values to 
within a few orders of magnitude. 

As a result, after some algebra we see that the equations of Group II 
of the chr.inv. -Maxwell-equations (3.143-3.144) vanish, while the equa- 
tions of Group I take the following form: 

-^L=4np, - = -2n 3 3 . 3.146 

The obtained Maxwell equations characterize the electromagnetic field 
that originates due to the charges and currents (the field sources). 
Should the right-hand side of the equations be zero, it would be a source- 
free electromagnetic field (existing independently of the sources). 

The field sources (the charge density p and the current vector j % ) 
are connected to each other through the law of conservation of electric 
charge. The law has the following general covariant formulation: 

V a f = 0. (3.147) 

The generally covariant law of conservation of the electric charge (3.147) 
is also known as the continuity equation. It means that the charge den- 
sity p and the currents j l , which are two chr.inv. -observable projections 
of the four-dimensional current vector j a are conserved within the four- 
dimensional volume of the field. Maxwell's equations are connected 
by the generally covariant Lorentz condition 

V ff A ff = 0, (3.148) 
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which means that the scalar electromagnetic potential ip and the vec- 
torial electromagnetic potential q l (the chr.inv. -observable projections 
of the four-dimensional potential A a of the electromagnetic field) are 
conserved within the four-dimensional volume of the field. 

In the general case of an arbitrary electromagnetic field, the conser- 
vation law V CT f = 0 (3.147) and the Lorcntz condition V CT A a = 0 (3.148) 
have the following chronometrically invariant formulation: 



*-^+pD + *V i j i -^F i j i = 0, (3.149) 



-*-^T + -D + *V i q i -irF i q i = 0. (3.150) 
c at c c l 

See Chapter 3 of the book [18] for details. The chr.inv. -differential 
operators *Vi = *Vi — ^Fi and *Vi can be found in Notations. 

Recall that we are considering a stationary vortex-free electromag- 
netic field. This means that the conditions (3.106) should be true 

*dif *dq n *dq t *dq k 

9^ = 0 ' ^r = 0 ' q * k = d^-^ = °- (3 - 151) 

In this particular case, and in the space which is free of deformation 
(Dik = 0), the conservation equation (3.149) and the Lorentz condition 
(3.149) are satisfied as identities. 

§3.7 Solving Maxwell's equations in the vortical electromag- 
netic field of a regular rotating star 

Consider a regular rotating star whose electromagnetic field is vortical. 
This means that the curl q ik of the three-dimensional vectorial chr.inv. - 
potential qi of the field is non-zero 

*dqi *dq k 

Assume, according to in Chapter 3 of the book [18] where we considered 
relativistic electrodynamics, that the four-dimensional electromagnetic 
field potential A a has the form 

dx a dx a dx? , 

A a = ^ — , 9ap =i. 3.153 

as as as 

In this case, the chr.inv. -projections of A a have the form 
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where 



V = 



<P 



v 2 = h ik vW k 



(3.155) 



is the relativistic scalar potential of the electromagnetic field [18]. We as- 
sume that the charges move within the star at small velocities (v 2 <C c 2 ). 
Then tp = tp. 

We also assume, according to our consideration of a rotating neutron 
star, that ip = const and q 1 =q 2 = 0. Then the chr.inv. -components of 
A a take the form 



p = const , 



c 



(3.156) 



where v 3 = while <fi is the equatorial coordinate of the spherical 
polar coordinates r, 9, (p. Assuming that the electromagnetic field curl 
is due to the angular rotation of the star, we have 

d<t> 



dr 



= U). 



(3.157) 



Then the non-zero components of the three-dimensional vector chr.inv.- 
potential and its curl q ik have the form 



o (put 
c 

q 3 = r z sm z 6 , 



931 = 



923 



c 

dq 3 = 
dr 
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2(fUJ . 2 

r sin 1 



2tpu 2 . 
r sin 0 cos ( 



(3.158) 
(3.159) 
(3.160) 
(3.161) 



In other words, we consider here only circular motion of the charge 
along the equatorial coordinate </>, which is the geographical longitude 
of the star (in spherical polar coordinates r, 6, (f>). 

Substitute now the formulae for An- (3.48-3.51) obtained for the 
space metric of a regular rotating star (3.28) into the definition of the 
magnetic strength tensor Hu~ (3.105). We obtain that the non-zero 
components of the tensor have the form 



#23 — — 



2ipujr 2 smO 



cos 9 



1 



3_Tg 

4a 



H. 



31 



2tpuir 



sin 2 9 — cos 9 



4a 



(3.162) 
(3.163) 
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Using the definition of the magnetic strength preudo- vector (3.105) 



we have, finally, 
H *i = _2^cj " 



1 ™^ rj _j_ Aran _ ^ r\*i 

JA in. n — ™ & Qrrin * & i 



1 - 



H *2 = 2^ 



2 a 3 

.,•2 



1 



2 a 3 



ostf+i 
sin 9 — cot ( 



1 



3r g 


3r g r 2 


4a 


2 a 3 


3r g 


3r s r 


4a 


2 a 3 



(3.164) 

(3.165) 
(3.166) 



while the covariant (lower-index) versions of the preudo- vector can be 
calculated as = hnH* 1 and = h 2 2H* 2 . In the framework of 
our approximation (with the higher-order terms withheld), we have 



H* 1 = - 



cos 9 + ^ 



H „ 2 = 2(pcj 
cr 



{sm9 - cot 9) . 



(3.167) 
(3.168) 



Concerning the chr.inv. -Maxwell equations: in the stationary elec- 
tromagnetic field, they have the form (3.143-3.144) 



1 



H tk A lk = 4np 



*V k H ik -\F k H ik = —f 
c z c 



\ I, 



1 



E* lk A, 



II. 



¥ V k E* ik - — F k E* lk = 0 



(3.169) 



(3.170) 



Substitute, into the equations, the obtained formula for the magnetic 
field strength, and also the formula for the electric field strength (3.105), 
which holds for the stationary electromagnetic field, has the form 



ikn : 



c z 2 a 3 



Equations of Group II (3.170) are satisfied identically. Equations of 
Group I (3.170) take the form 



3(pr g 
2 a 3 



o;</3cot# „ 
= 4*7>> L- a =-^3 ■ 



(3.172) 
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where p and j 3 are the charge density and the current of the vortical 
electromagnetic field. 

As is easy to see, these solutions are identical to the solutions (3.146) 
obtained in the vortex-free electromagnetic field. It is easy to obtain 
that the conservation equation V a j a = 0 (3.147) and the Lorentz condi- 
tion V a A a = 0 (3.148), whose chronometrically invariant formulations 
are (3.149) and (3.149) respectively, are satisfied as identities. 

This means that all the results obtained earlier in the vortex-free 
electromagnetic field of a regular rotating star are as well true in the 
present case where the electromagnetic field of the star is vortical. 

This happens because all the terms that appear in the equations due 
to the electromagnetic field curl vanish in the framework of the second- 
order approximation. In the parlance of physics proper, this means that 
the presence of a curl in the electromagnetic field does not change the 
field sources of a regular rotating star. The vortical electromagnetic field 
can be meaningful only in the case of exotic stars, whose characteristics 
differ from those of regular stars. We will see the difference in Chapter 5 
when considering rapidly rotating neutron stars (pulsars). 

§3.8 Geometrization of the electromagnetic field for a regular 
rotating star 

Using the geometric formula for the scalar electromagnetic potential 

(3.173) 

that follows from (3.119), we write down the non- vanishing chr.inv.- 
Maxwell equations (3.146), or (3.172) which is the same as 

3c 2 r 0 pn 3c 2 pn r a ,„„„, N 

p=p =^h = ^hi' (3A74) 

, ujc 2 pn cot 9 ,„ . 

3 = 3 =-2^Vg^' (3 - 175) 

3=3 = « = V^ff = ^ , (3-176) 

where p is the charge density and j 3 is the current of the vortical elec- 
tromagnetic field, while p and j 3 imply the vortcx-frcc field. The di- 
mcnsionless numerical coefficient n——^- (3.119) is within the range 
0 < n < j. To see why j, see formula (3.118). 
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The electromagnetic field sources are expressed herein through only 
the geometric characteristics of the star's space and the fundamental 
constants. This means that we have completely geometrized the sources 
of the stationary electromagnetic field of a regular rotating star. 

Now, express the electric and magnetic field strengths through the 
geometric formula for the scalar electromagnetic potential <p (3.173). 
Using the formulae for the non-zero components E\, H* , and H* 2 
obtained in the stationary electromagnetic field of a regular rotating 
star, we thus obtain 



H* 1 = -2luc^-(cos9 + ^J =-2c^(u J cos6 + n* 1 ), (3.178) 
H** = ^£ ( sin(? _ cot e) = 2c ^1 - fT 2 ) , (3.179) 

H*i = huH* 1 = H* 1 , (3.180) 

H t2 = h 22 H* 2 =r 2 H* 2 . (3.181) 

Here, according to our calculation (3.109-3.110) made in the framework 
of the space metric of a regular rotating star (3.28), 

ST 1 = = , n* 2 = UJC ° td , 0, 2 = cor cotO, (3.182) 
2 r 

n^n* j = lo 2 Q +cot 2 6^ . (3.183) 

Hence forth, we express the field density p em and the flow of momen- 
tum Jg m of the vortical electromagnetic field, which are the physically 
observable projections of the energy-momentum tensor of the field (see 
Einstein's equations). Using their general formulation made for any 
electromagnetic field, (3.102) and (3.103), we obtain 



2irG \Ac 2 1 *j j 2ir G 

2 

= Pcm + h 5 ' (3 - 184) 

-3 c 2 r g r nf uj sin6> 2 \ _ 3 c 2 r g uj n 
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As we can see, the observable characteristics of the electromagnetic 
field are expressed herein through only the geometric characteristics of 
the star's space and the fundamental constants. This concerns both 
the vortical field and the vortex-free field. In the sense of mathematics, 
this means that the electromagnetic field of a regular rotating star is 
completely geometrizcd. 

So, in the case of a regular rotating star, both Maxwell's equations 
and Einstein's field equations are satisfied with the inclusion of the elec- 
tromagnetic field. This means that they consist a self-consistent system 
of the Einstein-Maxwell equations that completely describes both grav- 
itational and electromagnetic phenomena inside regular rotating stars. 

Finally, we can conclude something pretty interesting for astro- 
physics by first writing the formula for the charge density p (3.174) 
in the form 

P=-M^^G~- ( 3 - 186 ) 

The first multiplier herein coincides with the formula for the "critical 
density" of substance in the Universe 

known from observational cosmology. Herein, H = — is the Hubble con- 
stant, while a is the radius of the observable Universe. In analogy to the 
Universe, the critical density can be formally introduced for any liquid 
star. Thus, we can express the charge density p of the electromagnetic 
field of the star as 

p = p CT yfn~G ^, (3.188) 

where n< \ and, numerically, -v/G = 2.6 xlO~ 4 cm 3 / 2 /gram 1 / 2 sec. 

If the charge density is p = p cr y/nG, the physical radius of the star 
coincides with the Hilbcrt radius a = r g . Because r g <C a for regular 
stars, we conclude that the charge density of the electromagnetic field 
inside any regular rotating star is much less than p cr y/nG, i.e. 

p< p cr y/nG. (3.189) 

A few words should be said at the end. As is known, the General 
Theory of Relativity is the geometric theory of space-time-matter. Its 
primary task is to express all physical phenomena as the manifestations 
of space (space-time) geometry. The gravitational field was initially 
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geometrized by Einstein, due to Einstein's field equations. However, 
the electromagnetic field was not geometrized: as was shown by Ein- 
stein, mathematically this problem in a general case is very non-trivial. 
Nevertheless, it is possible to solve this problem in a particular case 
where some particular conditions simplify the mathematics. Thus, as 
was shown above, we have completely geometrized the electromagnetic 
field in the internal field of a regular rotating star. 

§3.9 Conclusion 

This Chapter is complementary to the previous Chapter 2, wherein we 
considered regular stars including the Sun. Three primary tasks were 
achieved in this Chapter. 

First. In Chapter 2, when considering the internal space metric of 
a liquid star, we followed the historical path as Schwarzschild did when 
introducing the metric. Namely, — even when we introduced the com- 
plete form of the internal space metric of a liquid sphere (which contains 
singularities), we used the Schwarzschild notation. This notation comes 
from the general form of a spherically symmetric metric, and thus con- 
tains the coefficients e v and e A which are the functions of r and t. This 
is the common method for writing any spherically symmetric metric. 
But when we calculate the physically observable characteristics directly 
on the basis of the metric, we obtain them in the firm expressed through 
e v and e x which are unknown. Therefore, we obtain the physically ob- 
servable characteristics of the space in the incomplete form that needs 
further calculation of the coefficients e v and e A . This makes a huge 
additional trouble when solving particular problems in the framework 
of the space metric. Therefore, in this Chapter, we initially introduced 
the internal space metric of a liquid sphere in the final form, where the 
coefficients e v and e A are already expressed through the main charac- 
teristics of the sphere such as its physical radius and the Hilbert radius, 
and through the radial coordinate r and time t. As a result, we have 
obtained all the components of the fundamental metric tensor in ex- 
plicit form, without unknown coefficients. It was the subject of §3.1 
and §3.2. Therefore, once we (or someone else) further solve problems 
in the framework of the internal space metric of a regular liquid star, 
we initially will have implicit formulae for all the physically observable 
characteristics of its internal space. 

Second. We considered the space metric of a non-rotating liquid 
sphere. Nevertheless we know that the majority of stars rotate. Most 
probably, even all stars rotate, but many of them rotate slowly so that 
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the Doppler splitting of the spectral lines that is due to the rotation 
cannot be registered by modern methods of spectroscopy. In any case, 
as a matter of fact if we target a liquid star possessing an electromag- 
netic field, we should consider the internal space metric of a rotating 
liquid sphere. This metric was introduced in §3.3, then we subsequently 
introduced Einstein's field equations and Maxwell's equations in the 
form satisfying the metric. We showed that the electric component of 
the field primarily originates due to the gravitational field of the star, 
while the magnetic field component is primarily due to the field rota- 
tion. Also, we found that the vortical character of the electromagnetic 
field does not play a significant role in regular rotating stars. 

Third. Concerning the most important achievement of this Chap- 
ter: in §3.8 we showed that, in the case of the internal space metric of a 
rotating liquid star, all the physically observable characteristics of the 
electromagnetic field are expressed through only the geometric charac- 
teristics of the star's space and the fundamental constants. This fact 
means that, in the internal field of a rotating liquid star, the electro- 
magnetic field is completely geometrized. 



Chapter 4 



Stellar Wind 



§4.1 Finding the escape velocity condition for a star 

A flow of the particles of the stellar substance is permanently erupted 
from the surface of any star. A fraction of the flow consists of so rapid 
particles that they leave the gravitational field of the star forever, for the 
outer cosmos, thus producing a stellar wind* In terms of our mathemat- 
ical theory of liquid stars, this means that the particles of the surface 
layer of the star are faster than the escape velocity for the star. 

Why do the particles of the stellar substance leave the surface of a 
star? Can this process be likened to the boiling of water in a kettle, 
or is it entirely something else? Finding the answer to this question 
constitutes our research task for this Chapter. 

To answer this question we should study the motion of the particles 
of the stellar substance inside a star. To do it, we shall first find the 
formula for the escape velocity, expressed through the components of 
the space metric of a liquid star. Then we shall deduce the equations 
of motion of the particles of the stellar substance inside the star. Thus 
we shall obtain the physical conditions under which the particles of the 
surface layer are faster than the escape velocity for the star. Afterwards, 
we will be able to solve the equations of motion of the particles of the 
stellar substance. 

The said escape velocity, known also as the second cosmic velocity 
Vh, is the velocity at which a test-particle can "leave", forever, the 
gravitational field of the massive body. 1 ' 

Let us assume that the particles of the stellar substance travel, radi- 
ally, from within of the star onto its surface. Let the particles reach the 
surface then leave the star, forever, for the outer cosmos, thus forming 
stellar wind. We therefore refer to the formula for the velocity of the 
particles of the stellar substance, which is expressed through the escape 
velocity for the star, to as the escape velocity condition. 

*Wolf-Rayet stars differ from all the regular stars only by an extremely huge 
stellar wind: the flow is so powerful that a Wolf-Rayet star loses a substantial part 
of its mass with the stellar wind. 

tThis is analogous to the first cosmic velocity, known also as the orbital velocity, 
which allows the test-particle to be orbiting the massive body without falling down 
onto its surface. 
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Thus, for a spherically symmetric body whose mass is M, the escape 
velocity at the distance r from its center is 



v n = ^— — . (4.1) 
This formula comes from the mass-point metric (1.1), 

ds 2 = (l - y ) c 2 dt 2 - - r 2 (d9 2 + sm 2 0 d<f> 2 ) , (4.2) 

r 

where r g — 2 2 M , while M is the mass of the body (the field's source). 

As was shown in Chapter 2, the field of any liquid star has two pri- 
mary regions. They are described by two different metrics. The metric 
of a liquid sphere is valid from the center of the star (r = 0) to its surface 
(r = a). The mass-point metric is valid from the surface of the star to 
the outer cosmos. In other words, the particles of the stellar substance 
travel inside the star along those trajectories which are in accordance 
with the metric of a liquid sphere. If the particles leave the star (in 
the case that their velocity exceeds the escape velocity), they travel in 
the cosmos along those trajectories which are in accordance with the 
mass-point metric. 

Therefore, the velocity of the particles of the stellar substance trav- 
elling from the surface of the star for the outer cosmos results as the 
solution of the equations of motion of a mass-bearing particle accord- 
ing to the mass-point metric. Being expressed in terms of the escape 
velocity, it is the escape velocity condition for the star. 

We derive this formula as a solution of the chr.inv. -equations of non- 
isotropic geodesies [18, 19] 



dm m - m ^ , k 

^FiY 1 + —D lk v\ k = 0 



d( 



mv 



dr 



2m(D l k + Al) v k - mF l + mA^ v"v fc = 0 



(4.3) 



which are the equations of observable motion of a mass-bearing particle 
travelling at the observable velocity v\ The equations result as the 
observable projections of the well-known generally covariant equations 
of non-isotropic geodesies (see [18, 19] for detail). 

We solve the equations (4.3) for a particle of the stellar substance, 
which travels only along the radial direction r. Therefore, 

^ = ^0, v 2 =v 3 = 0. (4.4) 
dr 
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To solve these equations (4.3), we need to formulate the charac- 
teristics of the space of the mass-point metric (4.2). As is seen from 
the metric (4.2), the space is free of rotation and deformation (Aik = 0, 
Dik = 0)- Only the gravitational inertial force Fi and the Christoffel 
symbols A l nk remain non-zero. Calculating these quantities, and also 
the components of the chr.inv. -metric tensor hik according to their def- 
initions given in §1.3, we obtain, for the metric (4.2), 



Fi = - 



2r 2 



F = — 



c'r g 
2r 2 



(4.5) 



hn = j^ = l- r -f, h 22 = L.=r\ h 33 = -L § =r 2 S m 2 9, (4.6) 



r g 1 



2r 2 



A 1 

A l _ A 33 
^22 — . 2 a 
Sin 9 



(4.7) 



A 2 - A 3 - 

^12 — ^13 — 



A 2 

^33 



— sin 9 cos 9, Al 3 = cot9. (4.8) 



With these, we obtain that the chr.inv. -equations of motion (4.3) in 
the space of the mass-point metric have the form 



1 dm 
to dr 

LA 

to dr 



2r 2 
dr\ 



1 - 



dr 
1 



2r 2 



1 - 



dr) 2r 2 



(4.9) 



where 



m 



mo 



dr 
d^ 



(4.10) 



Here, denoting the relativistic mass of the particle on the surface (r = a) 
of the star as TO( 0 ) (this is the "start-mass" of the particle when leaving 
the star), and denoting the observable velocity of the particle when it 
leaves the star as r 0 , we have 
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Proceed to solve the chr.inv. -equations of motion (4.9). Substituting 
the scalar equation into the vectorial equation, we obtain the vectorial 
equation of motion written with respect to the radial distance r 

• 2 2 

+ Vf =0- (4-12) 

r 

Denote r = y, then 

/ , dy 

r = yy, y= Tr , 

and the equation (4.12) takes the form 

yy'-^^r + ^f = 0. (4.13) 

r 

Assuming u(r) = y 2 , we reduce it to the linear differential equation 



u ,_2^_«_ + ^ =0- (414) 

7" 1 _ l£ 



This equation has the following exact solution: 

u = c- F (u 0 + g{r)c F dr), u 0 = y 2 = r 2 , (4.15) 



where 



F(r)= f f(r)dr, f { r) = -L- , 9 {r) = ^. (4.16) 

Jr r 1 — r 

r 

Integrating the function / (r) , we obtain 

F(r)=ln|^), e' = | ^|, (4.17) 



1_ -JL 




a c 2 r„ (l - —\ dr 

' a) -^(l--)|l ^1. (4.18) 



Substituting (4.16-4.18) into (4.15), and neglecting the high powers 
of the term ^ (this ratio is tiny for regular stars), we obtain 

+( ?(H_H\ (4 . 19 ) 
a r I \ a r I 
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From here, we obtain the formula for the radial velocity of the particle 
of the stellar substance, which leaves the star with stellar wind. Since 
v„ (4.1) on the surface of the star (r = a) takes the form 



we obtain 



r=^=cJ r l±^-^+ 2 4(4- C ^-). (4.21) 
dr W c l r c z \c z r ) 

This is the escape velocity condition we were looking for. If r 0 = 0, the 
equation (4.21) manifests the obvious condition 



dr 
rf7 



'9 



< v„ . (4.22) 



According to this condition, the particle of the stellar substance cannot 
leave the gravitational field of the star, if its start- velocity on the surface 
of the star is zero. Therefore, when further considering stellar wind, we 
always assume f 0 ^ 0 in all the equations of the theory. 

Let us obtain the final simplification of the escape velocity condi- 
tion (4.21). Compare the estimated numerical values of all the terms 
contained in the radicand. We denote the last term of the radicand as 



2r n 2 



(4.23) 



For the Sun, which is a typical regular star, we have: v n = 617 km/sec, 
r g = 2.9 km, f 0 = 750 km/sec* and a = 7.0 xlO 5 km. Since q = 0 by r = a, 
assume Ooas for stellar wind. We obtain, after some algebra, 

* 2 2 

T ° + 9 V " ~ KT 5 , ^<4.1xl0- 6 , (Z <5.3xlO- 11 . (4.24) 

For a typical star of the Wolf-Rayet family (see Table 1.1), we have: 
v„ = 982 km/sec, r g = 150 km, r 0 = 2200 km/sec, and a = 1.4xl0 7 km. 
Therefore, for a typical Wolf-Rayet star, we obtain 

* 2 2 

r ° + 0 V " ~ 6.4 xl0~ 5 , I2.<l.lxl0- 5 , ,7<1.2xl0- 9 . (4.25) 
c 2 r 



*ro — 750 km/scc is typical for the particles of the fast component of the solar 
wind, whose composition is that of the photosphere. In contrast, the slow component 
of the solar wind has a composition close to that of the corona. Its particles travel 
from the Sun at a velocity of about 400 km/scc. 



102 



Chapter 4 Stellar Wind 



As is seen, the term q has so small a numerical value (four orders less 
than the other terms in the formula) that can be neglected for stellar 
wind, which comes from both a regular star and a Wolf-Rayet star. 
Therefore, the final formula for the escape velocity condition has the 
form 



As follows from the final formula, on the surface of the star (r = a) 
the velocity of the particle of the stellar substance is f 0 . 

§4.2 Light-like (massless) particles inside a regular star 

We now consider how particles of the stellar substance and particles of 
light behave inside the star. (Stars are filled not only with substance 
but also with light.) 

First, consider light-like (massless) particles inside a regular star. 
Such particles travel along isotropic geodesic lines. The chr.inv. -equa- 
tions of isotropic geodesies have the form [18, 19] 



These are the equations of observable motion of a light-like particle — 
a photon whose frequency is w, — which travels with the observable 
velocity of light c*. These chr.inv. -equations emerge as the observable 
projections of the well-known generally covariant equations of isotropic 
geodesies (see [18,19] for details). 

As previously, we assume that regular stars do not rotate or deform 
(Aih = 0, Djfc = 0). Also, we consider only a photon which travels radi- 
ally (along the direction x 1 = r) from the center of the star to its surface. 
Therefore, the isotropic geodesic equations (4.27) inside a regular star 
take the particular form 




(4.26) 



du> 




d(ujc l ) 
dr 



du 
~d^ 




0 



d^c 1 ) 
dr 



} , 



(4.28) 



- uF 1 +ujA\ 1 c 1 c 1 = 



0 



where the observable (light) velocity of the photon is c 1 = 
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Consider the scalar geodesic equation of (4.28). Substituting iq 
(3.6), obtained for the metric of a liquid sphere, we have 



1 dw 

lo dr 



dr 



r 2 dr 



(4.29) 



Re-write this equation in the following form, which can easily be inte- 
grated: 



d In lu = 



3\/l 



1 



r 9 r 
n3 



= din 



(4.30) 



We study only the photons inside the star. We therefore will look for 
the solution for the interval r g ^ r ^ a. We obtain, after integration, 



D 



uj = 



(4.31) 



where B is the integration constant. 

Assume that the photons originate from the Hilbert surface (r 0 = r g ). 
Then 



B = w 0 |3 A /l-^ 



(4.32) 



where ojq is the initial value of the photon's frequency (on the Hilbert 
surface of the star). Since r g <^a for regular stars, we neglect the high- 
power terms of ^ . Finally, the solution of the scalar geodesic equation, 
which manifests the photon's frequency (4.31), takes the form 



w 0 3 \ 1 



3a/1 - ^ 



(4.33) 



1 



Now, consider the vectorial geodesic equation of (4.28). With our 
assumption of the radial motion of the photon, it has the form 



dr 2 



1 duj dr 
ui dr dr 



+ - — - + AM 



(!) 



F 1 = 0. 



(4.34) 
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Denote r=£| and r=^. Substitute (4.29), (3.8), and F 1 

(3.7). Thus the vectorial geodesic equation (4.34) transforms into the 
non-linear differential equation of the second order with respect to r 

r„r f 2 
r h-—/ , , , + 



a3 (V^-^?)\A 



+ ¥ _^ + £V_j^- _ 0 . (4 . 35) 



a3 



In this form the equation is simply non-solvable. Therefore, we simplify 
it by the formula for f 2 taken from the obvious relation hikC l c k = c? , 
which in the present case has the form 



2 

= c 2 . (4.36) 



_ r g r- 
a 3 

As a result, the initial equation (4.35) takes the form 

r + = 0 . (4.37) 

This is the equation of harmonic oscillation with the frequency 



<l = -JUL = Z!L = J™M i (4.38) 
a V a a V a A 

which is obviously dependent on the escape velocity calculated for the 
star, v„ (4.20). 

In general, the frequency 0 (4.38) only depends on the mass M 
and radius a, which are the integral characteristics of the star. We 
therefore refer to it as the proper frequency of the star. Table 3.1 gives 
the numerical values of the proper frequency fi for the typical members 
of the known families of stars. 

The proper frequency reaches its ultimate-high magnitude S! max = ^ 
by r g = a. This is the case of gravitational collapsars (black holes), which 
is also applicable to the Universe as a whole. According to observational 
estimates, the Universe's radius is a — 1.3 xlO 28 cm, and it is the same 
as its Hilbcrt radius r g : the Universe is a huge gravitational collapsar. 
Therefore, the proper frequency of the Universe is 

fW = - = 2.3xl(T 18 sec" 1 (4.39) 
a 
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Object 


Mass M, 


Radius 


Proper frequency 




gram 


a, cm 


fl, sec -1 


Wolf- Ray et stars 


1.0 xlO 35 


1.4 xlO 12 


7.0xl0~ 5 


Red super-giant* 


4.0 xlO 34 


7.0 xlO 13 


1.6xl0~ 7 


White super-giant* 


3.4 xlO 34 


4.8 xlO 12 


6.4xl0~ 6 


Sun 


2.0 xlO 33 


7.0 xlO 10 


8.8xl0~ 4 


Jupiter (proto-star) 


1.9 xlO 30 


7.1 xlO 9 


8.4 xl0~ 4 


Red dwarfs 


6.7 xlO 32 


2.3 xlO 10 


2.7xl0~ 3 


Brown dwarf* 


4.1 xlO 31 


7.0 xlO 9 


7.4 xl0~ 2 


White dwarf § 


2.0 xlO 33 


6.4 xlO 8 


1.0 


Universe 


8.8 xlO 55 


1.3 xlO 28 


2.3xl0~ 18 



*Betelgeuse. 'Rigel. *Corot-Exo-3. § Sirius B. 



Table 3.1: The proper frequency fi for the typical members of the known 
families of stars, and for the Universe. 



which matches the numerical value of the Hubble constant, which is 
H = ^ = (2.3 ± 0.3) xl0~ 18 sec -1 . In this case, the integral mass of the 
Universe should be, according to (4.38), 

o2 a 3 

M = —— = 8.8 xlO 55 gram (4.40) 
2 G 

which coincides with the observed range of the average density of sub- 
stance in the Universe, which is from 10~ 28 to 10 -31 gram/cm 3 . 

Let us return to the vectorial equation of motion of the light-like 
particles inside liquid stars. 

The vectorial geodesic equation in its final form (4.37) solves as 

where B\ and B 2 are the integration constants. Assuming r and r at 
the initial moment of time r 0 = 0 to be r 0 = r g and r 0 = c , we obtain 

Bi=r g , B 2 = a.f^. (4.42) 

V r 9 

As a result, we obtain the final solution for r, which is 

/a civ 
— sin fir, SI = - * — , (4.43) 
r„ a V a 
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i.e. the harmonic oscillation equation r = A\ cosOr + A2 sin fir. Differ- 
entiating (4.43), we obtain the oscillation velocity of the photon 



As is seen from the solution (4.43), the light-like matter of each 
single star oscillates at the frequency fl (4.38) depending on the mass 
and radius of the star, and with two primary amplitudes: 

a) The oscillation with the amplitude A\ — r g . In fact this means 
that the surface of the Hilbert core of the star, wherein the stellar 
energy is released, oscillates at the proper frequency of the star. 
Thus the released light-like matter oscillates at the same frequency 
and amplitude as its source (the Hilbert core); 

b) The oscillation with the amplitude A 2 — y 7 a 3 /r g , which coincides 
with the outer space breaking of the star's field (see Chapter 2 for 
details). The space breaking is located outside the star, in the cos- 
mos. For the Sun (a = 7.0*10 10 cm, r g = 2.9><10 5 cm), we obtain 
A 2 = 3.4xl0 13 cm = 2.3 AU which is the distance from the Sun 
to the maximal concentration of asteroids (in the asteroid strip). 
This means that the spherical surface the outer space breaking of 
the field of each single star oscillates at the proper frequency of 
the star, as well as the Hilbert core. 

Hence, we arrive at the following fundamental conclusions: 

1. The surfaces of both the inner space breaking and the outer space 
breaking of the star's field oscillate at the same frequency ft, which 
is the proper frequency of the star; 

2. This frequency and the amplitudes thereof only depend on the 
mass M and radius a of the star; 

3. In fact, this common oscillation of light-like matter comprising the 
star is due to the gravitational field of the star (originating in the 
star's mass M). 

How does this oscillation affect the photon's frequency? To answer 
this question, consider the obtained solution for the photon's frequency 
to (4.33) in the framework of the two ultimate cases which correspond 
to two oscillation amplitudes: r = A\ = r g and r = A 2 = \Ja z jr g . Thus 
the frequency takes the following numerical values: 



r = c cos CIt — — sin Or 

a 




(4.44) 



r = A 1 =r g : 



UJ = UJq 




= U 0 , 



(4.45) 
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U! = OJq 



■4 



a 



i 



(4.46) 



a 



As is seen from these formulae, the primary oscillation of the gravita- 
tional field of the star does not affect photons which are close to the 
Hilbert core (at the center of the star). The oscillation affects only 
photons at large distances from the Hilbert core. 

§4.3 Particles of the stellar substance inside a regular star 

Such particles travel along non-isotropic geodesies. The chr.inv. -equa- 
tions of non-isotropic geodesies (see [18,19]) have the form (4.3): 



We assume a regular star to be a liquid sphere, which is free of rota- 
tion and deformation (Aik = 0, Dik = 0). For a particle of the stellar 
substance, which travels inside the star radially from the center to the 
surface, the observable velocity is v 1 = ^ while v 2 = v 3 = 0. In this case, 
the chr.inv. -equations of non-isotropic geodesies (4.47) take the form 



They have the same structure as the chr.inv. -equations of isotropic 
geodesies (4.28). They solve in the same way. But mass-bearing par- 
ticles do not possess the light speed condition hikc l c k ~c 2 (4.36) we 
have used for the isotropic geodesic equations. Therefore, the chr.inv. - 
equations of non-isotropic geodesies (4.48) will have another solution 
than that which we have obtained for the chr.inv. -equations of isotropic 
geodesies (4.28). 

Substitute, into the scalar equation of (4.48), Fi (3.6) which we have 
obtained for the metric of a liquid sphere. We thus obtain the scalar 
geodesic equation in the form 



dm 




d(mv l ) 




(4.48) 



108 



Chapter 4 Stellar Wind 



This equation can be re-written in the form 



dm TO : 



= din 



r 9 r 
n 3 



(4.50) 



which is easy to integrate. With r g < r < a (we study only the particles 
inside the star) we obtain, after integration, 



D 



to = 



3a/1-^ 



(4.51) 



1 



where B is the integration constant. 

Let the particles start from the Hilbert surface (r 0 =r g ). Then 



B = TO 



(0) 



where 



TO. 



(0) 



3, 



too 



1 



(4.52) 



(4.53) 



2 



is the initial value of the relativistic mass of the particle on the Hilbert 
surface of the star. Since r g <§C a for regular stars, we neglect the high- 
power terms of With all these taken into account, the solution of 
the scalar geodesic equation, which is (4.51), takes the form 



n»(0) 



3\/l 



(4.54) 
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where 

m (o) = I w ( 4 - 55 ) 

in the framework of our approximation mentioned above. 

Now, consider the vectorial geodesic equation of (4.48). With our 
assumption that the particle of the stellar substance travels radially, 
from the center of the star to its surface, the equation has the form 

d 2 r 1 dm dr A , / dr \ 2 . n 

d^ + ^d-r +A "{d- T )- F =°- ^ 

With (4.29), A^ (3.8), and F 1 (3.7) substituted, and with the 

notations f = and r = ^ , this equation transforms into the non- 
linear differential equation of the second order with respect to r 

.. _ rgr_ 

2 + ilf , ^ Z ^ == = o . (4.57) 



a 3 



1 _ r _a]_ a 3 „ f ~ f. 7~?> 

1 ^ 3 y 1 ---y 1 - — 



It is identical to the equation (4.35) we have obtained for the photon, 
and is non-solvable as well. To simplify the equation, we express f 2 
from the obvious relation hnrr — r 2 , which takes the form 



where 



mo 

m = 



(4.59) 



It follows, from (4.54), that 



toq 

TO 



(4.60) 
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Therefore, from (4.58), we obtain 



•2 2 

r = c 



1-- 



(3^¥-^¥) 2 (i-l 



(4.61) 



Substituting this formula for r 2 into the initial differential equation 
(4.57), and neglecting the high-power terms of we obtain the vecto- 
rial geodesic equation (4.57) in the solvable form 



„ | (c 2 + fl)r g r _ q 
2 a 3 



(4.62) 



This is the equation of harmonic oscillation with the frequency 



(c 2 + fp r g 
2 a 3 



(4.63) 



It concerns the particles of the stellar substance. As is easy to see, this 
formula transforms into the formula for the photon's frequency il (4.38) 
by the ultimate condition r = c. 

The vectorial geodesic equation (4.62) solves as 



r = Qi cos fir + Q2 sin Or, 



(4.64) 



where the integration constant Q\ and Q2 results from the conditions 
ro = r g and fg at the initial moment of time tq = 0. We obtain 



Qi 



•a i 



Q2 



foCtv2a 



Therefore, the final solution for r has the form 

r 0 a^/2a 



r g cos fir + 

v (c 2 + ^)r s 



sin O7 



(4.65) 



(4.66) 



which is the harmonic oscillation equation r = A\ cosOt + A2 sinfir. 
Differentiating (4.66), we obtain the velocity of the particle 



5 sinSZr + r 0 cos ilr . 

2a A 



(4.67) 



The obtained solution (4.66) manifests that the liquid substance of 
each single star oscillates at the frequency VL (4.63) depending on the 
mass and radius of the star, and with two primary amplitudes: 
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a) The amplitude Ax = r g , which is the same as that for photons (see 
above). The Hilbert core of the star, consisting of liquid substance, 
oscillates at the proper frequency of the star; 

b) The amplitude 



„ . (4.68) 

V( C + r o) r 9 

It depends on the initial velocity of the particles of the stellar 
substance, fo- If r o = c j A 2 = y a 3 /r s , which coincides with the 
amplitude of stellar photons (see above). According to (4.68), 
the initial velocity for the particles whose oscillation amplitude 
reaches the surface of the star (A 2 = a) is 



C V^ yn v„ 2GM 
ro = 7= = , ^ -7= , Vn = \ , (4.69) 



y / 2a- r g ^2 - iL ' 



where v„ (4.20) is the escape velocity for the star (with which the 
particle leaves the gravitational field of the star forever). Thus, 
by the condition A 2 ^ a in (4.68), we obtain the velocity which 
is necessary for a particle of the stellar substance in order to be 
erupted from the surface of the star 



ro>^. (4.70) 

Let us transform the obtained formula for the proper frequency fl 
(4.63) in order to express it through the orbital velocity v L for the star* 



^ = -\hr\ l + ^^^F\ l + ^ = -\ l + ^ • ( 4 - 71 ) 

a V 2a V c 2 av^ V c 2 a V c 2 v 7 
Using this formula, we express r (4.66) in the form 

r = r g cosflT-\ T ° a sin Or, (4.72) 



which is r = Ai cos Or + A2 sin Or. So we have 



A 1= r g , A 2 = . (4.73) 

1 + ^ 



*The orbital velocity vi, known also as the first cosmic velocity, allows the test- 
particle to be orbiting the massive body without falling down onto its surface. 
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Thus, r (4.67) transforms to 



r = - — \/l + ^ sinOr + r 0 cosftT. (4.74) 
a V 

Consider now the amplitude A 2 (4.73) for some particular cases, 
where it has different numerical values: 

1) If r 0 = 0, we have A 2 = 0 according to the definition of A 2 (4.73). 
Hence the particles of the stellar substance oscillate at the ampli- 
tude r g . In other words, if r 0 = 0, the particles cannot leave the 
surface of the star; 

2) If fo = v I , the particles also cannot leave the star. This is because 
they oscillate at an amplitude less than the physical radius of the 
star 

A 2 = < a ■ (4.75) 

3) If ro= v ii> then the particles leave the star. This is because in the 
case of fo = v n we have 

A 2 = , aV ^ = ~ (l- ^r) aV2~aV2> a; (4.76) 



2 c 



4) If A 2 = a: the amplitude equals the physical radius of the star. 
Then we obtain, from the definition of A 2 (4.73), that 



fo = ^=c(l-^)v I <v I . (4.77) 



1 + 

The particles are a little slower than the orbital velocity for the 
star. This means that, if the amplitude equals the physical radius 
of the star (A 2 = a), the particles may jump up from the surface 
of the star and yet they do not leave the star for the orbit (they 
always fall back down on the star). 

Thus, the new mathematical theory of liquid stars provides a solid 
theoretical ground to stellar wind as that consisting of two components. 
One is a little slower than the orbital velocity for the star, while the 
other travels faster than the escape velocity. This conclusion matches 
observational data. For example, the solar wind consists of two compo- 
nents. The slow solar wind travels at a velocity of about 400 km/sec that 
is slower than the orbital velocity for the Sun, which is v : = 440 km/sec. 
The fast solar wind travels with a velocity of about 750 km/sec, which 
exceeds the escape velocity, v„ = 617 km/sec. 
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§4.4 Conclusion 

Let us summarize the main conclusions which we have obtained on the 
origin of stellar wind. The conclusions are as follows: 

1. Given any star, the light-like matter that fills it oscillates at the 
frequency 



Each single star has its own frequency fl according to its mass M 
and radius a. Therefore, it is the proper frequency of the star; 

2. The oscillation occurs with two primary amplitudes. The am- 
plitude A\ — r g coincides with the surface of the Hilbert core of 
the star, wherein the stellar energy is released. The amplitude 
A 2 = \J o? ' /fg coincides with the surface of the outer space break- 
ing of the star's field, which is located in the cosmos. For the Sun, 
A 2 = 3.4 x 10 13 cm = 2.3 AU coincides with the maximal concentra- 
tion of asteroids (in the asteroid strip). This common oscillation 
of the light-like matter of the star is due to the gravitational field 
of the star (its source is the star's mass M). In other words, it is 
the own "breathing" of the star; 

3. Particles of the stellar substance oscillate at the frequency 



and with two primary amplitudes. The frequency depends on the 
initial velocity of the particles, ro, and can be expressed through 
the escape velocity v„ and the orbital velocity v : for the star; 

4. The amplitude A\ = r g is the same as that for photons. This 
means that the physical surface of the Hilbert core oscillates at 
the frequency 17 that above. The other amplitude A 2 depends on 
the initial velocity of the particles 



5. Stars emit light (photons) and erupt the particles of the stellar 
substance (stellar wind) not by the order of special physical condi- 
tions, but automatically. The equations of motion of the particles 
(both the photons and the substance), which travel radially inside 





A 2 = 
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a liquid star, are the harmonic free-oscillation equation 
r + il 2 r = 0, fl 2 = -- 



2F 1 c'r q 



where F\ = — c ^J i r is the linearized form (in the sense of r g <C a) 
of the force of gravity acting inside any liquid star. This is a non- 
Newtonian gravitation (the force is proportional to the distance), 
which is the cause of the oscillation of both the stellar light-like 
matter and the stellar substance. Once the oscillation amplitude 
exceeds the physical radius of the star, the particles come out the 
star for the cosmos. Therefore, the cause of the emission of stars 
is the internal structure of these self-gravitating bodies, which are 
the liquid spheres in the weightless state in the cosmos; 

6. According to the theory, stellar wind should consist of two com- 
ponents: slow stellar wind and fast stellar wind. Particles whose 
oscillation amplitude reaches the star's surface (A 2 = a) have the 
initial velocity 

ro = ; = 1 - tto v, < v, . 



1 + ^ 



2c 2 



which does not exceed the orbital velocity Vj for the star. Particles 
which are as fast as the escape velocity for the star (r 0 =v n ) have 
the oscillation amplitude 

c 2 

This means that slow stellar wind is composed of particles whose 
oscillation amplitude is in the range of a ^ A 2 < a These par- 
ticles leave the surface of the star, but not forever. They always 
fall back down on the star. Fast stellar wind is composed of par- 
ticles whose oscillation amplitude is A 2 ^a Particles of fast 
stellar wind leave the gravitational held of the star, forever, for the 
outer cosmos. Naturally, solar wind consists of slow solar wind, 
which travels at ^400 km/sec (slower than v IQ = 440 km/sec), 
and of fast solar wind travelling at ~ 750 km/sec (faster than 
Vno = 617 km/sec). 



Chapter 5 



Neutron Stars and Pulsars 



§5.1 Introducing the space metric of a rotating neutron star 

This Chapter is most short, and most complicate in the math among 
the other Chapters of this book. We will apply our model of liquid stars 
to neutron stars and pulsars. The high level of complexity is due to the 
fact that, once we introduce, in the space metric, the rotation around 
even a single coordinate axis, the further calculations become highly 
problematic. Anyhow, let us begin. 

Neutron stars and pulsars are attributed to Type II of our classifica- 
tion of stars according to the General Theory of Relativity (see Table 1.1 
in §1.2). This means that the physical radius a of such a star is a little 
larger than its Hilbert radius r g : the star is almost a collapsar, but still 
has a possibility to shine as a regular star. In §1.2 we showed that the 
space metric of a liquid sphere transforms into de Sitter's metric of a 
vacuum sphere under the condition of gravitational collapse (a = r g , i.e. 
the liquid sphere is a collapsar). The metric has the form (1.16) 



The physical parameters of neutron stars and pulsars are close to the 
parameters of collapsars, but are not the same (see Table 1.1). Therefore 
the metric (5.1), which includes the collapse condition, is close to the 
true metric of a neutron star or a pulsar, but is not. 

How to modify the space metric of a collapsar, (5.1), in order to 
obtain the metric of a neutron star or a pulsar? To leave the collapse 
condition, but be near it in the same time. Easy. 

Remind, the particular condition of gravitational collapse (goo = 0) 
comes from the general condition of gravitational collapse according to 
which the physical observable time r (1.30) stops on the surface of the 
object 



If the local space of the object does not rotate (all goi = 0), the afore- 
mentioned particular condition of collapse (goo = 0) occurs. 




(5.1) 




(5.2) 
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Once the object rotates (at least one of all three quantities goi is non- 
zero), the condition 500 = 0 can remain true on the surface of the object 
but does not mean gravitational collapse. This is due to the second 
term of the complete condition of collapse (5.2) which is non-zero in this 
case. Therefore, once the rotation is introduced into the metric (5.1), 
the metric describes the sphere which is out the state of gravitational 
collapse. The faster the sphere rotates, the more its state differs from 
the state of a collapsed sphere. 

If we will find, in addition to the modified metric which contains the 
rotation, Einstein's field equations in the form containing the strong 
magnetic field and also, in the same time, satisfying this metric, we will 
have the complete description of a neutron star or a pulsar. This is our 
research plan for this Chapter. 

First, we add the space rotation to the metric (5.1), according to the 
theory of chronometric invariant: see formulae (1.45) of §1.3. Assume 
that the object — the liquid sphere of the radius a — rotates, with an 
angular speed oj, along its equatorial axis (the axis <f> in the spherical 
coordinates r, 9, <j>). In this case, the initially metric of the collapsed 
liquid sphere (5.1) takes the following form 

,2 1 r 2 \ 0,2 2wr 2 cos8 , ,, 

ds 2 = - 1 ^ c 2 dt 2 + cdt d6 - 

4 V a 2 J c 

dr 2 

- r 2 (d0 2 + sin 2 0 d<j> 2 ) , (5.3) 



which means that the sphere is not a collapsar (due to the rotation, see 
the explanation above). 

The linear velocity of such a rotation is determined by goi of the 
space metric according to the general formula (1.45). In the present 
case of the metric (5.3), it has the form 

2uoar 2 cos 6* ,_ A . 

v 1 =v 2 =0, v 3 = - - . (5.4) 



The ultimate magnitude of the rotation speed of the neutron stars, 
registered in the astronomical observations, is about 1, 000 km/sec. We 
therefore neglect the terms p-, where v 2 ~h lk ViVk <C c 2 . The condition 
v 2 -C c 2 also means that the rotating body remains to be a sphere. 

The three-dimensional observable chr.inv. -metric tensor hik (1-34) 
of the space, whose metric is (5.3), has the components 

hl1 = ]~n = ' h22 = ]^22= r2 > ^3-^=r 2 sin 2 0, (5.5) 
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while the determinant of the chr.inv. -metric tensor hik, and its non-zero 
logarithmic derivatives along the spatial coordinates are 

, ii , a 2 r A sin 2 # . . 

/i = det ||/Nfc|| = — o 5-' ( 5 - 6 ) 



*51nv^ 2 a 2 



dr r (a 2 — r 2 ) 



(5.7) 
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= cot 9. (5.8) 



Also, due to the assumed condition v 2 -C c 2 (the non-relativistic ro- 
tation of the object), the chr.inv. -differential operator along the spatial 
coordinates (1.41) coincides with the regular differential operator. 

With goo an d goi of the metric (5.3), we now deduce the formulae 
for the chr.inv. -vector of the gravitational inertial force Fj acting in the 
space, and for the chr.inv. -tensor of the angular rotation of the space, 
Aik . According to the definitions of these quantities which come from 
the chronometrically invariant formalism (see §1.3), we obtain 

2 2 

Fi = - 2 j , F* = — , (5.9) 

2wa 3 r cosd 2u>acos9 



,13 



(a 2_ r 2)3/2' rv^^sin 2 ^ 



V 'a 2 — r 2 r 2 \J a 2 — 



sm t 



(5.10) 



Two other chronometrically invariant (physical observable) quanti- 
ties will be needed for our further calculations of Einstein's field equa- 
tions. These are the chr.inv. -Christoffcl symbols of the second kind A l kn 
and the chr.inv. -curvature tensor Cjfcy. After some algebra according to 
the general formulae of these quantities (see §1.3 for detail), we obtain 
that, in the space of the metric (5.3), the chr.inv. -Christoffcl symbols 
A\. n have the following non-zero components 

A i r Ai r(a 2 -r 2 ) 

A ll = ^2— 2- A22 = -2 ( 5 - U ) 

^— ^r 2 ^ , A 2 2 = A?3 = i, (5.12) 
A§ 3 = -sin0 cos 9, A^ 3 = cot^, (5.13) 
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while the non-zero components of the chr.inv. -curvature tensor Cikij are 

r 2 

C1212 = — 2 o . ( 5 - 14 ) 



r 



2 



Ci3i3 = ~5 n-sin^, (5.15) 



r 



4 



C 23 23 = — ? sin 2 0. (5.16) 



Respectively, the contraction d k — h mn Ci m kn (the chr.inv. -analogy of 
Ricci's tensor) has the following non-zero components 

2 2r 2 2r 2 
Cn=— 5 2, C 22 = r , C 33 = r sin 2 0. (5.17) 

With these characteristics of the space of the metric (5.3), we are 
able to deduce Einstein's field equations in the form satisfying the metric 
(that is the next step of our research of neutron stars and pulsars). 

§5.2 Einstein's field equations and the conservation law equa- 
tions satisfying the metric 

Consider the chr.inv. -Einstein equations in the general form (1.85-1.87). 
In a stationary space (that means that the space is free of deformations) , 
such as the space of the metric (5.3) that we suggest to neutron stars 
and pulsars, the chr.inv. -Einstein equations take the simplified form 

A jt AV + (*V, - 1 Fj^j F==-^ (pc 2 + U)+ Ac 2 , (5.18) 

-| FjA ij - *Vj A 11 = xJ\ (5.19) 

2A ij Al + I (*V,F fe + *V fe F,) - 1 F t F k - c 2 C tk = 

= 77 {pc 2 h lk + 2U ik - Uh ik ) + Xc 2 h ik . (5.20) 

Herein, p, J 1 , and U lk are the chr.inv. -projections (1.84) of the energy- 
momentum tensor of the continuous matter that fills the space. These 
are the observable density of mass, the observable density of momentum, 
and the observable stress-tensor, respectively. While U = h lk Ui k is the 
trace of the observable stress-tensor. Note that the energy-momentum 
tensor has now an arbitrary form. So, the sort of the distributed matter 
is not specified for yet. 
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Substitute, into the equations, the chr.inv. -characteristics of the met- 
ric (5.3). While doing so we should take into account the fact that the 
initially (non-rotating) metric (5.1) was deduced under the conditions 
a 2 = j > 0 and A > 0 (see §1.2 for detail). As a result, we transform the 
chr.inv. -Einstein equations (5.18-5.20) to the form 

8u; 2 a 4 cot 2 6» 2u 2 a 2 x, , „. 
(a2 _ r2)2 +#Z?=20"? + U), (5-21) 

2 / UJaCOt6 =-xj\ (5.22) 
r 2 \fa 2 —r 2 sin (9 

8w 2 a 4 cot 2 6> 2 rn _ xt/u(a 2 -r 2 ) 



(a 2 -r 2 ) 2 2 
4w 2 a 4 r cot 



(a 2 -r 2 ) 2 



= -xU 12 , (5.24) 



2^ 2 a 2 x 2 xE/gg 

2c^ 2 a 2 8^ 2 a 4 cot 2 fl x 2 xl/33 
a 2 -r 2 (a 2 -H) 2 2 r 2 sin 0 

By the calculation of U = h ik U lk = h 11 Un + h 22 U 2 2 + h 33 U 33 from 
the three respective tensorial equations of these, we obtain the relation 
which connects the quantities p and U 

16w 2 o 4 cot 2 6» 4w 2 a 2 x ,„ 2 > 

Tl ^ + - 2 = o( 3 ^-^- 5 - 27 

(a 2 — r 2 ) 2 a z — r z 2 

Summarizing (5.21) and (5.27), we obtain the formula for the density 
of the distributed matter that fills the space 

12u; 2 a 4 cot 2 6> 3oj 2 a 2 2 

"75 2^2— + 1 2 = ^-° C • 5 ' 28 

(a 2 — r 2 ) 2 a z — r z 

Multiplying (5.21) by 3, then subtracting (5.27) from the obtained 
product, we obtain the formula for U 

4Lu 2 a 4 cot 2 e co 2 a 2 TT 

-7-2 2^2- + - 2= KU - 5 - 29 

(a z — r z ) z a z — r z 

Comparing the obtained formulae (5.28) and (5.29), we see that p 
and U of the distributed matter which fills the space of the metric (5.3) 
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are connected to each other by the relation 



(5.30) 



Finally, we transform the four tensorial equations of the obtained 
chr.inv. -Einstein equations (5.21-5.26) so that they express the non- 
zero contravariant components of the stress-tensor: U 11 = h lm h ln U mn , 
U 12 = h lm h 2n U mn , U 22 = h 2m h 2n U mn , U 33 = h 3m h 3n U mn . Taking the 
obtained relation U = \pc 2 (5.30) and also the obtained formula for p 
(5.28) into account, we obtain 

8w 2 a 2 cot 2 0 xpc 2 (a 2 - 



hU 1x = 



r 2 ) 



3 a 2 



12 



4w 2 o 2 cot6> 



r(a 2 — r 2 ) 



xU 22 = ^ 



xU 33 = 



2uj 2 a 



2„2 



r 2 sm 2 0 



2uj 2 a 2 8w 2 fl 4 



>tpc 



(a 2 — r 



2^2 



(5.31) 
(5.32) 
(5.33) 
(5.34) 



Now, we have to check whether the obtained chr.inv. -Einstein equa- 
tions (i.e. the given particular type of the distributed matter) satisfying 
the metric (5.3) or not. 

How to do it? The terms consisting Einstein's field equations are of 
two sorts. These are the characteristics of the particular space and the 
characteristics of the matter which fills the space (the latter are the com- 
ponents of the energy-momentum tensor of the matter) . Suppose that 
we have received, in another way, the formulae for the components of 
the energy-momentum tensor of the given matter as expressed through 
the characteristics of the given space. Then, substituting one into the 
other, we will see: if the equations become identities, they satisfy the 
particular space; however if not, then not. 

To find how p, J 1 , and U lk of the obtained chr.inv. -Einstein equa- 
tions are expressed through the geometric characteristics of the space, 
we consider the conservation law equations (1.89-1.90). The equations 
are the extended chr.inv. -notation of the conservation law V a T c " J = 0 
for the energy-momentum tensor of the distributed matter. 

In the space, which is free of deformations, such as the space of the 
metric (5.3) we suggested to neutron stars and pulsars, the conservation 
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law equations (1.89-1.90) take the simplified form 

*^+*V i J i -^F i J i =Q, (5.35) 
at & 

*8 J k ~ 

2A i k J i + *VtU tk - P F k =0, (5.36) 



dt 

where *V{ — *Vi — -^Fi (see Notations). According to the obtained 
chr.inv. -Einstein equations, we have only J 3 ^ 0 of all the three compo- 
nents J 1 of the observable density of momentum in the rotating liquid 
sphere. Also, as was shown in §4.2, only Fi^O therein. Therefore, 
concerning the scalar conservation law equation (5.35), we have 

*V, J 1 -\FiJ* = *V 3 J 3 — \f 3 J 3 = 

^ + J 3 A^-1F3J 3 )-^3J3=0. (5.37) 

As a result, the scalar conservation law equation (5.35) transforms into 
the condition 

f=0, (5.38) 

which means that the observable density of the matter (the liquid sub- 
stance and the fields) that fills the sphere is stationary. 

Of the three vectorial conservation law equations (5.36), the equation 
with the index k — 3 vanishes. The rest two vectorial equations (with 
k = 1, 2) take the form, respectively 

2A 31 (a 2 -r 2 ) j3 dU^_ dU^_ ( dlnVh \ 

a 2 J + dr + 89 + \ 86 ) u + 

+ A\ 2 U 22 + A 33 C/ 33 + (a\, + - 1 F^j U 11 = P F\ (5.39) 

2A, 2 ~ dU 12 dU 2 



J 3 + 




dr dO 

+ A 2 3 U 33 + ( 2A 2 2 + - 1 F! j U 12 = 0. (5.40) 

Consider these two conservation law equations (5.39, 5.40) which 
remain non- vanished for yet. Substitute, into the equations, the charac- 
teristics of the space of the rotating liquid sphere and the characteristics 
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of the matter which fills it. The formulae for U %k (5.31-5.34) and the 
formula for J 3 (5.22) which come from the chr.inv. -Einstein equations. 
The formulae for the logarithmic derivatives (5.7, 5.8). The obtained 
formula for p (5.28). The formulae for the acting gravitational inertial 
force Fx (5.9) and for the non-zero components A13, A23, of the tensor 
of the angular rotation of the space (5.10). When all the formulae have 
been substituted into the remaining conservation law equations (5.39, 
5.40), after some algebra we see that the equations vanish as well. 

So, the common solution of Einstein's field equations and the con- 
servation law equations in the space of the rotating liquid sphere showed 
that the suggested equations are valid in the space. In other words, the 
space metric (5.3) we have suggested to neutron stars or pulsars satisfies 
Einstein's field equations (and vice versa). 

§5.3 Introducing the electromagnetic field 

As is known, every neutron star or pulsar bears the strong magnetic 
field. Therefore, we go to the next stage in this research. We need to 
introduce such an energy-momentum tensor that describes the electro- 
magnetic field and satisfies the relation U — \ p<? (5.30) which follows 
from the obtained chr.inv. -Einstein equations. The energy-momentum 
tensor which satisfies the relation U = 3 pc 2 , satisfies the space metric 
we suggested to neutron stars and pulsars. Once the energy-momentum 
tensor will be obtained, the equations of the electromagnetic field will 
be able to be deduced. Then we conclude how the electromagnetic field 
is distributed inside a neutron star or a pulsar, according to our theory. 
This is our plan for now. 

The energy-momentum tensor of an arbitrary electromagnetic field 
has the following general form 



where F a p is the electromagnetic field tensor (the Maxwell tensor). It 
is the curl of the four-dimensional electromagnetic potential A a 



The physical observable chr.inv. -projections of the vector A a are the 
scalar potential ip and the vector potential q l of the electromagnetic 




(5.41) 



F af} = V a A p - Vp A a = 



dA 0 dA, 



(5.42) 



dx a dx$ ' 



field 



A {) 



q l = A 1 . 



(5.43) 



<P = 
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The theory of the electromagnetic field expressed in terms of the 
chronometrically invariant formalism is well developed in our book [18]. 
(See Chapter 3 therein.) So forth we follow with the theory, and refer 
everyone who is curious in the details to the book. 

The physical observable components of the electromagnetic field ten- 
sor F a fj (5.42) have the form 

Too EiE' 1 + H^H* 1 

Pom = = , 5.44 

.9oo 8ttc^ 

JL = 4^ = e tkm E k H* m , (5.45) 

V500 47TC 

U l e i = c 2 T lk = Pcm c 2 h lk - -L {E l E k + H* i H* k ) . (5.46) 

47T 

Herein, E l and H* 1 are, respectively, the three-dimensional chr.inv. - 
vector and chr.inv. -pseudo- vector of the electric and magnetic strengthcs 
of the field, while e %mn is the unit completely antisymmetric three- 
dimensional chr.inv. -pseudo-tensor. They are expressed as [18] 



-ty — — € H/n i J-Jn — ~ — — — r n 

dx n c dt c 2 

rr*i _ 1 Amn rr rr J^n _ _ A 

- 2 £ Umn ' Umn ~ d X n dx m c 



> . (5.47) 



We see that the observable electric and magnetic strengthes depend 
on not only the electromagnetic field itself (the scalar and vectorial 
potentials, tp and q l ), but also on the acting gravitational inertial force 
Fi and the angular velocity of the space rotation, A^. 

Pulsars are the massive objects which possess the strong electro- 
magnetic field and the rapid rotation. Therefore, the factors of Fi and 
Aik are significant in our study. We will, however, neglect the temporal 
variations and spatial non-uniformities of the electromagnetic poten- 
tials by assuming tp = const and qi = const. With these assumptions, 
the observable electric and magnetic strengthes take the form 

Ei = -^Fi, H n = -^n*\ (5.48) 
c 2 c 

where CI* 1 is the three-dimensional chr.inv. -pseudo- vector of the angular 
velocity that the star rotates 

q*z _J_ imn a q jj_ Aran /r aq\ 
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The components of the chr.inv. -tensor of the angular speed of the space 
rotation, An- , are determined by the metric of the space of the rotating 
star. They, calculated for the metric of a neutron star or a pulsar, are 
presented in formula (5.10). 

Note. The formulae for the electric strength and the magnetic 
strength (5.48) which we finally suggest to neutron stars and pulsars, 
manifest that the electromagnetic field of such a star is due to its grav- 
itation and rotation. Namely, — the electric strength E % of the field 
is manifested only due to the gravitational field of the star (even if the 
electromagnetic field potential (p is presented in the star) . The magnetic 
strength H* % is manifested only if the star rotates. 

Using these formulae for E l and H* 1 (5.48), and all the aforemen- 
tioned assumptions we suggested to neutron stars and pulsars, we trans- 
form the physical observable components of the electromagnetic field 
tensor F a p (5.44-5.45) to the form 

/?om= 2&(^ + ^"^)' (5 ' 50) 
JL = ^£ ikm F k n* m , (5.51) 

v£ = -£*[-hr + ^ ) hlk ~ =3 -r-r- + ™* k , (5.52) 



27rc 2 \ 4c 2 / ^c 2 \ 4c 2 

which corresponds to the vortex-free electromagnetic field of a neutron 
star or a pulsar. From here, we obtain the formula for U cm = 

U em = ^ {^T- + = Pemc 2 . (5.53) 

As is seen from this formula, (5.53), in the framework of the assumed 
conditions of the particular electromagnetic field, we have U = pc 2 . 
However, as we obtained earlier for the space metric of a liquid neu- 
tron star or a liquid pulsar, there should be U = \p<? (5.30). In other 
words, according to the metric, we should have 

C/cm = ^ PcmC 2 , (5.54) 

where 



Therefore, our task now is to find such a physical condition under which 
the electromagnetic field satisfies (5.55), and hence (5.54). 
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Let us find the condition. With the use of the obtained relationship 
U cm = PcmC 2 (5.53), we re-write the (desirable) formula U em = -ifi^fi*- 5 
(5.55) in the form 

or, because k= in the equivalent notation 

c 2 n* j = c - Fj F j . (5.57) 

-, _ 4Gy 2 

Note that the quantity is dimensionless. The scalar electromagnetic 
potential is constant, ip = const, according to our initially assumptions. 
Therefore, and because the magnetic strength is H* % = — ^-fl* 1 (5.48), 
the stationary rotating star is a permanent magnet. 
Denote 

^T-=n, (5-58) 
where n<\, while c and G are the fundamental constants. Therefore, 



c 



.2 



Lp = — = < 1.74x10 



24 



1 /2 rIT1 l/2 



gram ' cm 



(5.59) 



2sfG 

Having the scalar electromagnetic potential ip within this scale of the 
magnitudes, the electromagnetic field satisfies the space metric of a 
neutron star or a pulsar. 

As a result, we re-write the obtained formula (5.57) in the form 

1 

I -An' 1 ' ' " ~ 4 



c 2 n„jn* j = —FiF\ n<-. (5.60) 



With this particular condition, which connects the acting force of gravi- 
tation and the angular velocity of rotation of the space, the electromag- 
netic field satisfies the space metric and Einstein's field equations which 
we suggested to neutron stars or pulsars. 

§5.4 The distribution of the magnetic strength 

To find how the magnetic strength is distributed along the surface of a 
neutron star or a pulsar, we consider Maxwell's equations. The general 
formulation of the two groups of Maxwell's equations is 

V CT F^ = f, V a F*^ =0, (5.61) 
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where F*^ a ~e^ CTQ/3 F Q( g is the pseudo-tensor which is dual to the elec- 
tromagnetic field tensor F Qj g, while j 11 is the four-dimensional current 
vector. 

This formulation of Maxwell's equations means an arbitrary elec- 
tromagnetic field. Transform it with taking into account our assump- 
tions which are particular to neutron stars and pulsars. As previously, 
we neglect the temporal variations and spatial non-uniformities of the 
electromagnetic potentials by assuming ip = const and qi = const. With 
these assumptions, the current vector is zero (j^ = 0) thus Maxwell's 
equations (5.61) take the particular form 



V a F» a =0, 



V a F*^ = 0. 



(5.62) 



Write down the particular Maxwell equations (5.62) according to the 
chronometrically invariant formalism. The chr.inv. -Maxwell equations 
have the form (see Chapter 3 of the book [18] for detail) 



'V,^ - -H lk A. 



0 



*V k H ik - \ F k H ik - - 



*dE l 



+ DE' 



(5.63) 



A ik = 0 



*ViH n --E 
c 



*V k E nk - — F k E* ik - - 



"dH* 



+ DH* 



0 



II, 



(5.64) 



where E*' lk = — e lkn E k is the pseudo-tensor which is dual to the electric 
strength tensor Ei, while D — h %k Di k is the rate of the deformation of 
the space. Because the space of the rotating liquid sphere which is under 
our consideration does not deform, and also, according to our initially 
assumptions, the electric and magnetic strengthes are stationary, the 
chr.inv. -Maxwell equations take the simplified form 



1 



- - H lk A ik = 0 



*V k W k --F k H tk = 0 



(5.65) 



E* lk A; 



ik 



"V k E* lk - —F k E* lk = 0 



> II. 



(5.66) 
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Substitute, into the chr.inv. -Maxwell equation (5.65, 5.66) which 
are already adapted to the space metric of a neutron star or a pul- 
sar, the respective formulae for E % and H lk (5.48) (and those for their 
dual preudo-tensors), and also the respective characteristics of the space 
which we have obtained in §4.1. 

The first (scalar) equation of Group I (5.65) takes the form 

V _ 2 =40. J -n* 3 - (5-67) 

Two of the vectorial equations of Group I vanish, while the third vec- 
torial equation takes the form 

W C ° td 0, (5.68) 



r 3 {a 2 -r 2 )Va 2 -r 2 sin(9 

where u>, according to the space metric of the star (5.3), is the angular 
speed of the rotation of the star along the equatorial axis (f>. Both the 
scalar and vectorial equations of Group II (5.66). Therefore, the dry rest 
which we have from the chr.inv. -Maxwell equations adapted to neutron 
stars and pulsars, are only the equations (5.67) and (5.68). 

Due to the obvious assumption that stars are not point-like objects 
(so, a > 0), and that the radial coordinate is positive (r > 0), we arrive 
at the solely valid solution of the equation (5.68): 

e = ±\- (5.69) 

The solution means that the vectorial equation of Group I of the chr.inv. - 
Maxwell equations is applicable only to the poles of a neutron star or a 
pulsar. 

The vectorial equation of Group I describes the chr.inv. -function 
*VkH lk , which means the observable three-dimensional distribution 
of the magnetic strength H %k of the electromagnetic field of the star. 
Therefore, the solution (5.69) we have obtained means that the mag- 
netic field of a neutron star or a pulsar manifests itself only at the South 
Pole and North Pole of the star. 

Calculate the magnetic strength H n = — 2 ^-W l (5.48) for this case. 
The components of the unit antisymmentric chr.inv. -pseudo-tensor £ tkm 
are explained in detail in Chapter 2 of the book [18]. Thus, after algebra 
we obtain the components of the angular velocity chr.inv. -pseudo- vector 
n n (5.49) of the star 

n* 1 = ^ = u , n»! = A 23 Vh = , (5.70) 

Vh a 2 -r 2 
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\Jh r\a A — r l ) a z — r z 

According to the obtained solution of the chr.inv. -Maxwell equa- 
tions, which is 6* =±77 (5.69), we have cot0 = O in the case. Therefore, 
f2* 2 = f2*2 = 0. This means that the magnetic field of a neutron star or 
a pulsar has the solely non-zero components H* 1 = — tt* 1 , which thus 
is the radial r-components directed from the centre of the star to its 
South Pole and North Pole, then — to the respective polar directions 
from the star into the outer cosmos. 

This result, in common with the solution 6 = ±^, were obtained on 
the basis of our mathematical theory of the liquid neutron stars and 
pulsars. These purely theoretical results completely coincide with the 
well-known observational data about the pulsars. 

§5.5 The frequency and the magnetic strength of a pulsar 

The electromagnetic radiation of a pulsar (rapidly rotating neutron star) 
is the same as the rotational frequency of the star itself. Let us calculate, 
on the basis of our theory of the liquid neutron stars and pulsars, the 
frequency of the electromagnetic radiation of a typical pulsar. 

Calculate Q*j£l*i at the South Pole and North Pole of a rotating 
neutron star (a pulsar), where 0 = ±§. We obtain 

n^n* j = n^n* 1 = — — - . (5.72) 

a 2 - r 2 

Then the condition (5.60), which connects the angular velocity of rota- 
tion of the space and the force of gravity acting in it, takes the form 



iu 2 a 2 



1 — 4n a 2 (a 2 — r 2 ) 



(5.73) 



The magnetic strength of the electromagnetic field of a neutron star 
or a pulsar is expressed as H* 1 = — W % (5.48). It is due to the rotation 
of the star. Therefore, studying of the obtained relation (5.73), we can 
make a conclusion about the electromagnetic radiation of the star. 

The relation (5.73) has a breaking at the surface of the star (r = a). 
We therefore assume r^a. The relation (5.73) thus takes the form 

2 l-4l! (J 2 tt 4 

r 2 = — , 5.74 

n c z 

where r, with taking the previous solution 6 = ±^ (5.69) into account, 
is the radial distance along the polar axis of the rotation of the star. 
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If the electromagnetic radiation is produced near the surface of the 
star, inside the surface layer, we have r ~ a. Thus, after some trivial 
transformations, we obtain the formula for the frequency of the oscilla- 
tion of the magnetic field of the star 



W = W0 VT^' Wo = ^' (5 - 75) 

where wo is the ultimate high rotational frequency of the star, at which 
the star rotates with the light speed. 

Assume a = 10 6 cm which is the typical radius of a neutron star. 
With this radius, we obtain 

uj 0 = 3xl0 4 sec" 1 . (5.76) 

It follows, from (5.75), that n is expressed as 

,,,2 

(5.77) 



■4w 2 ■ 



The observed frequencies of the radio-pulsars are in the range between 
w m in = 0.53 and w max = 448.57 sec^ 1 . This means that uj 2 <^ujq. We 
therefore can neglect u> in the denominator of (5.77). We obtain 



2 2 



n=^ = —. (5.78) 



Therefore, for the real radio-pulsars, the number n lies in the range 

3.1 xl(T 10 < n < 2.2 xl(T 4 . (5.79) 

Also, according to formula (5.58) deduced in the framework of our 
theory, the scalar potential of the electromagnetic field of a pulsar is 

2 !~n 



<p = <?\Ig- (5 - 80) 

Consequently, for the real radio-pulsars, we have 

6.1 xlO 19 < Lp < 5.2 xlO 22 , (5.81) 

in [gram 1 / 2 cm 1 / 2 scc -1 ]. This interval of the magnitudes of (p satisfies 
the upper theoretical limitation on the potential, which, according to 
our theory, is tp < 1.74 xlO 24 gram 1 / 2 cm 1 / 2 sec -1 (5.59). 

Finally, we now will calculate, on the basis of our theory, the ex- 
pected range of the magnitudes for the magnetic field strength H* 1 of 
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the pulsars. According to our theory of the liquid neutron stars and 
pulsars, H* l = —^-fl n (5.48). With the calculated range of the magni- 
tudes of the scalar electromagnetic potential ip, and with the esteemed 
range of the rotational frequencies u of the pulsars, we obtain 

2.1 xlO 9 < H* 1 < 1.5xl0 15 , (5.82) 

in [gram^cm^^sec -1 ]. This range of the numerical values very cor- 
responds to the magnitudes of the magnetic field of the radio-pulsars, 
which are known due to the radio-astronomical observations. 

§5.6 Solving Maxwell's equations in the stationary vortex- 
free magnetic field of a neutron star 

Previously, in §5.4— §5.5, we studied Maxwell's equations in the elec- 
tromagnetic field of a neutron star or a pulsar under the assumption 
that the four-dimensional current vector j a was zero within the field 
(j a = 0). See (5.62) and so forth. In other words, we assumed that the 
electromagnetic field was free of currents. 

Further, this assumption generates the following problem. Look at 
the formula for the observable momentum of the electromagnetic field 
J* m (5.51), which is the Poynting vector of the field. It has the form 

Jim = -^~ c £ lkm E k H stm = e lkm F k Q* m . (5.83) 

By the current-free assumption j a — 0 made in §5.4, we obtained that 
only the component H* 1 ~ — ^O* 1 of the magnetic field strength H* 1 
was non-zero at the South Pole and North Pole of the star. In this case 
the circular momentum of the field, Jg m , which should generate the mag- 
netic component H* 1 , would have been zero: J^ m = ^ e zl2 EiH*2 = 0. 
This causes some trouble, because a model that satisfies astronomical 
observations of the pulsars should obviously show H* 1 = — ^ Q* 1 =/= 0 
andJi = ^£ 312 i?iff,2^0. 

Recall that we previously arrived at the difficulty that iJ* 1 ^0 but 
Jg m = 0 as a result of our assumption according to which the electro- 
magnetic field was free of currents (j a — 0). Therefore, we now will solve 
Maxwell's equations in common with the condition j a ^=0. 

First, we will resolve this problem in the vortex-free electromagnetic 
field. In §5.7, this problem will be solved in the vortical field. 

The space (space-time) metric of the rotating neutron star has the 
form (5.3). This metric means that the liquid sphere is not a collapsar 
due to its rotation (see the necessary explanation in the beginning of 
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this Chapter) . Physical and geometric characteristics of the space were 
calculated and presented in §5.1. In addition to these, we only should 
add that the pseudo-vector of the angular velocity of space rotation, 



il* 1 — \ e lmn A mn , has the following components 



f2* — to , 



2 2wa 2 cot# 2uja 2 r cot 9 

r(a z — r z ) a z — r z 



(5.84) 



Respectively, the square of the pseudo-vector of the angular velocity is 

W V / 4aW0\ 

^ = ^2~T2 ^ + -^2-^2- ) ■ ( 5 - 85 ) 

Assume now that the scalar electromagnetic potential of the field re- 
mains unchanged, <p = const, while the vectorial electromagnetic poten- 
tial gi is vortex-free. Then the components of the electric and magnetic 
field strengths (5.47) take the form 



c z cr 

Tj*i \ imn tt tt \ 

2 c 



(5.86) 



Herein, using the definition of Q**, which is fi** = ^e lmn A mn , we rewrite 
the formula for H* % in the form 

H* i — — Q**, H*i = — . (5.87) 

c c 

Using the formula for F\ (5.9), then calculating fJ* 1 and ft* 2 from (5.84), 
we obtain the substantial components of E % and H* % . They are 

ff.i = ~^4, F- = -^, (5.89) 
c(a z — H) c 

4yua 2 rcot6> 2 4</^a 2 cot6> 

c(a z — r z ) cr\a l — r z ) 

Let us find how the magnetic strength is distributed along the surface 
of the sphere. Consider Maxwell's equations in their full form (5.61) 

V CT F^ = f, V a F*^ = 0, (5.91) 
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which implies the presence of the field current (j a ^0)- According to 
the chronometrically invariant formalism, their physically observable 
(chronometrically invariant) projections — the chr.inv. -Maxwell equa- 
tions — have the form (see Chapter 3 of the book [18] for detail) 



*VjEi --H ik A ik =4Trp 



ik 



ik 



1 (*dE l 
c 



dt 



DE< )=^f 



> I, (5.92) 



1 



E* %k A lk = 0 



1 1 / *8 H* 1 

*V k E* tk - — F k E nk - - ( — — + DH" 



> II. (5.93) 



Herein, E* lk = — e tkn E k is the pseudo-tensor which is dual to the electric 
strength tensor E i7 while D = h lk D ik is the rate of space deformation. 
Because the space of the rotating liquid sphere which we consider does 
not deform, and also, according to our initial assumptions, the electric 
and magnetic strengths are stationary, the chr.inv. -Maxwell equations 
take the simplified form 



1 



VjE 3 - -H ik A ik =Anp 



■ FkH 



ik 



4-7T 



> I, 



(5.94) 



*V 4J ff" - -E* lk A lk = 0 



II. 



*V k E 



*ik 



■F k E 



(5.95) 



The first equation of Group I (5.94) takes the form 



2„2 



4<pw a 
c 2 (a 2 



1 + 



4a 2 cot 2 1 



y(3a 2 - 2r 2 ) 
a 2 (a 2 — r 2 ) 



= 47rp. (5.96) 



It follows from the second equation of Group I that j 1 = j 2 ~0 in the 
framework of our model, while the equation for j 3 takes the form 



(paja' J 



cote 



r 2 (a 2 -r 2 )Va 2 - r 2 sin 6 



= -T 



(5.97) 
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and the absolute value of the chr.inv. -current vector f is 

(puia 3 cot 8 



J = VU? = — T ^-— T == . (5.98) 
irr (or — r z ) Va — r z 

Equations of Group II (5.95) satisfy themselves as identities. So, these 
formulae for p, j 3 , and j (5.96-5.98) are the exact solutions to the 
Maxwell equations we have just considered. 

So, we have obtained the exact solutions to Maxwell's equations in 
the internal electromagnetic field of a neutron star or a pulsar, where 
the field originates due to its sources: the distributed charges p and the 
currents f . 

The conservation law of the electric charges and the currents sets 
up a connexion between the sources of the electromagnetic field. The 
generally covariant form of the conservation law, which is also known 
as the continuity equation, has the form 

V CT r = 0. (5.99) 

It means that the distributed charges p and the currents j l , which are 
the respective physically observable (chronometrically invariant) pro- 
jections of the four-dimensional current vector j a , are conserved within 
the four-dimensional volume of the field. Also, Maxwell's equations are 
connected by the Lorenz condition 

V ff A ff = 0, (5.100) 

which is imposed on the four-dimensional electromagnetic potential A a . 

In a general case, the conservation law V CT j a = 0 and the Lorentz 
condition V CT A a = 0 written in terms of the chronometrically invariant 
formalism, have the form (see Chapter 3 of the book [18]), respectively, 

*-^ + pD + *V i j i -i>F i j i = 0, (5.101) 

-*-^f + -D + *V t q l -^F l q^0 1 (5.102) 

where *V, = *V, — \Fi (see Notations for definition of *V,). 

It is easy to show that, under the particular conditions of the prob- 
lem we are considering, the chr.inv. -continuity equations (5.101) and 
the chr.inv. -Lorentz conditions (5.102) are satisfied as identities. 

Now, on the basis of the obtained exact solutions (5.96-5.98) of the 
Maxwell equations, we look for the Poynting vector J* m that is the 
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observable momentum of the electromagnetic field. We need to know 
how the Poynting vector is distributed along the surface of the sphere, 
which is the surface of a neutron star or a pulsar. 

The Poynting vector J* m is the second of the three physically ob- 
servable projections of the electromagnetic field tensor F a p (5.42), which 
are p em (5.44), J» m (5.45), and (5.46). We look for the Poynting 
vector Jg m (5.45) in the framework of the particular conditions accord- 
ing to which the scalar potential ip of the electromagnetic field remains 
constant, while the vectorial potential % of the field is vortex-free, i.e. 

*dqi *dq k 
ip = const , — - — — = 0 . 

OX K ox % 

Substituting, into (5.44-5.46), the substantial (non-zero) compo- 
nents of the electric strength E % and the magnetic strength H* 1 , which 
are (5.88-5.90), we obtain 



Pem — 



ip 2 ( FjF j 
Ac 2 
uj 2 a 2 



2ttc 4 

2 



27TC 4 



4u 2 q 4 cot 2 6> 



c 2 r 2 



a 2 — r 2 (a 2 — r 2 ) 2 4a 2 (a 2 — r 2 ) 



t "k^ &*m . r— 

2nc 4 Vh 



2ttc 4 



(p 2 uia 



cote 



ire 2 (a 2 — r 2 ) 3 / 2 sin 9 ' 

ip war cot 0 



ttc 2 (a 2 — r 2 ) 3 / 2 ' 



, (5.103) 

(5.104) 
(5.105) 



Looking at these equations, we can conclude something about the 
neutron stars and pulsars whose electromagnetic field is vortex-free (the 
case of the §5.6). So, we see that the electromagnetic field density p cm 
is due to the gravitational inertial force, which is the non-Newtonian 
force of repulsion F{ acting within the sphere, and due to the sphere's 
rotation. The electromagnetic field density p em can be non-zero by the 
separate condition Fi ^ 0 or A ik ^ 0, and the common condition Fi^O 
and Aik ^ 0. The density of the field momentum J l em is non-zero only 
by the common condition Fij^O and ^ 0. 

As follows from (5.103), the density p em of the vortex-free electro- 
magnetic field of a rotating neutron star (a pulsar) is zero at the equator 
of the star (9 = 0). Then the field density p em increases with the geo- 
graphic latitude 9 to the South Pole and North Pole, where 6 = ^ so 
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the density of the vortex-free electromagnetic field takes the ultimately 
high magnitude p em = (p em ) 

max- 

Contrarily, the density of the electromagnetic field momentum Jg m 
(5.105) is ultimately high at the equator, where 8 = 0. Then the magni- 
tude of the field momentum Jg m falls down with the geographic latitude 
9 to the South Pole and North Pole, where it is J^ m = 0. 

In addition to these, we can also conclude something about the 
charge density p and the currents j l within the vortex-free electromag- 
netic field of a neutron star or a pulsar. 

We can initially re- write the respective formulae for the charge den- 
sity p (5.96) and the current j 3 (5.97), obtained from Group I of the 
chr.inv. -Maxwell equations, as follows: 

■ ?3 = ~nr 2 (a 2 -r 2 ) 3 / 2 sm9 ' (5A ° 7) 



where 



Vj Fi = ^n 2 2? >°' ( 5 - 108 ) 



3 = 



rr — - 3 - 3 tpua 3 cot 6 

VWmjfm - 7rr{a 2_ r 2 )3 /2 ' ( 5A09 ) 

As a result, we see that the electromagnetic charge density within a 
neutron star or a pulsar is positive p > 0 (that should be according to 
the physical sense of a physical field) if 

<>. ,<>• ' /• '. (5.110) 

Now re- write this inequality with the formula for p (5.96). We obtain the 
following condition which is necessary according to the physical sense: 

4oj 2 a 2 ( 4cot 2 <9\ 3a 2 -2r 2 

> 5 . (5.111) 



Compare the formulae for the electromagnetic field current j 3 (5.107) 
and its power j (5.109) to the obtained formulae for the density of 
the electromagnetic field momentum J 3 m (5.104) and the power of the 
momentum J om (5.105). We have 

c2j L = -^-f, c 2 J cm = ^-j. (5.112) 
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Taking (5.57) into account, we express the scalar electromagnetic 
field potential ip (which is ip = const, according to our initial assump- 
tions) through the dimensionless constant n = (5.58) (the constant 
is n < j, see in the end of §5.3, for detail). So, we have 

9 fn 9 nc 4 1 

y = c 2 ^j-, ^=-q^ n< l- ( 5 - 113 ) 

With these, we obtain 



3 _ EiH^ _ 4nc 3 ua cot 9 (kii^ 
cm ~ sfh ~ G ( a 2-r 2 )3/ 2 sin0 ' [ ' 



/ n,.fi*i_ y,-^' , (5.116) 



3 ~~ C V G r 2 (a 2 - r 2 ) 3 / 2 sinfl ' ( } 

("») 

We see that the greater the scalar electromagnetic potential ip (5.113) 
of a neutron star or a pulsar, the stronger the three-dimensional circu- 
lar current j 3 and the three-dimensional circular momentum J 3 m of the 
electromagnetic field. Moreover, the current and momentum flow of 
the electromagnetic field exist in the star only if the star rotates on 
the equatorial plane (x 1 ^ 3 ), i.e. only if f2»2 7^0. If the neutron star 
does not rotate (Cl^jfl*^ =0), the electric charge density of its internal 
electromagnetic field would be negative (p < 0). 

So, we have arrived at a non-satisfactory result. Both the circu- 
lar electromagnetic field current j 3 (the current j l that goes along the 
longitude coordinate <ft) and the electromagnetic field momentum J 3 m 
(that is the Poynting vector of the field) are zero at the South Pole and 
North Pole of the star, where the geographical latitude is 9 = ^, and 
reach the ultimate power at the equator (9 — 0). 

Herein, we have assumed that the electromagnetic field of the ro- 
tating neutron star is vortex-free. The final feat to match with the 
observational data will be done with the vortical electromagnetic field 
of a rotating neutron star (a pulsar). We will do it next, in §5.7. 
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§5.7 Solving Maxwell's equations in the stationary vortical 
magnetic field of a neutron star 

In analogy to §3.7, consider a rotating neutron star (a pulsar) whose 
electromagnetic field is vortical. The curl of the three-dimensional 
vectorial chr.inv. -potential qi of the field is non-zero 



qik = 



*dqi *dq k 



?0. 



dx k dx' 

The four-dimensional electromagnetic field potential 



dx ° 

A = V ~T~ ' 

as 

has two chr.inv. -projections 
A 0 



9a/3 



dx a dx 13 
ds ds 



= 1. 



= <P, 



where 



'9oo 



2 i i k 

v = h ik v v , 



Ai i <P i ; dx 1 
A l = q l = -v\ v 1 = — 
c dr 



v 2 < c 2 , = 



= ¥>■ 



(5.119) 



(5.120) 



(5.121) 



(5.122) 



According to our initial assumptions, ip = const and q 1 = q 2 = 0 in the 
field. Thus v 3 = ^ =u, and the non-zero components of the vectorial 
electromagnetic potential q l and the field curl q ik have the form 



q 3 = —, 

c 

IfUJ 2 . 2 

qs = r sin i 

c 



931 



923 



dqz = 

dr 
dq 3 
86 



2<pu> 



r sin 9 , 



2tpU 2 . n 

r sin ^ cos ( 



(5.123) 
(5.124) 

(5.125) 

(5.126) 



Using the definition of the field strengths (5.47), we calculate the 
non-zero components of the magnetic strength of the vortical field 



H23 — — 



2(pwr 2 sm9 



a 



H31 — 



2ipu> 



sin 2 9 - 



2 a 3 cos 9 



{a 2 — r 2 ) \J a 2 — r 2 



(5.127) 
(5.128) 
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Using the relation 



H » = _ ^mn Amn = _ £ *mn qmn _ ^_ Q „ ^ (g 12g) 

we re-write the components H23 (5.127) and -H31 (5.128) in the form 
2tpu> 



H* 



1 + 



\J a 2 — r 2 



cos 8 



(5.130) 



H 



*2 



2tpuj ( yja 2 — r 2 sin 9 2 a 2 cot ( 



cr 



(5.131) 



while their covariant (lower-index) versions can be calculated as follows: 
H,i = /in * x and #* 2 = /i 22 H* 2 . 

Let us then find the solution to the chr.inv. -Maxwell equations. 
In the case where the electromagnetic field is stationary, the chr.inv.- 
Maxwell equations have the form (5.94-5.95) 



1 



*V k H* k --F k H tk = — f 



1 



4-7T 

c 



> I, 



E* ik A lk = 0 



*V k E* lk - -F t £ nft = 0 



> 11. 



(5.132) 



(5.133) 



After substituting the electric and magnetic strengths of the vortical 
electromagnetic field, we see that equations of Group II (5.133) are 
satisfied as identities. Equations of Group II take the form 



4(puj 2 



1 



4 a 2 cot 2 6» 



a cos 6 



v^ 2 " 
^(3a 2 -2r 2 ) 
a 2 (a 2 — r 2 ) 



4irp , 
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cote 



3 (fiUl 
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(5.134) 
(5.135) 



where p and j 3 are the change density and the current of the vortical 
electromagnetic field. The physical sense of these equations readily looks 
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more understandable when rc 



ritten in the form 



cote 



ipuJ 1 acos9 
ttc 2 yj a 2 — r 2 

3<puj 



7rr 2 (a 2 — r 2 ) 3 / 2 sin0 2-na 2 



(5.136) 



(5.137) 



where uj = fl* 1 . Express now the charge density p and the current j 3 
of the vortical electromagnetic field through the same characteristics p 
(5.106) and j 3 (5.107) we have calculated in the vortex-free field 



P = P- 



(fui 2 a cos 9 



TTC 



p = .3 _ 

J J a 9 



(5.138) 
(5.139) 



As is seen from the equations (5.138) and (5.139), given a rotating 
neutron star (a pulsar) whose electromagnetic field is vortical, the charge 
density and currents of the field differ from those of the vortex-free 
electromagnetic field by those terms depending on the star's rotation. 

Respectively, the field density p cm (5.44) and the circular momentum 
flow J 3 m (5.45) of the vortical electromagnetic field have the form 



[ — FF J 

2irc 4 Uc 2 0 



n^n* 3 ] + 



2ttc 4 

•2 



auj cos o 
\/a 2 — r 2 



om 2^c 2 a 2 
Or, in the other notation, 



2cva 3 



cot0 
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( a 2_ r 2)3/2 sin6 ) 



r 2 sin 2 ( 



(5.140) 
(5.141) 

(5.142) 

., ,. (5.143) 

To understand the meaning of these resulting formulae, recall that, 
as follows from the formulae for A$\ (5.10) 

2u;a 3 r cosd 2ujacot9 
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uj 2 a cos 9 
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( a 2_ r 2)3/2 



A 31 = 



•ija 2 — r 2 sin( 



(5.144) 
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this component and hence the pseudo- vector fi* 2 = ^e 231 A 3 i depend on 
the geographical latitude 6, while fi* 1 = ^e 123 A 2 3 ~cj does not. 

The obtained formulae for the current vector j 3 (5.139) and the 
Poynting vector J 3 m (5.143) of the vortical electromagnetic field contain 
that very term which does not depend on the geographical latitude of 
the star that possesses the field. This means that, contrary to the 
vortex-free electromagnetic field, the current vector j 3 and the flow of 
momentum J 3 m of the vortical electromagnetic field are non-zero at the 
South Pole and North Pole of the star. 

The obtained result, J 3 m ^0 at the South Pole and North Pole, 
means that a rotating neutron star whose electromagnetic field is vor- 
tical can emit electromagnetic radiation along its polar axis, while that 
possessing a vortex-free electromagnetic field — cannot. 

Also, we have to make one more important conclusion in the frame- 
work of our mathematical theory of pulsars. Look at the definition of 
the magnetic field strength B* % =-\ e lmn H mn (5.47), which is 



H *i = 1 £ imn ( *jhm _ *jhn _ ^ A 

2 \ dx n dx m 



ran 

C 



1 „imn ( * d 1m *dq n \ 2ip_ ^ 



2 V dx n dx m 

= \e mm q m n-—W ;i . (5.145) 
2 c 

As is seen from this definition, the pseudo- vector of the magnetic field 
strength H* 1 is the sum of the pseudo- vector of the electromagnetic field 
curl q* 1 = i e lmn q mn and the pseudo- vector fl* 1 of the angular velocity of 
the star's rotation. The magnetic field strength H* 1 coincides with the 
pseudo- vector of the star's rotation, fi", only if the electromagnetic field 
curl q mn is zero. If q mn ^= 0 which is true in a vortical electromagnetic 
field, H* % deviates fi". The stronger the curl q mn of the electromagnetic 
field is, the more the magnetic axis deviates from the axis of the star's 
rotation. 

Astronomers inform us that the electromagnetic field of observed 
pulsars is very strong. They imply that electromagnetic radiation can 
leave such a star only at the polar regions, where the latitudinal and 
longitudinal electromagnetic field components are not so strong as at 
the equatorial latitudes. Also, according to the oscillating behavior of 
the signal registered from pulsars, astronomers conclude that the axis 
of emission and the axis of rotation of the pulsar does not coincide. All 
these facts from observational astronomy match with our theoretical 
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results about pulsars. 

As a result, our mathematical theory of pulsars leads us to the fol- 
lowing conclusions that match observational data: 

• A rotating neutron star can be a pulsar only if its electromagnetic 
field is vortical. Moreover, the curl of the electromagnetic field 
means that the magnetic axis does not coincide with the axis of 
the star's rotation. Otherwise, in the neutron star whose electro- 
magnetic field is vortex-free, electromagnetic radiation docs not 
come from the South Pole and North Pole of the star. 

All these theoretical results have been obtained in the framework 
of our assumption that the scalar and vectorial electromagnetic field 
potentials of the star do not depend on time. Of course, some temporal 
variations of the potentials should pose an effect on the Poynting vector 
(the flow of momentum) of the field, and thus on the electromagnetic 
radiation emitted by the pulsar. But now we neglect these effects. 

§5.8 Geometrization of the vortical electromagnetic field of 
a neutron star 

Geometrization of the electromagnetic field is one of the primary tasks 
in the General Theory of Relativity. As was shown already by Ein- 
stein, this problem in a general case is very non-trivial from the side of 
mathematics. So, it is still not resolved in general. Nevertheless, ge- 
ometrization of the electromagnetic field is possible in particular cases, 
under some particular conditions that simplify the mathematics. 

We now show that in the particular case of a pulsar the electromag- 
netic field is geometrized. In the language of mathematics this means 
that once we have Einstein's field equations and Maxwell's equations, 
the characteristics of the electromagnetic field can be expressed through 
only the geometric characteristics of the space. 

Consider Einstein's field equations (5.18-5.20) and Maxwell's equa- 
tions (5.132-5.133) we have obtained in the internal field of a rotating 
neutron star. Note that in the case of the de Sitter-like metric we have 
derived for neutron stars, the A-term describes physical vacuum in the 
state of inflation X = np (sec Chapter 1 for detail). Also, as we showed 
in §5.2, this form of Einstein's equations satisfies the equations of con- 
servation in the space. 

We will consider the vortical electromagnetic field. This is because 
we have shown that only the vortical field gives the result that matches 
with astronomical observations of pulsars, i.e. the fact that a pulsar 
emits electromagnetic radiation from only its polar regions. 
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We have the scalar electromagnetic field potential ip (which remains 
unchanged for each particular star, according to our initial assump- 
tion) expressed through the fundamental constants as if — c 2 (5.58), 



G 

where n< \ (see in the end of §5.3). With it, we obtain the electric and 
magnetic strengths of the vortical field (see §5.7) in the form 



uj\Ja 2 — r 2 



-2c A /- I ft* 1 + — — cos6M , (5.148) 
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a 



2 



= h 11 H* 1 = -s . IT 1 , (5.150) 



H* 2 = h 22 H* 2 = r 2 H* 2 . (5.151) 



Herein, according to the internal space metric of a rotating neutron star 



or a pulsar, we have ft* 1 (5.70), ft* 2 (5.71), ft^ft*J (5.85): 



ft* 1 = ^=^, ft,x = A 2S Vh = , (5.152) 

yfh a 2 -r 2 

ft* 2 = ^^^, n^A^^, (5.153) 
Vh r(a z — r 2 ) a z —r 2 

As is seen from these formulae, both the electric and magnetic field 
strengths are expressed here through only the geometric characteristics 
of the internal space of the pulsar. 
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Also, in the same way, the charge density p (5.136) and the current 
vector j 3 (5.137) of the vortical electromagnetic field from Maxwell's 
equations are expressed as: 
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P = ~ 

7T 



= p- 



n oj a cos tt 
G y/a 2 — r 2 



n w a cost 
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(5.155) 
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Respectively, the density p em (5.140) and the flow of momentum 
Je m (5.141) of the vortical electromagnetic field — the characteristics 
that come from Einstein's equations — are expressed as: 
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We see that all the characteristics of the vortical magnetic field are 
uniquely expressed through only the geometric characteristics of the 
space inside the pulsar. Therefore, the vortical electromagnetic field of 
a rotating neutron star (a pulsar) is hereby geometrized. 

This fact also means that the system of Einstein's equations and 
Maxwell's equations in the internal space of a pulsar is a self-consistent 
system of equations. This self-consistent system of Einstein-Maxwell 
equations completely describes both gravitational and electromagnetic 
phenomena inside the pulsar. 

However, if the electromagnetic field of a rotating neutron star is 
vortex-free, Einstein's equations and Maxwell's equations do not com- 
prise a self-consistent system. The electromagnetic field is not geometr- 
ized inside the star as such. As was shown in §5.7, such a neutron star 



144 



Chapter 5 Neutron Stars and Pulsars 



cannot emit electromagnetic radiation from its polar regions. Therefore, 
it cannot be a pulsar. 

§5.9 Boundaries of the physical space a pulsar 

Consider an observer whose reference frame is connected to the internal 
space of a star. Where, from his point of view, does the observable 
physical space of the star end? At which distance from the star? 

These questions are answered in the framework of the theory of 
physically observable (chronometrically invariant) quantities in General 
Relativity. In terms of physical observables, the real physical space that 
is allowed to be registered by an observer "ends" at that distance from 
him where the physical observable time stops: dr = 0. The physical 
observable time r is calculated according to the metric of the observer's 
space. Therefore, the real physical boundaries of the observer's ob- 
servable space are determined by the stopped time condition dr = 0 
according to the space metric. 

Let us calculate the boundary of the observable space of a pulsar. 
This is the distance from the center of the pulsar at which, according 
to the space metric of the pulsar (an observer whose reference frame is 
connected with the pulsar), the observable time stops (dr = 0). 

As follows from the chronometrically invariant formalism, the ob- 
servable time interval is formulated as (1.30) 

dr = Jg^dt^ ^ 0t dx 1 = J~goo dt — \-Vidx l . (5.159) 

It consists of two terms. The first term is due to the gravitational field 
potential w = c 2 (l — ^goo) (1-44). The second term is due to the fact 
that the space rotates, and is dependent on the linear velocity of the 
rotation v i =-c- a M= (1.45). 

Therefore, the condition dr = 0 by which the observable time stops 
in the space of a gravitating and rotating body is expressed as 

y/goo dt = \ Vidx\ (5.160) 



c 2 

The space (space-time) metric of a rotating neutron star has the 
form (5.3). See §5.1 for details. In the metric, 

5oo = J (l- J), (5.161) 

2uiar 2 cos# ,„ „ „„. 

v 1 =v 2 =0, v 3 = = ■ (5.162) 

V» — r z 
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In this case, the stopped time condition (5.160) takes the form 



2 



1 / dx 3 \ 
9oo = ^[v 3 —j, (5.163) 

where ^ = ff=w. Substituting the formulae for g Q0 (5.161) and v 3 
(5.162) into the stopped time condition (5.163), we obtain the formula 
for the distance r at which the observable time stops 



r = = . (5.164) 

, 2 C 

73" 



This formula, (5.164), shows the boundary at which the physical 
space of the rotating neutron star ends. As is seen from this formula, the 
boundary r is the same as the neutron star's radius a at the South Pole 
and North Pole: the geographical latitude is 6 = ^ therein, so cos 8 = 0 
and thus r = a according to (5.164). Then the boundary r decreases 
near the equator where it takes the ultimately low numerical value 

(5.165) 



sjl + 



4a 2 to 2 



which depends only on the neutron star's radius a and the angular 
velocity of its rotation uj. 

The more rapid the neutron star rotates, the more oblate the phys- 
ical space of the star at its equator becomes. According to our formula 
(5.165), the oblateness manifests itself only at relativistic speeds of ro- 
tation, i.e. in pulsars. 

Consider PSR J1748-2446ad, that is the fastest-known pulsar dis- 
covered in 2004 [30]. It rotates at a period of 0.00139595482(6) sec 
which means the angular velocity w = ^f- = 4, 501 sec -1 . Its radius a is 
estimated to be less than 16 km. Proceeding from these observational 
data, we can calculate the oblateness of the physical space of the pulsar 
at its equator: 

' 0.90. (5.166) 



yjl + 



4a 2 uj 2 



§5.10 Conclusion 

So, the complete mathematical theory of the liquid neutron stars and 
pulsars is presented here, in this Chapter. We now repeat the most im- 
portant conclusions we made on the basis of the theory. They are: 
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1. As follows from our mathematical theory, the electromagnetic field 
of a rotating neutron star (a pulsar) is due to its rotation and grav- 
itation. The faster the star rotates, the stronger is the magnetic 
strength H* 1 of the field; 

2. The magnetic field strength H* 1 of the pulsar is directed strict 
along the polar axis of its rotation. Electromagnetic radiation is 
emitted only from the poles of the star, then comes into the outer 
cosmos strictly along the axis of the rotation of the star; 

3. The electric strength Ei depends on the spatial distribution of 
the scalar potential and on the temporal variation of the vecto- 
rial potential of the electromagnetic field. The magnetic strength 
H* 1 depends on the curl of the vector potential the field, and on 
the angular velocity of the star. Thus the temporal and spatial 
variations of the electromagnetic field potentials should affect the 
outcoming electromagnetic impulses emitted by the pulsar; 

4. The Poynting vector (the electromagnetic field momentum) is non- 
zero at the South Pole and North Pole of a rotating neutron star 
only if its electromagnetic field is vortical. Therefore, a rotating 
neutron star is a pulsar, thus emitting electromagnetic radiation 
from the polar regions, only if possesses a vortical electromagnetic 
field. Also, the presence of the field curls means that the magnetic 
axis does not coincide with the axis of the star's rotation. A rotat- 
ing neutron star whose electromagnetic field is vortex-free cannot 
emit electromagnetic radiation along its polar axis, so it cannot 
be a pulsar. 

All the conclusions are valid only for a rotating star whose physical 
radius is close to its Hilbcrt radius. These are rotating neutron stars 
and pulsars, not the regular stars such as the Sun etc. 
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§6.1 Non-rotating liquid collapsars. The main characteristics 

Now, we are going to study the collapse condition of a non-rotating 
sphere of perfect liquid (that is, a collapsed rotation- free star), in terms 
of our new model of liquid stars. At first sight, this problem statement 
sounds like non-sense: perfect liquid is incompressible, hence such a 
liquid body cannot be compressed. Yes, it would be non-sense, if one 
would consider collapse as the process of compression of a liquid cos- 
mic body. We do not do it that way: in contrast, we do not discuss 
cosmogony. We merely consider a collapsar as an already existing ob- 
ject. This amounts to the occurrence of the physical conditions, not the 
evolutionary compression of a liquid cosmic body. 

Hence, a cosmic body is a collapsar if the parameters of its field on 
its physical surface correspond to the condition of gravitational collapse. 
Namely, — the field of gravity is so strong on the surface of the body that 
light signals cannot leave the body for the outer cosmos. In terms of the 
General Theory of Relativity, this means that the physically observable 
time stops on the surface. 

According to the theory of physically observable quantities (chrono- 
metrically invariant formalism), the physically observable time interval 
dr (1.30) is formulated in terms of the gravitational potential w and the 
linear velocity of space rotation Vi as follows: 

dr = y/goo dt H dx l = 

= { l -^) dt -^v i dx i . (6.1) 
Thus the general condition of gravitational collapse has the form 

dr = J^dt+-?^dx l = 0. (6.2) 

In a rotation- free space (wherein ^ = 0), the general condition of grav- 
itational collapse is as simple as 



dr = yfg^ dt = 0 



(6.3) 
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or merely 

2 

5oo 



(l-J) 2 = 0. (6-4) 



Thus, a cosmic object of rotation-free space is a collapsar, if the three- 
dimensional gravitational potential w on its surface takes the value 



w = c 2 . (6.5) 



Consider the collapse condition in the case of a non-rotating star 
consisting of perfect liquid. As is seen from the space metric of such a 
liquid star, which is (2.76) 



„2_ 1 ( Q /, "Poa 2 I xpor 2 \ o 



^=4 \*r-^—r-^T- cdt ~ 



-i _ itpgT- 
1 3 

or, in terms of the Hilbert radius r g (2.78), 



' / '"' -r 2 (d0 2 +sin 2 0# 2 ), (6.6) 



4 



<U " r 2 (d6 2 +sm 2 6dct> 2 ), (6.7) 



1 3 



the collapse condition (goo = 0) in this case has the form 



Kp 0 a 2 xpor 2 



3 A /1_Z^_ A /1_Z^ = 0 , (6.8) 
or, similarly, 

Thus, we obtain the radial coordinate r, by which a non-rotating liquid 
star whose radius is a meets the state of gravitational collapse: 



la* 



r c = J9a 2 . (6.10) 

V r 9 

Since we keep in mind real cosmic objects, the numerical value of r c 
should be real (as well as a and r g ). This requirement is satisfied by 
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the following obvious condition: 



a 1.125 r g . 



(6.11) 



If this condition holds not (a ^ 1.125 r g ), the non-rotating liquid body 
(star) cannot be in the state of gravitational collapse. 

The general collapse condition (6.11) includes the particular con- 
dition a = r g . Given this particular case of a collapsed non-rotating 
liquid star, we see that the physical radius a of the star's surface, the 
Hilbcrt radius r g of the star, and the radius of the outer space breaking 
fbr = y/a 3 /r g of the star's field coincide: 



The said collapse condition, a — r g , is only a particular case of the gen- 
eral collapse condition (6.11). The general collapse condition (6.11) 
includes three particular cases, concerning the location of the physical 
surface of the collapsed liquid star: 

1) The collapsed liquid star is larger than the Hilbert radius calcu- 
lated for the star (a>r g ) but less than 1.125 r g ; 

2) The surface of the collapsed liquid star coincides with its Hilbert 
radius (a — r g ); 

3) The collapsed liquid star is completely located within its Hilbert 
radius (a<r g ). 

It is obvious that r c is imaginary for r g -C a, so the state of gravitational 
collapse is impossible for such a star. For example, considering the Sun 
(a = 7xl0 7 cm, M = 2 xlO 33 gram, r 9 = 3xl0 5 cm), we see that r c (6.10) 
takes an imaginary numerical value. The same is as well true for other 
regular stars, ranging from super-giants to white dwarfs. Hence, regular 
stars cannot collapse. 

In fact, the particular collapse condition r c — r\, r ^r g = a (6.12) for- 
mulates the collapse radius r c as follows* 



For example, if a collapsed liquid sphere should consist of regular water 
(po= 1-0 gram/cm 3 ), its radius would be r c = 4.0xl0 13 cm = 3.1 AU, i.e. 
be located within the asteroid strip (the asteroids are located, approxi- 
mately, from 2.1 AU to 4.3 AU from the Sun). 



r br = r g = a. 



(6.12) 




(6.13) 



* x — 8ttj2 _ ig q xio 28 cm/gram is Einstein's constant of gravitation. 
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One more example: to be a collapsar whose size is that of neutron 
stars and pulsars (their radius is a= (8-16) *10 5 cm = 8-16 km), a liquid 
body should have po — 2.5 xl0 15 -6.3 xlO 14 gram/cm 3 , according to the 
obtained formula for r c (6.13). 

Hence forth, we can calculate the mass of a non-rotating liquid col- 
lapsar, proceeding from the formulae M = | ira 3 po and a~r c ~ y/Z/ttpv 
(6.13). We obtain the following dependencies: 

M = — = 6.8xl0 27 a gram, (6.14) 

K 

A^tt 2.7xl0 41 ,„ . 

M = 3 = -=- gram. 6.15 

For example, once a collapsed liquid sphere should have the size of 
neutron stars and pulsars, a = (8-16) xlO 5 cm = 8-16 km, its mass would 
be M= (5.4-11) xlO 33 gram (i.e. 2.7-5.5 masses of the Sun). 

§6.2 The Universe as a huge liquid collapsar 

Here is another example: the Universe itself. Astronomers estimate the 
average density of substance in the Universe to be in the range of 10 -28 
to 10 -31 gram/cm 3 . Also, according to the observational estimates of 
the Hubble constant H = ^ = (2.3 ± 0.3) xlO -18 sec -1 , the radius of the 
Universe is a — 1.3 xlO 28 cm. At the upper boundary of the estimated 
density p 0 = 10 -28 gram/cm 3 , the collapse radius r c (6.10) meets the 
field of real numerical values. Respectively, we obtain, according to 
observational estimates, 



a = 


1.3 


xlO 28 cm 


Po = 


io- 


28 gram/cm 


M = 


9.2 


xlO 56 gram 


r g = 


1.4 


xlO 28 cm 


Tbr = 


1.3 


xlO 28 cm 


r c = 


1.5 


xlO 28 cm 



(6.16) 



This is a reason to suggest that the Universe can be considered as a 
sphere of perfect liquid, which is in the state of gravitational collapse. 
We will refer to it as the liquid model of the Universe. In this case, 
we should have r c = rbr = r g = a (6.12). Proceeding from this condition 
and the numerical value of the radius of the Universe, a = 1.3 xlO 28 cm, 
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M, gram 


Po, 


g/cm 3 


a 


cm 


r g , cm 


r br ! 


cm 


r c , cm 


Astron. 




















esteems 


9.2 xlO 56 




1Q -28 


1.3 


xlO 28 


1.4 xlO 28 


1.3 x 


10 28 


1.5 xlO 28 


Liquid 


8.8 xlO 55 




xlO" 31 




xlO 28 


1.3 xlO 28 




10 28 


1.3 xlO 28 


model 


9.6 


1.3 


1.3 x 



Table 5.1: The model of the observable Universe as a non-rotating liquid 
sphere in the state of gravitational collapse. The calculated parameters of 
the liquid model are compared to the observational esteems. 



obtained from the Hubble constant, we calculate the mass and density 
which should be attributed to the Universe in the framework of the 
present liquid model (according to a=r g = and M — |7ra 3 p 0 )- We 
obtain 

a = 1.3 xlO 28 cm 

po = 9.6 xl0~ 31 gram/cm 3 

M = 8.8 xlO 55 gram 

r g = 1.3 xlO 28 cm 

r br = 1.3 xlO 28 cm 

r c = 1.3 xlO 28 cm 

The obtained numerical values (6.17) are compared to the estimates 
of observational astronomy (6.16) in Table 5.1. Since these observa- 
tional estimates are known very approximately, we can conclude that 
the observable Universe is a huge collapsar, so all the world we observe, 
including us, is located within a huge black hole. 

In particular, this conclusion meets another made in 1965 by Kyril 
Stanyukovich [31]. He neither studied the geometric properties of a 
liquid sphere nor introduced a particular space metric. His analysis 
was based on the properties of elementary particles. Following this 
way, Stanyukovich obtained that the Hilbert radius of the Universe is 
the same as the observed event horizon: the observable Universe is a 
collapsar. Thus, despite employing another theoretical basis than our 
own, he arrived at the same conclusion. 

§6.3 Pressure and density inside liquid collapsars 

Pressure and density inside non-rotating liquid collapsars. . . The regular 
formula (2.130) we have obtained for the pressure p inside a sphere of 



(6.17) 



152 



Chapter 6 Black Holes 



perfect liquid, 



P = Poc 2 — r -- \ ° 2 , (6-18) 



under the collapse condition a — y/3/xpo takes the simplest form 

p = — poc 2 = const, (6.19) 



where p 0 = const by definition inside a sphere filled with perfect liquid. 
This formula is the equation of state of the liquid. This state is known as 
inflation: at positive density of the substance the pressure is negative, 
so the inner pressure of the substance tries to expand the body from 
within (despite a liquid body is incompressible). 

As is seen, the pressure is constant as well as the density. This means 
that the liquid substance, which fills a liquid collapsar, is in the state of 
inflation, and has the same pressure and density throughout the entire 
volume of the collapsar, from its center to the surface. 

§6.4 The inner forces of gravitation. The inner redshift 

The formula for the force of gravitation acting inside a non-rotating liq- 
uid collapsar can be found from the formula for the force acting inside 
a non-rotating liquid sphere, once the sphere is in the state of gravita- 
tional collapse (in this case, its physical radius is a = r g = y/3/xpo). 

Following this way, on the basis of the obtained formulae of the 
covariant component F\ (2.123, 2.125) and the contravariant component 
F 1 (2.124, 2.126) of the gravitational force, we obtain 

xp 0 c 2 r 1 
tx = 5 2 

1 3 

3 ~ a? ' 

Since r <a inside the sphere, F\ > 0. Therefore, this is a force of repul- 
sion. The force increases with distance r, from zero at the center of the 
liquid collapsar to its ultimate-high value on the surface. 

If the observable Universe is really a huge liquid collapsar (at least 
astronomical data evidence it, as was shown above), the repulsive radial 
force acting inside the collapsar may cause a frequency shift in travelling 
photons. To investigate this problem we consider the chr.inv. -equations 



c 2 r 



1 _ 11 



(6.20) 



(6.21) 
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of isotropic geodesies [18, 19] 

c l c k = 0 

(6.22) 



-F t c l + —D lk c l c k = 0 



^ + 2 W (^ + 4 1 )c l - W F +W A: i c" C ^0 

which are the equations of observable motion of a light-like (massless) 
particle which travels with the observable velocity of light c 1 (a photon 
whose frequency is oj) . The equations of isotropic geodesies result as the 
observable projections of the well-known generally covariant equations 
of isotropic geodesies (see [18,19] for details). 

In a rotation- free and non-deforming space (Aik — 0, Z),fc = 0), such 
as the space of a non-rotating liquid collapsar, the equations (6.22) take 
the form 

dui LJ 



2^C l =0 

CIT C Z 

d(uic l \ , 
V ; -wf + uA^cV =0 



(6.23) 



dr 

Let a photon travel only along the radial direction x 1 — r. Consider 
the chr.inv. -scalar geodesic equation (equation of energy) of the photon 
with the obtained formula for F\ (6.20) substituted. We also take into 
account that the photon's observable velocity is the observable velocity 
of light along the radial direction, c 1 = ^ . We obtain 



1 dto r dr 



oj dr a 2 — r 2 dr ' 
This equation solves as dlnu = — ^d\n\a 2 — r 2 \, or 



(6.24) 



d\nuj = dln 1 (6.25) 
y/a 2 — r 2 

Herefrom we obtain the function 

lo (r) = = : , Q = const. (6.26) 
V a 2 — r 2 

The integration constant Q is found from the obvious boundary condi- 
tion u) (r = 0) = ljo- It is Q = a 2 u 0 - Finally, we arrive at the solution 



154 



Chapter 6 Black Holes 



At distances (of the photon's travel) small to the physical radius of the 
collapsar (r <§; a) , this formula becomes 

W ~ W ° i 1 + h>) • (6 ' 28) 

This causes a square redshift (we also refer to it as a parabolic redshift, 
due to the parabolic square function) 

W ~ Wo 1 1>0 (6.29) 



w 0 



in the photon's spectrum: the force of repulsion F\, acting along the 
radial coordinate from the observer (in the observer's reference frame), 
deccelerates photons travelling from a distant object to him. At small 
distances, of the photon's travel (r <C a), the redshift is 



_, (6.30) 



2 a 

or, formulating this result through the Hubble constant H 

H 2 r 2 



2c 2 



(6.31) 



Thus, the observable parameters of the Universe manifest that it is 
a huge collapsar. These data match the calculations according to the 
theory of non-rotating liquid collapsars presented here. 



§6.5 The state of the collapsed liquid substance 

We now discuss the state of the substance that fills non-rotating liquid 
collapsars. As easy to see, once a liquid star is in the state of gravi- 
tational collapse (r g =a), the space metric of the star (6.7) takes the 
form 

ds2 = 7 f 1 - ^) c2dt2 - - r2 ( d ° 2 + sin ^ d( f) ■ ( 6 - 32 ) 

4 V a ) 1 — — 

This metric, under the particular condition a 2 = j > 0 (thus A > 0), has 
the same form as de Sitter's metric (1.5), 

ds 2 =(l-^-) c 2 dt 2 - - r 2 (dd 2 + sm 2 9 # 2 ) , (6.33) 
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which describes a spherical distribution of physical vacuum (the A-ficld 
in Einstein's field equations). 

This means that liquid collapsars consist of perfect liquid whose state 
is similar to the state of physical vacuum. The only difference is that 
the liquid that fills the collapsars possesses positive density, while the 
density of physical vacuum is negative with A > 0 (see §5.2 and §5.3 of 
our book [18] for details). Also, regular liquid collapsars have a small 
size and high density (in contrast to the Universe as a whole) . Therefore, 
the liquid that fills the regular (compact) collapsars is in a state similar 
to the state of high-density physical vacuum. 

What is physical vacuum, known also as the A-field? It is due to the 
general formulation of Einstein's field equations 

R a p ~ ^ 9af)R = -xT a/3 + \g a p (6.34) 

containing the A-term on the right-hand side. The right-hand side de- 
termines distributed matter which fills the space, while the left-hand 
side determines the space geometry which is Riemannian according to 
the formulation. Re-write the field equations in the form 

R a p ~ 7; 9apR = -xTap , (6.35) 

where the common energy-momentum tensor T a p = T a p + T a p charac- 
terizes both distributed substance and physical vacuum (A-field). 
The energy-momentum tensor of physical vacuum 

T a p = -— g a p (6.36) 
has the following physically observable projections 

p = = — — = const < 0 , (6.37) 
J 1 = = 0 , (6.38) 



U ik = c 2 f ik = - c 2 h ik = -pc 2 h ik . (6.39) 

which are calculated as well as the observable projections (1.84) of any 
energy-momentum tensor. 

The scalar chr.inv. -projection p = — = const implies that physical 
vacuum is homogeneously distributed in the space, i.e. is a homogeneous 
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medium. The vectorial chr.inv. -projection J 4 = 0 manifests that the 
physical vacuum is free of energy flow, i.e. is a non-emitting medium. 

Let us then find the equation of state of physical vacuum. According 
to the chronometrically invariant formalism [18,23], the chr.inv. -stress 
tensor U tk is expressed through the pressure inside a distributed medium 
as follows: 

Uik = Pohik - a ik = ph ik - (3 ik , (6.40) 

where po is the equilibrium pressure known due to the equation of state, 
p is the true pressure inside the medium, is the chr.inv. -viscous stress 
tensor, Pik — ctik — ^cthik is its anisotropic part which reveals itself in 
anisotropic deformations, while a = h a^k is the trace of the viscous 
stress-tensor aik- Since a spherically symmetric space is isotropic by 
definition, we have = 0 in the present case. Also, by the initial as- 
sumption, the medium is non- viscous (a ife = 0). Therefore, for physical 
vacuum, we have 

Uik = phik — ~pc 2 h ik ■ (6.41) 

Respectively, with the formula of the trace of the observable stress- 
tensor U — h lk Uik, we obtain the equation of state of physical vacuum 

p = -pc 2 , (6.42) 

which, with negative density p = — £ < 0, manifests the state of deflation 
(the inner pressure of the medium tries to compress the sphere). 

Hence forth, we deduce the covariant and contravariant components 
of the force of gravitation acting inside a vacuum (de Sitter) collapsar. 
Following the same way of deduction as that for the force acting inside 
a liquid collapsar (6.20, 6.21), we obtain 

F 1 = ^ L ^— , F^ = ^, (6.43) 

3 -I Ar 2 ' 3 J 

1 3 



while for the frequency shift of a photon we obtain 

A-- 
6 



W ° ' Co ( 1 + ^) • (6.44) 



ojq 1 Ar 



2 



> 0 (6.45) 



To understand the results, let us recall that we were able to trans- 
form the space metric of a collapsed liquid sphere (6.32) to de Sitter's 
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space metric (6.33) only by the particular condition a 2 = j > 0. Hence, 
we have assumed A > 0. With A > 0 we have obtained a negative density 
of physical vacuum p = — — <0 (6.37), the state of inflation p = — pc 2 
(6.42), the repulsing force fi > 0 (6.43), and the redshift (6.45). 

These are the same results as those obtained for a liquid collapsar, 
except for the negative density p——^<0 (and, hence, the positive 
pressure p — ~pc 2 >0 which gives the state of deflation). 

If we should assume a negative value of the A (i.e. A < 0), in order 
to obtain a positive density of physical vacuum, the collapsar's radius 
a would be imaginary, which is non-sense for the observed Universe. 

However, there is another way to remove this difficulty with respect 
to the theory. Consider Einstein's field equations (6.34) in a modified 
form wherein both the energy-momentum tensor of distributed sub- 
stance and the A-term are taken with the same sign, i.e. 

Rap ~ ^ g a pR = -xT a p - \g a p . (6.46) 

In this case, the energy-momentum tensor of physical vacuum is 

rr _ A 
J-a/3 — — g a /3 , 

while the physically observable projections of it are 

Too A 

p = = — = const > 0 , 

5oo x 

V5oo 

fjik = c 2 ftk = _ A c 2 h ik = _ p c 2 hlk ^g_ 50 ^ 

Given this case, physical vacuum (the A-field) is as well in the state 
of inflation (p = — pc 2 ), however its density is positive p= ^ > 0. Thus, 
the modified form (6.46) of Einstein's field equations removes the afore- 
mentioned contradiction between the theory of liquid collapsars and the 
observed positive density of substance in the Universe. 

Hence, physical vacuum (the A-field) is a homogeneous, non- viscous, 
non-emitting medium in the state of inflation. 

Concerning the deduced redshift formula (6.45), it depends only on 
the formulation of the force of repulsion which is deduced from goo of 
de Sitter's metric (6.33). Since we did not change the space metric, the 
redshift formula (6.45) remains unchanged as well. 



(6.47) 



(6.48) 
(6.49) 
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§6.6 Time flows in the opposite direction inside collapsars 

In a rotation-free space such as the space of a non-rotating liquid star, 
the observable time interval dr (1.30) has a simplified formulation which 
is dr^y/goo dt. Therefore, dr in the field of a non-rotating liquid star, 
according to goo of the metric (6.6), has the form 




xp 0 a 2 xpor 2 



dr = ±- 3 A l-^- A (6.51) 



Under the particular condition a = r g = ^/3/xpo characterizing the star 
in the state of gravitational collapse, this formula transforms to 



dr = T yi-^fdt. (6.52) 

We see that the sign of the observable time interval dr in a liquid 
star whose state is regular (out of collapse) is opposite to that in the 
star being a collapsar. In other words, the observable time in the field 
of regular stars flows in the opposite direction than the observable time 
inside a collapsar. Just one illustration: we regularly assume that ob- 
servable time flows from the past to the future. If so, the observable 
time inside collapsars flows from the future to the past. 

§6.7 The boundary conditions of a liquid collapsar 



With the condition a = r g = yj'i/ ' xpo characterizing liquid collapsars, 
the non-zero components of the Riemann-Christoffcl curvature tensor 
Ra/3-yS (2.113-2.116) obtained in §2.3 take the form 



#0101 = ^ttt = T~? = const ' ( 6 - 53 ) 
12 Aa z 



#1212 - - C 1212 = = , (6.54) 



xpo 


r 2 


r 2 


1 


3 


1 3 


~ a 2 




xpo 


r 2 sin 2 0 


r 2 


sin 2 (9 


3 


1 3 


~ a 2 


1-4 

a 2 



#1313 = -C1313 = — = — , (6.55) 



4 

i?2323 = - C 23 23 = ^ r 4 sin 2 6 = T — sin 2 9 . (6.56) 
3 a z 

Since i?oioi = — const and i?oioi > 0 at the positive density po > 0 of 
the liquid, the "inner" space of a liquid collapsar is a four-dimensional 
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positive constant curvature space. This is in contrast to our result of 
§2.3 where we showed that the space of a liquid sphere has a variable 
four- curvature which is negative. This means that: 

• The state of gravitational collapse is a "bridge" connecting the 
world of varying four-dimensional negative curvature (world of 
regular stars) and the world of four-dimensional positive constant 
curvature (inside those stars in the state of gravitational collapse) . 
Concerning the three-dimensional observable curvature of the space 
inside non-rotating liquid collapsars, we calculate Cn (2.104), C22 
(2.105), and the observable curvature scalar C = h lk Cik under the con- 
dition a = r g = yJZpHp® characterizing liquid collapsars. We obtain 

2x Po 1 2 1 

Ci1 - -——^z = Z ' (6 - 57) 

3 a 2 

C 33 2xp 0 r 2 2r 2 

C 22 = — — a = 7. = T ' 6 - 58 

sir 6 3 a 1 

6 

C = - 2^p 0 = — o = const < 0 • (6.59) 
a z 

It is a three-dimensional negative constant curvature space as well as 
the space of regular liquid stars (as a matter of fact that regular stars 
are out of the state of collapse). 

Hence forth, we express the force of gravitation acting in the "inner" 
space of a non-rotating liquid collapsar through the three-dimensional 
observable curvature of the "inner" space. From the respective formulae 
for Fi (6.20) and F 1 (6.21), we obtain 

2 2 

Fi = -^Cu, F 1 — C22 • (6.60) 
I Ir 

we see that both the three-dimensional observable curvature and the 
force of gravitation possess space breaking: 

Cu^-oo, Fi=^oo (6.61) 

by the boundary condition r = a on the surface of the collapsar. This 
result is, however, trivial. 

§6.8 Rotating liquid collapsars 

We now reveal rotating liquid collapsars. Let us rotate the space metric 
of a liquid collapsar (6.32) with an angular velocity uj around the line 



160 



Chapter 6 Black Holes 



orthogonal to the equatorial plane. Among the goi-th components of 
g a p, the non-zero component g 03 characterizes the rotation of space 
(while goi = go2 = 0). It formulates as 

2ujr 2 cos0 

.903 = , (6.62) 

so the linear velocity of space rotation Vi (1.45) is 

2ur 2 cos# „ ,„ „„. 

v 3 = — , vi = v 2 = 0 . (6.63) 

As a result, we obtain the space metric of a rotating liquid collapsar 
,2 1 ( r 2 \ 2 , 2 dr 2 

d,S 2 = - 1 - C 2 dt 2 



4 

2lot 2 cos 9 



1 - — 

a? 

cdt d(j) - r 2 (d9 2 + sm 2 6 d(j) 2 ) . (6.64) 



It is possible to prove that this space metric satisfies Einstein's field 
equations containing the energy-momentum tensor for perfect liquid 
(2.4). This means that, once we substitute the particular components 
of g a p taken from the metric (6.64) into the field equations, the left- 
hand side and right-hand side of the equations are the same: the field 
equations become identities and are thus satisfied. 

The general condition of gravitational collapse means that physical 
observable time stops (dr = 0). The definition of dr (1.30), 

dr = Jqlxidt H ^= dx l = JTi^dt Vidx % , (6.65) 

takes both g 0 o and goi into account. Therefore, with Vi ^ 0, the collapse 
condition is not dr = ,ygoo dt = 0 as that for non-rotating collapsars, but 
takes, in the present case under consideration, the complete form 

Vgoo - -4 = 0, (6.66) 
&■ 

where u 3 = ^ =lo. Substituting g 00 from the metric (6.64), w 3 (6.63), 
and u 3 — uj, we obtain the radius of the collapse surface of a rotating 
liquid collapsar 

r c = < a , (6.67) 

v 1 + — ? — 
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hence 

/ 2cj 2 a 2 cos6>\ 
r c ~a(l - )=a-Aa. (6.68) 

Assuming uj = 10 3 sec -1 and a = 10 6 cm for example, we obtain Ao~ 
22 cos 6, i.e. A a ~ 22 meters at the equator of the star, and Aa = 0 at 
the South Pole and North Pole. 

We see that the collapse surface matches the sphere's radius a only 
at the poles of the rotation (where the latitude 6 is ±-|, so we have 
cos 8 — 0). In other words, rotating liquid collapsars are not spheres but 
have an elliptic form, which is flattened on the equatorial plane (which 
is orthogonal to the axis of rotation) . 

Once the collapsar does not rotate (lj — 0) , its form is spherically 
symmetric (r c — a). Contrarily, at an ultimate relativistic speed of the 
rotation, the collapsar's elliptic form is highly flattened on the equatorial 
plane. In the ultimate case, where the collapsar rotates at a speed very 
close to the speed of light (w 2 a 2 — > c 2 ), its form is set up by the equation 

(6.69) 



Vl + 4cos6> ' 

The other parameters of rotating liquid collapsars we have obtained 
in the framework of the theory do not change the principal results ob- 
tained in §5.1 for non-rotating liquid collapsars. The only difference is 
that correction for the angular velocity of the collapsar's rotation to. We 
therefore omit these results from consideration. 

§6.9 Conclusion 

Finally, let us recall all the theoretical results of liquid collapsars that 
have been obtained above: 

1. The radial coordinate r c (6.10) by which a non-rotating liquid 
sphere of radius a meets gravitational collapse, is 




For regular stars, r c is in the range of imaginary numerical values. 
Therefore, regular stars ranging from super-giants to dwarfs and 
white dwarfs cannot collapse; 

By the requirement that the collapse radius r c should be real for 
real objects, the physical radius a of a collapsar should be 

a < 1.125 r g . 
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If its radius is a^l.l25r g , the non-rotating liquid body (star) 
cannot be in the state of gravitational collapse. 

The density of substance is the primary characteristic of non- 
rotating liquid collapsars. The physical radius a of such a collapsar 
is reciprocal to the square root of its density po (6.13) 

4.0xl0 13 

■ cm; 




Po 

4. The mass M of a non-rotating liquid collapsar is proportional to 
its physical radius a (6.14) 

M = — = 6.8xl0 27 a gram, 
x 

and is reciprocal to the square root of its density po (6.15) 

AVStt 2.7x10 41 
M = = — — gram; 

5. The observable Universe is completely located inside its collapse 
radius. In other words, it is a gravitational collapsar: all the stars 
and galaxies, including us, exist within a huge black hole. Its 
parameters, calculated according to the liquid model, are 

a = 1.3 xlO 28 cm, 

po = 9.6 xl0~ 31 gram/cm 3 , 

M = 8.8 xlO 55 gram; 

6. The liquid which fills the collapsars is in the state of inflation 

p = — poc 2 = const, 

i.e. at the positive density of the substance the pressure is nega- 
tive, so the inner pressure tries to expand the body from within 
(while the collapsar does not expand, because a liquid body is in- 
compressible). The pressure and density remain unchanged from 
the center of the collapsar to its surface; 

7. The gravitational inertial force acting inside a non-rotating liquid 
collapsar is a force of repulsion. It increases with distance, from 
zero at the center of the collapsar to its ultimate-high value on 
the surface; 
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8. The inner force of repulsion produces a square (parabolic) redshift 
in the photons travelling within the collapsar, to its center; 

9. The state of the liquid filling regular (compact) collapsars is sim- 
ilar to the state of high-density physical vacuum (high-density A- 
field), which is a homogeneous, non- viscous, non-emitting medium 
in the state of inflation; 

10. The observable time flows in different directions inside and out- 
side collapsars: once we assume that the observable time of our 
world flows from the past to the future, the observable time inside 
collapsars flows from the future to the past; 

11. The state of gravitational collapse is a "bridge" connecting the 
world of varying four-dimensional negative curvature (world of 
regular stars) and the world of four-dimensional positive constant 
curvature inside gravitational collapsars (black holes); 

12. Rotating liquid collapsars are not spheres but possess an elliptic 
form which is flattened on the equatorial plane. The radius r c 
of a rotating liquid collapsar is formulated through the sphere's 
radius a, the latitude 0, and the angular velocity of the respective 
rotation oj as 



a 




Notations 



Ordinary differential of a vector: 

BA a 

dA a = ^—dx a . 
dx° 

Absolute differential of a contravariant vector: 

DA a = V,3 A a dx 13 = dA a + T^A^dx 13 . 

Absolute differential of a covariant vector: 

DA a = Vp A a dx? - dA a - T^A^dx?. 

Absolute derivative of a contravariant vector: 

v« A a - BA " - + r a a» 
^ A ~ dxP - 8x0 +l ^ A ■ 

Absolute derivative of a covariant vector: 

v a a - DAa - -T* A 

V/3 a ~ dx& ~ dx? 

Absolute derivative of a 2nd rank contravariant tensor: 



Absolute derivative of a 2nd rank covariant tensor: 

BF 

Absolute divergence of a vector: 

8A a 

V7 A a — - i- V a A a 

VaA ~ dx° + aa ' 

Chr.inv. -divergence of a chr.inv. -vector: 

dx l dx l dx % ^ 

Physical chr.inv. -divergence: 
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D'Alembert's general covariant operator: 

□ = g a(i V a V 0 . 
Laplace's ordinary operator: 

A = -g lk V l V k . 

Chr.inv. -Laplace operator: 

*A = h tk *V l *V k . 

Chr.inv. -derivative with respect to the time coordinate and that with 
respect to the spatial coordinates: 

*d _ 1 d *d _ d 1*8 

dt y/goo dt ' dx 1 dx 1 c 2 1 dt ' 

The square of the physically observable velocity: 

v 2 = VjV 4 = h ik v\ k . 
The linear velocity of the space rotation: 

Vi = - c — — , v = - eg" ^/g^j , Vl = h ik v . 



9oo 

The square of «j. This is the proof: because of g a ag ati = fff , then under 
a — (3 — 0 we have goag a0 = <>o = 1, hence v 2 = v k v k = c 2 (l — .goo5°°), i-e.: 

v 2 = h lk v l v k . 

The determinants of the metric tensors g a p and h a p are connected as: 

V - 9 = \/h\/9oo- 

Derivative with respect to the physically observable time: 

d__*d_ fc *d 
dr dt dx k 

The 1st derivative with respect to the space-time interval: 

d 1 d 



ds f-, V* dr 
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The 2nd derivative with respect to the space-time interval: 
i 1 d 2 If- V+v^ + i^v^v-V 



i 2 dr 2 ( c 2_ v 2^ 2 \ dr 2 dx m J dr' 

The chr.inv. -metric tensor: 

h lk = -g lk + \ v iVk , h tk = - g ik , h\ = 5 k . 
cr 

Zelmanov's relations between the Christoffel regular symbols and the 
chr.inv. -characteristics of the space of reference: 

Di + Ai = -^(r 0k - 90kT ™ 



goo V 5oo 

r 2 r k 

Ja^feppm uiQuks\m j^k 00 
9 9 1 a B - h h fl , s . h 

5oo 

Zelmanov's 1st identity and 2nd identity: 



0. 



*8Aj k 1 (*dF k *dFj\ 
dt 2\ dx l dx k J 

} hm + + + \ {FiA-km + F k A mi + F m A ik ) = 0 . 



dx l dx k dx m 2 
Derivative from v 2 with respect to the physically observable time: 

) ) =2AfcW +^vv v +2v fe — . 

The completely antisymmetric chr.inv. -tensor: 

-ifem — jpOikm e 0tkm -Ebifcm rr 
£ — V900 & — J=- j ^ifcm — — eoifcmVft. 

V'l V5oo 
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produced by the thermonuclear fusion is that energy which the stars shine: the tiny 
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then transferred to the physical surface of the star due to the thermal conductiv- 
ity. A new classification of stars is introduced according to the space breaking of 
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